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McCON  -  A  GENERAL  CONTOURING  PROGRAM  FOR  PERSONAL  COMPUTERS 


Introduction 


1.  The  need  to  prepare  contour  maps  and  profiles  along  cross-sections 
often  arises  during  the  course  of  engineering  and  geologic  investigations. 
Although  topographic  applications  are  the  most  widely  known,  contour  maps  and 
plots  are  prevalent  for  stress  analysis,  groundwater  analysis,  and  for  almost 
any  type  of  data  for  which  trends  are  to  be  examined,  estimated,  or  predicted. 
Over  the  past  twenty  years,  numerous  computer  programs  have  been  developed  on 
numerous  electronic  platforms  to  aid  in  the  preparation  of  contour  maps.  The 
user  of  any  of  the  various  programs  (including  the  one  described  herein)  is 
cautioned  to  be  aware  of  the  application  for  which  the  program  was  written  and 
to  be  aware  of  the  methods  employed  to  estimate  the  contourable  values  in  the 
vicinity  of  the  given  data  sets.  The  recent  advent  of  high  quality  graphics 
capabilities  on  personal  computers  and  workstations  has  resulted  in  the 
availability  of  a  large  number  of  contouring  programs  at  reasonable  cost  from 
a  variety  of  suppliers.  The  program  McCON  Was  developed  such  that  the  "raw" 
data  could  be  obtained  either  from  self -prepared  data  files  or  from  data  files 
commonly  produced  by  other  computerized  applications.  It  was  also  developed 
to  permit  the  inclusion  of  zones  across  which  the  data  is  discontinuous  and  to 
permit  the  preparation  of  profiles  along  cross-sections  at  any  orientation. 

2.  McCON  is  a  general  purpose  contouring  program  which  executes  rapidly 
on  IBM  compatible  personal  computers.  The  minimal  data  needed  to  execute  the 
program  is  a  free  field  data  file  of  the  location  coordinates  (x,y)  and  the 
contour  values  (z)  of  the  data  points.  As  configured,  the  program  will  accept 
up  to  999  data  (x,y,z)  triplets.  The  data  points  do  not  need  to  follow  any 
regular  spacing.  Negative  values  of  x,y,  or  z  are  permitted  and  the  (x,y) 
coordinates  may  have  an  arbitrary  range  of  values.  The  coordinate  ranges  are 
internally  examined  and  scaled  to  yield  a  contour  plot  which  fills  the  video 
screen. 

3.  Various  Options  are  provided  by  McCON: 

a.  A  noncontoured  subregion  (i.e.,  an  area  within  which  the  user 
does  not  wish  contour  lines  to  be  drawn)  may  be  specified  provided 
that  this  region  is  of  a  convex  shape. 

b.  Contouring  may  be  restricted  to  be  entirely  contained  within  a 
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given  convex  subregion. 

c.  Contouring  of  areas  having  physical  boundaries  that  are  not 
convex  may  be  prepared  by  subdividing  the  overall  area  into  convex 
"composite"  regions. 

d.  Internal  boundaries  across  which  contours  are  discontinuous 
(e.g.,  geologic  fault  lines,  seepage  barriers,  etc.)  may  be  includ¬ 
ed  . 

e.  "Templates"  (sets  of  isolated  or  connected  lines)  to  be  drawn 
onto  the  contour  plot  to  enhance  the  artwork  of  the  drawing  are 
included. 

f.  Profiles,  or  cross-sections,  along  any  plan  orientation  may  be 
created  and  separately  plotted. 

4.  McCON  requires  590K  of  RAM  memory  for  the  executable  program.  The 
source  larguage  is  Microsoft  (R)  FORTRAN.  MicroGCS  (the  graphics  compatibili¬ 
ty  system  developed  by  the  US  Army  Corps  of  Engineers)  is  used  to  support  the 
graphics.  All  calls  to  subroutines  beginning  with  the  letter  "U"  in  the 
accompanying  source  listing  (Appendix  B)  are  MicroGCS  subroutines.  The 
graphical  output  may  also  (by  option)  be  sent  to  a  variety  of  pen- type  flatbed 
plotters  via  routines  provided  by  MicroGCS. 

5 .  The  program  generates  nonintersecting  triangles  which  connect  each 
and  every  data  (node)  point.  Triangle  generation  ceases  when  all  of  the  node 
points  are  connected  to  at  least  one  triangle  and  the  resulting  mesh  of 
triangles  encompasses  all  of  the  node  points  in  a  convex  fashion  (i.e.,  the 
outer  edges  of  the  triangle  mesh  form  a  convex  shape) .  The  resulting  mesh 
will  contain  no  areas  that  are  not  included  within  a  triangle  (i.e.,  the  mesh 
will  contain  no  "holes").  A  scheme  was  developed  to  create  the  mesh  in  a  one- 
pass  "spawning"  process;  therefore,  as  each  triangle  is  created,  sufficient 
information  exists  to  immediately  draw  the  contours  within  the  new  triangle. 
Because  of  this  one-pass  triangle  generation  process,  the  execution  time  is 
quite  rapid.  Typically,  a  set  of  100  nodes  (on  an  82386  machine  with  math  co¬ 
processor  and  EGA  card)  will  require  approximately  10  sec  to  generate  the 
triangle  mesh;  a  set  of  400  nodes,  56  sec;  and  900  nodes,  165  sec.  The  time 
devoted  to  contour  line  drawing  is  dependent  upon  the  resulting  number  of 
contour  lines  to  be  drawn.  A  relatively  heavy  set  of  contour  lines  may 
require  30-45  sec  of  execution  time  for  the  contours  to  be  drawn  on  the 
screen;  however,  the  execution  time  is  not  significantly  related  to  the 


complexity  of  the  contours  nor  to  the  spatial  distribution  of  the  nodes. 
Considerably  more  time  for  drawing  is  required  if  the  output  is  sent  to  a 
plotter . 

6.  The  contour  lines  are  drawn  as  a  series  of  connecting  straight  line 
segments  and  circular  segments.  This  combination  gives  an  aesthetically 
pleasing  appearance  to  the  resulting  contour  plot. 

Uniqueness  of  Contours 

7 .  The  process  of  defining  the  locus  of  points  of  equal  contour  value 
(the  contour  line)  is  always  dependent  upon  the  manner  in  which  the  desired 
contour  values  are  established  (i.e.,  interpolated)  between  the  given  data 
points.  Figure  1  shows  two  differing  manual  interpretations  as  regards 
contouring  the  elevations  given  at  4  points.  Figure  la  shows  the  result 
obtained  by  interpolating  the  elevation  data  along  the  right  descending  dashed 
line;  whereas  Figure  lb  is  the  result  of  interpolation  along  a  right  ascending 
line.  It  is  obvious  that  the  two  schemes  give  radically  different  interpreta¬ 
tions  . 

8.  However,  neither  interpretation  for  establishing  triangles  is 
inherently  superior.  In  this  case,  more  data  are  required  to  properly  estab¬ 
lish  the  character  (or  trend)  of  the  contours.  If  it  were  established  that 
the  elevation  at  the  mid-point  was  equal  to  40,  then  the  interpretation  shown 
in  Figure  la  would  be  proper.  Figure  lb  would  be  applicable  if  the  mid-point 
elevation  was  equal  to  65.  Resorting  to  polynomial,  power,  or  spline  fits 
does  not  necessarily  improve  the  interpretation;  each  technique  would  only 
lead  to  a  different  estimate.  All  contouring  techniques  rely  on  estimating 
the  behavior  of  the  contourable  value  within  the  neighborhood  of  a  few 
(typically  3  or  4)  given  field  values.  Within  that  neighborhood,  many 
interpretive  schemes  are  possible.  However  each  scheme  implies  a  knowledge  of 
the  character  of  the  "landform."  The  interpolation  scheme  used  for  McCON  is 
quite  similar  to  the  manual  interpolation  methods  discussed  above  (i.e.,  a 
linear  interpolation  of  the  contour  value  gradient  from  node  to  node  along 
triangle  edges) .  The  way  and  order  in  which  the  triangles  are  generated  may 
influence  the  contour  plot.  In  any  case,  one  interpretation  is  as  proper  as 
any  other  in  light  of  the  lack  of  other  data.  There  is  no  good  reason  to 
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a.  Result  obtained  by  inter-  b.  Result  of  interpolation 

polating  elevation  data  along  a  different  line 

along  dashed  lines 


Figure  1.  Converting  elevations  to  contours 


state  that  an  alternate  interpretation  is  better  unless  additional  information 
is  known. 


Mesh  Generation 

9.  The  scheme  to  generate  the  triangle  mesh  is  summarized  as  follows: 

a.  Create  a  single  "seed"  triangle  from  any  three  suitable  node 
points  (i.e.,  the  seed  triangle  must  not  contain  any  other  node 
points  within  its  boundaries) . 

b.  Spawn  additional  triangles  from  each  edge  of  the  seed  triangle 
such  that  the  spawning  and  the  spawned  triangle  share  the  common 
edge.  The  criteria  for  determining  which  node  point  to  use  to 
create  the  new  triangle  are  given  below. 

c.  In  the  order  in  which  they  were  created,  use  the  previously 
spawned  triangles  to  spawn  new  triangles.  After  each  triangle  is 
spawned,  all  of  the  previously  spawned  triangles  are  examined  to  see 
if  the  new  triangle  is  sharing  an  edge  other  than  the  edge  of  the 
parent  triangle.  Each  newly  spawned  triangle  will,  at  most,  form 
two  new  edges  (the  third  edge  is  already  shared  with  its  parent) 
from  which  to  subsequently  spawn  additional  triangles.  The  new 
triangle  may  furnish  only  one  (or  no)  new  edge  for  further  spawning 
since  the  new  triangle  may  also  share  edges  with  previously 
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generated  triangles.  As  triangles  are  spawned,  a  spawning  sequence 
number  is  assigned.  Triangles  may  be  spawned  only  from  edges  that 
are  not  already  shared  by  two  triangles. 

d.  Step  c  is  repeated  until  after  some  spawn  (which  will  be  the 
last  spawn)  it  is  detected  that  the  order-of-sequence  number  of  the 
spawning  triangle  is  equal  to  the  total  number  of  generated  trian¬ 
gles.  When  this  occurs,  the  mesh  of  generated  triangles  will 
connect  and  encompass  all  of  the  node  points.  The  result  will  be  a 
triangle  mesh  with  its  outer  edges  forming  a  convex  shape.  Spawning 
will  continue  as  long  as  the  order-of-sequence  number  of  the  spawn¬ 
ing  triangle  is  less  than  the  number  of  generated  triangles. 

10.  Figure  2  shows  the  order  in  which  the  10  triangles  connecting  the  9 
node  points  were  spawned.  Triangle  1  (the  seed  triangle)  spawned  triangles  2, 
3 ,  and  4 .  Triangle  2  spawned  triangles  5  and  6 .  Triangle  3  spawned  triangles 
7  and  8.  (Triangle  7  also  shares  another  of  its  edges  with  the  by  then 
existing  triangle  4.)  Triangle  4  only  spawned  triangle  9  since  triangle  7  had 
already  been  spawned  by  triangle  3.  Triangle  5  spawned  triangle  10.  Trian¬ 
gles  6,  7,  8,  9,  and  10  spawned  nothing.  The  outer  edge  of  the  mesh  (de¬ 
scribed  by  node  points  2,  6,  8,  7,  9,  and  4)  is  convex. 

11.  All  triangles  (except  the  seed  triangle)  are  spawned  by  an  edge  of 
a  parent  triangle.  The  parent  triangle  provides  two  vertices  for  the  next 
formed  triangle.  There  are  three  criteria  to  decide  which  node  to  select  as 
the  third  vertex  of  the  new  triangle. 

CRITERION  A.  Every  node  point  which  lies  on  the  side  of  the  spawning 
edge  which  is  opposite  to  the  side  containing  the  remaining  vertex  of 
the  spawning  triangle  becomes  a  candidate  node. 

This  criterion  prevents  new  triangles  from  intersecting  (overlapping)  their 
parent.  As  triangles  are  formed,  the  node  point  numbers  of  the  vertices  of 
the  trian  le  are  stored  in  an  ordered  array.  This  array  is  ordered  such  that 
the  vertices  are  listed  in  a  counterclockwise  fashion  about  the  triangle.  For 
example,  referring  to  Figure  2,  the  vertex  .tray  for  triangle  1  must  be 

(3.1.5)  (or  (1,5,3)  or  (5,3,1))  since  these  numbers  go  counterclockwise  about 
the  triangle.  Similarly,  triangle  8  is  described  as  (5,4,9)  or  (9,5,4)  or 

(4.9.5) .  The  formula  for  the  area,  A,  of  a  triangle  in  terms  of  its  vertex 
coordinates  is: 
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(1) 


A  =  [  Xi  (y2-y3)  +  x2  (y3-yi)  +  x3  (yx-yz)  ]  /  2 
or  symbolically, 

A  =  Area(l ,2,3) 

where  xa  and  yx  are  the  coordinates  of  vertex  1,  etc.  This  formula  yields  a 
positive  value  if  the  points  1,2,3  are  ordered  in  a  counterclockwise  sense 
about  the  triangle  and  a  negative  value  if  ordered  clockwise.  Therefore, 
knowing  the  fashion  in  which  the  spawning  edge  is  ordered  (say  3,5  for 
triangle  1)  only  those  nodes  which  yield  a  positive  (counterclockwise)  area 
(i.e.  Area  (3,5,6),  Area  (3,5,8),  or  Area  (3,5,7))  are  on  the  proper  side  of 
spawning  edge  3,5.  All  other  nodes  lie  on  the  same  side  of  the  spawning  line 
as  the  spawning  triangles  remaining  vertex  (node  1)  and  the  application  of  the 
area  formula  gives  negative  or  zero  areas  for  these  nodes.  Equation  1  is 
applied  often  during  the  course  of  mesh  generation.  A  mesh  formed  by  250 
nodes  typically  will  require  over  60,000  area  calculations;  999  nodes,  over 
225,000  calculations. 

12.  For  each  node  point,  if  any,  that  passes  Criterion  A,  the  following 
criterion  must  be  met. 

CRITERION  B.  The  angle  at  the  vertex  of  the  "trial"  new  triangle  which 
is  opposite  the  spawning  edge  must  be  the  largest  angle  formed  by  the 
set  of  candidate  nodes  from  Criterion  A.  (Provided  that  the  candidate 
node  also  passes  Criterion  C  below.) 

This  criterion  will  ensure  that  no  other  node  points  may  lie  within  the  new 
triangle.  Referring  to  Figure  2,  presume  that  triangle  1  is  attempting  to 
spawn  a  triangle  from  edge  3,5.  In  order,  the  candidate  nodes  would  be  6,7, 
and  8.  (Node  9  is  rejected  since  Area  (3,5,9)  is  zero.)  Node  point  8 
supplies  the  largest  vertex  angle.  If  it  is  imagined  that  some  additional 
node  point  lies  within  triangle  2,  it  is  obvious  (for  either  acute  or  obtuse 
triangles)  that  the  imaginary  point  subtends  a  larger  angle  than  the  angle  3- 
8-5  and,  as  such,  would  be  a  better  choice  to  spawn  a  triangle. 

CRITERION  C.  Reject  all  candidate  nodes  that  yield  a  negative  (clock¬ 
wise)  area  (per  the  formula  given  in  Criterion  A)  when  associated  with 
any  line  previously  emanating  from  the  first  end  of  the  spawning  edge 
and  reject  those  candidate  nodes  which  yield  a  positive  area  when 


8 


associated  with  any  line  emanating  from  the  second  end  of  the  spawning 
edge . 

The  area  formula  (Equation  1)  is  applied  such  that  Xj.y!  is  the  coordinate  of 
the  spawning  edge  endpoint  under  consideration,  x2,y2  is  the  coordinate  of  the 
node  point  at  the  end  of  the  line(s)  emanating  from  the  end  of  the  spawning 
edge,  and  x3,y3  is  the  coordinate  of  the  candidate  node  point.  For  the 
purposes  of  applying  this  criterion,  it  is  only  necessary  to  perform  checks  on 
lines  which  emanate  above  the  horizon  of  the  spawning  edge.  The  purpose  of 
this  criterion  is  to  ensure  that  no  triangles  intersect  each  other. 

13.  To  illustrate  Criterion  C,  examine  the  conditions  shown  in  Fig¬ 
ure  3.  Suppose  that  the  triangles  shown  by  the  solid  lines  were  created  in 
some  proper  fashion  and  that  the  calculations  have  proceeded  to  the  point 


where  side  1-2  of  triangle  P  is  attempting  to  spawn.  According  to  Criterion 
A,  nodes  a,b,c,d,e,  and  g  are  all  possible  candidates  because  they  lie  above 
the  1-2  "horizon"  (i.e.,  Area  (2,1,?)  is  positive  for  all  labeled  nodes  except 
f ) .  By  inspection,  it  is  seen  that  only  choices  a,e,  or  g  would  result  in  a 
non- intersecting  new  triangle.  However,  applying  Criterion  B,  the  maximum 
apex  angle  would  occur  for  node  c  (and  the  next  best  angle  at  node  b) . 
Criterion  C  will  reject  the  candidacy  of  nodes  b,c,  and  d  because  they  form 
clockwise  ordered  triangles  with  the  line  emanating  from  node  1  to  node  a. 

That  is,  Area  (l,a,b),  Area  (l,a,c),  and  Area  (l,a,d)  are  all  ordered 
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clockwise.  Nodes  g  and  e  are  not  rejected  since  Area  (l,a,g)  and  Area  (l,a,e) 
are  ordered  counterclockwise.  Area  (l,a,a)  yields  an  area  of  zero  (which,  for 
the  purposes  of  Criterion  C,  is  non-negative)  and  is  retained  as  a  candidate 
node.  For  the  example  shown,  the  line  emanating  from  the  opposite  end  of  the 
spawning  edge,  line  2-g  does  not  cause  any  rejections  because  all  node  points 
under  consideration  give  clockwise  triangles  in  accordance  with  the  second 
half  of  Criterion  C.  Therefore,  node  a  (which  yields  the  maximum  vertex 
angle)  will  become  the  correct  choice  to  spawn  a  new  triangle.  During  the 
course  of  the  calculations  a  linked  list  of  all  of  the  node  points  which  are 
connected  (via  an  edge  of  a  triangle)  to  a  given  node  is  maintained.  This 
list  provides  the  means  for  applying  Criterion  C. 

14.  The  node  point  which  passes  all  of  the  criteria  (and  there  is  at 
most  one  node  that  can  pass)  is  used  to  form  a  new  triangle.  In  the  event 
that  no  node  passes  the  criteria,  no  triangle  is  spawned  from  the  edge  under 
consideration  and  the  calculations  proceed  to  the  next  edge  to  be  considered. 
Since  Criterion  B  guarantees  that  no  nodes  may  be  trapped  within  a  triangle, 
then  each  node  must  be  used  to  form  at  least  one  triangle  simply  because  all 
nodes  must  lie  outside  at  least  one  of  the  edges  of  every  triangle  (except  the 
nodes  which  form  that  given  triangle).  Therefore,  the  spawning  process  cannot 
stop  until  all  nodes  have  been  used.  The  completed  mesh  of  triangles  must 
form  a  unique  convex  shape  about  all  of  the  nodes  since  this  is  the  only  shape 
for  which  no  other  nodes  will  be  available  for  spawning.  Voids,  or  interior 
areas  not  included  within  a  triangle,  are  not  possible  since  the  void  would  be 
adjacent  to  two  or  more  triangles  with  unshared  edges.  During  the  spawning 
process  these  edges  would  have  had  an  opportunity  to  spawn.  The  criteria  for 
spawning  require  that  the  spawning  edge  always  spawns  unless  no  node  exists 
outside  of  that  edge  or  if  the  spawn  would  result  in  an  intersecting  triangle. 

Data  Input  Enhancements 

15.  McCON  provides  several  means  to  "massage"  the  input  data.  (Appen¬ 
dix  A  gives  a  more  complete  discussion  of  the  capabilities  for  data  enhance¬ 
ment.  )  In  general,  the  input  data  file  may  consist  of  up  to  20  columns  of 
data  values.  The  user  identifies  the  columns  which  contain  the  x  and  y 
coordinates  and  the  z  (contourable)  value.  The  z  value  may  be  the  sum  or 
difference  of  any  of  the  columns  of  the  data  file.  Therefore,  a  data  file  may 
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contain  sufficient  information  to  yield  many  contour  plots. 

16.  The  data  input  activity  of  McCON  permits  the  user  to  identify  the 
elements  of  the  original  input  data  file  and  to  enhance  the  data  in  the 
following  ways: 

a.  "As  is . "  Contouring  of  the  original  data  points  will  take 
place  without  any  data  smoothing  and  without  any  "border"  data 
point  creation. 

b.  "Original  points  with  boundary."  Additional  data  points  are 
created  which  form  a  rectangular  boundary  around  the  original 
points.  The  contourable  values  (elevations)  of  the  boundary 
points  are  computed  by  Equation  2  (given  below) . 

c.  "Grid  points  with  or  without  apron."  A  grid  which  overlays 
the  original  data  points  is  created.  The  contourable  values 
(elevations)  of  the  grid  points  are  computed  by  Equation  2  (given 
below) .  The  original  data  points  are  optionally  superimposed  on 
the  grid  points. 

d.  "Triple  triangle  refinement."  Following  the  generation  of  a 
mesh  of  triangles,  an  option  will  subdivide  each  "old"  mesh  trian¬ 
gle  into  three  sub- triangles .  The  contourable  value  (elevation)  at 
the  added  point,  (mid-point  of  the  old  triangle),  may  be  computed 
by  Equation  2  or  by  linear  interpolation. 

The  last  three  data  enhancements  listed  above  will  generally  result  in  a 
smoothing  of  the  contour  plots;  however,  increased  contour  lines  "waviness" 
may  result.  The  contourable  values  (elevations)  at  "created"  points  are 
computed  by  the  following  equation  for  inverse  power  distance  averaging: 


Zi  =  £  for  all  dt  j  *  o 

l  i 

Zi  =  Zj  for  any  dt  j  =  0 

where  Zi  =  the  elevation  at  the  created  point(s)  i, 

Zj  =  the  elevations  at  the  given  original  points,  j, 
dj  j  =  the  distance  between  point  i  and  j , 

Q  =  an  exponent  to  be  chosen  (usually  2) . 


(2) 


Equation  2  will  give  a  distance  weighted  average  for  the  elevations  to  be 
assigned.  Gridding,  which  always  requires  inverse  power  distance  averaging, 
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will  generally  result  in  contour  lines  that  are  more  circular  in  shape  than 
"as  is"  contours  and  which  often  isolate  local  highs  and  lows  with  closed 
contours.  An  excessive  number  of  grid  points  may  result  in  undue  contour  line 
waviness.  If  no  gridding  is  used,  the  contours  generally  exhibit  more  linear 
trends.  "Triple  triangle  refinement"  provides  the  only  method  for  input  data 
enhancement  for  "as  is"  contouring.  Again,  inverse  power  distance  averaging 
for  "triple  triangles"  may  result  in  many  closed  contours  (isolation  of  local 
highs  and  lows).  Linear  averaging  for  "triple  triangles"  is  generally 
recommended . 

17.  The  results  from  various  data  input  manipulations  are  shown  in 
Figure  4.  The  location  and  contourable  value  (elevation)  of  32  nodes  is  shown 
in  Figure  4a.  The  "as  is"  contour  plot  is  shown  in  Figure  4b.  The  results  of 
"gridding"  the  data  and  not  including  the  "original"  data  points  are  shown  in 
Figure  4c.  Notice  that  the  contours  are  more  circular.  When  the  original 
data  points  for  "gridded"  data  are  honored,  the  results  are  shown  in  Fig¬ 
ure  4d.  "Triple  triangle"  plots  are  shown  in  Figure  4e  and  in  Figure  4f. 

Contour  Drawing 

18.  Since  the  spawning  process  outlined  above  is  a  single-pass  process, 
the  contour  lines  can  be  drawn  as  the  triangles  are  spawned.  Figure  5  shows  a 
triangle  being  spawned.  The  numbers  at  the  nodes  represent  contourable  value 
(or  elevation).  Presume  that  it  is  required  to  draw  the  10  and  20  contours. 
Regardless  of  the  particular  values  of  the  elevations  at  the  nodes  surrounding 
any  triangle,  any  given  elevation  for  which  a  line  is  to  be  drawn  either  does 
not  intersect  the  triangle  at  all  or  it  intersects  two  sides  of  the  triangle. 
(The  special  case  where  two  nodes  of  the  triangle  have  equal  values  results  in 
the  contour  line  being  that  triangle’s  edge.)  The  contour  lines  are  drawn  as 
follows  as  the  triangles  are  spawned. 

a.  In  the  order  of  all  contour  lines  to  be  drawn,  determine 
whether  the  line  will  pass  between  the  ends  of  the  shared  edge  (as 
do  the  10  and  20  contour  lines).  If  not,  do  not  draw  lines;  if  so, 
go  to  the  next  step. 

b.  Compute  (by  linear  interpolation)  the  locations  where  the 
contour  line  intersects  the  shared  edge  and  also  where  it  inter¬ 
sects  the  second  edge  of  both  the  spawning  and  the  spawned 
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Figure  U .  Contours  resulting  from  various  input  data  enhancements 
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triangles.  Mathematically  construct  two  lines  connecting  these 
points.  Then  compute  the  location  of  the  mid-point  (designated  by 
the  "+"  symbols  in  Figure  5)  of  each  line  in  each  triangle.  If 
the  contour  line  is  to  be  drawn  with  "sharp"  corners,  the  two  lines 
going  from  the  mid-points  to  the  shared  edge  are  the  contour  lines 
to  be  drawn  on  the  screen  (or  plotter) .  The  other  "half"  of  the 
contour  line  will  not  be  drawn  until  a  triangle  is  spawned  sharing 
the  now  "unshared"  edge .  The  reason  for  leaving  half  of  the 
triangle  unfinished  is  to  accommodate  the  introduction  of  circular 
contour  line  segments  which  will  smooth  the  contour  lines.  In  the 
event  that  an  unfinished  side  does  not  share  an  edge,  (as  would  be 
the  case  for  all  triangles  adjacent  to  the  outer  convex  edge  of  the 
mesh)  the  program  will,  after  all  spawning  is  complete,  seek  all 
triangle  edges  that  are  unshared  and  finish  drawing  the  contour 
lines  out  to  the  mesh  boundary. 

19.  Circular  segments  (usually  desired)  for  the  contour  lines  are  drawn 
in  the  vicinity  of  the  shared  edge.  Figure  6  shows  a  detail  of  the  10  contour 
line.  The  points  o,  a,  and  b  are  the  same  points  as  shown  in  Figure  5.  All 
circular  segments  are  constructed  by  calculating  the  radius  and  center 
coordinates  of  the  circle  S  that  will  be  tangent  to  the  lines  oa  and  ob  at  a 
distance  of  a  from  the  point  o.  The  program  initializes  £  to  a  value  of  25 
screen  units  (the  width  of  the  entire  contour  plotting  area  is  1040  screen 
units)  and  does  not  permit  the  value  of  a  to  exceed  25  screen  units.  The 
variable  a  is  reduced  to  d  (where  d  is  the  length  of  the  shorter  of  oa  or  ob) 
if  d  is  less  than  25  screen  units.  The  distance  a  is  also  reduced  sufficient¬ 
ly  to  assure  that  the  circle  to  be  drawn  crosses  the  shared  edge.  Figure  7 
illustrates  a  situation  where  the  contour  line  would,  without  a  reduction  of 
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£,  not  cross  the  shared  line.  It  is  necessary  to  ensure  that  all  contour 
lines  are  drawn  entirely  within  the  two  triangles  under  consideration; 
otherwise,  problems  with  crossing  contour  lines  could  arise. 


Noncontoured  Region 

20.  At  times,  it  is  desirable  NOT  to  draw  contours  within  some  finite 
area  on  the  contour  plot.  For  example,  the  region  beneath  a  building  is 
generally  not  contoured  in  conjunction  with  the  surrounding  area.  McCON 
allows  the  user  to  exclude  from  contouring  a  single  internal  convex  area.  As 
previously  stated,  the  program  generates  triangles  that  always  completely 
encompass  the  range  of  data  points  without  "holes."  The  triangle  generation 
ceases  only  when  the  convex  set  of  triangles  connecting  the  exterior  nodes  is 
accomplished.  The  ability  to  create  a  noncontoured  region  requires  a  two-step 
spawning  process.  The  user  must  specify  which  nodes  form  the  outline  of  the 
noncontoured  region.  (Those  nodes  should  specify  an  internal  convex  region. 

If  the  specified  nodes  do  not  form  a  convex  set,  the  program  will  create  a 
convex  region  defined  by  those  points  that  may  not  conform  to  the  user's 
wishes.)  The  first  stage  of  the  triangle  generation  process  uses  only  those 
specified  nodes  to  create  a  mesh.  The  result  of  this  stage  will  be  a  convex 
mesh  that  has  an  external  boundary  coincident  with  the  convex  set  formed  by 
the  specified  nodes.  A  list  of  the  edges  of  the  triangles  which  form  the 
convex  outer  shell  of  the  noncontoured  region  is  maintained  (i.e.,  a  list  of 
edges  that  are  not  shared)  for  later  use.  No  contour  line  drawing  is  done 
during  this  stage.  Stage  two  is  accomplished  by  allowing  all  of  the  remaining 
nodes  to  be  used  for  spawning  additional  triangles.  The  spawning  is  renewed 
by  beginning  with  the  first  triangle  (the  seed  triangle)  created  at  stage  one 


and  stepping  through  each  triangle  in  order  of  its  creation.  No  new  triangles 
will  be  spawned  until  a  triangle  with  an  unshared  edge  is  reached  (i.e.  the 
triangles  along  the  noncontoured  boundary) .  The  new  triangles  (which  will  all 
be  outside  the  convex  shell  from  stage  one)  will  then  spawn  other  triangles  to 
encompass  all  of  the  nodes.  Contour  lines  are  drawn  during  stage  two. 
Following  the  spawning  of  all  the  triangles,  edge  contours  are  drawn  within 
triangles  having  an  unshared  edge.  The  unshared  edges  exist  for  all  triangles 
along  the  second  stage  outer  convex  shell  and  the  convex  shell  from  stage  one. 

21.  The  user  must  specify  the  nodes  that  describe  the  noncontoured 
region.  Regardless  of  the  location  of  the  specified  nodes,  a  first  stage 
convex  shell  will  be  formed  which  is  composed  of  only  these  nodes.  Referring 
to  Figure  8,  suppose  the  noncontoured  region  is  specified  by  nodes  5, 6, 2, 9, 
and  7.  (These  nodes  were  identified  during  the  course  of  program  execution  by 
repeatedly  positioning  the  screen  "cursor"  in  the  vicinity  of  each  specified 
node  on  the  noncontoured  boundary.)  Stage  one  mesh  generation  will  form 
triangles  using  only  the  n<">dcs  specified  as  describing  the  noncontoured  region 
plus  any  nodes  lying  on  the  lines  which  connect  the  specified  nodes  (e.g., 
node  8  lies  on  line  7 - 5) . 

22.  The  specification  of  a  noncontoured  region  will  not  yield  the  same 
result  as  simply  not  drawing  (or  erasing)  the  contour  lines  that  would  other¬ 
wise  be  within  the  region.  The  data  points  within  the  noncontoured  region  are 
not  used  for  any  purpose;  therefore,  they  are  not  used  to  form  mesh  triangles 
or  for  contour  interpolation. 

23.  An  example  of  a  contour  plot  with  a  noncontoured  region  is  shown  in 
Figure  9.  Six  points  were  required  to  define  the  region. 


Figure  8.  Describing  a  non¬ 
contoured  region 


Figure  9.  Example  of  a  non¬ 
contoured  region 
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UsinE  Successive  Shells  to  Improve  Efficiency 


24.  McCON  was  originally  written  for  mainframe  processing.  After  the 
program  was  installed  on  a  personal  computer,  it  was  found  that  excessive 
execution  time  occurred  for  problems  with  large  numbers  of  nodes.  For 
example,  approximately  30  min  (1800  sec)  was  required  to  generate  a  mesh  for 
400  nodes.  This  execution  time  was  reduced  to  56  sec  by  using  the  successive 
shell  procedure  to  be  discussed.  The  reason  for  the  excessive  execution  time 
was  the  application  of  the  three  criteria  to  choose  which  node  to  use  when 
spawning  triangles.  As  each  triangle  is  spawned,  all  nodes  must  be  evaluated 
as  candidates  for  each  new  triangle.  This  will  cause  the  execution  time  to 
increase  as  the  square  of  the  number  of  nodes.  In  addition,  following  each 
successful  spawn,  it  is  necessary  to  check  each  and  every  previous  triangle  to 
see  if  the  new  triangle  is  sharing  an  edge  other  than  the  edge  of  the  parent 
triangle.  The  time  for  this  check  also  increases  as  the  square  of  the  number 
of  generated  triangles.  Therefore,  a  procedure  to  reduce  the  number  of  nodes 
to  be  considered  as  candidate  nodes  and  to  reduce  the  number  of  triangles  to 
be  compared  after  a  new  spawn  would  hold  the  promise  of  substantially  improv¬ 
ing  execution  time. 

25.  The  procedure  of  successive  shells  is  basically  an  extension  of  the 
concepts  employed  to  accommodate  noncontoured  regions.  Imagine  some  (any) 
closed  convex  shape  to  be  superimposed  onto  a  set  of  nodes.  If  only  the  nodes 
which  lie  within  or  on  this  convex  figure  are  used  to  form  a  triangle  mesh,  it 
is  assured  that  the  convex  shell  around  the  resulting  mesh  will  be  entirely 
enclosed  by  the  imagined  convex  shape.  It  is  also  assured  that  none  of  the 
remaining  nodes  will  lie  between  these  two  convex  boundaries.  Then  imagine 
another  convex  shape  that  surrounds  the  first  to  be  superimposed  over  the 
nodes.  If  triangles  are  then  spawned  from  the  first  shell  (exactly  as  was 
done  for  expanding  the  mesh  about  a  noncontoured  region)  the  result  will  be  a 
second  convex  shell  that  is  encompassed  by  the  newly  imagined  convex  shape. 
This  procedure  may  be  repeated  with  successively  larger  imagined  convex  shapes 
until  all  of  the  nodes  are  encompassed.  Figure  10  illustrates  the  concept. 

The  triangles  labeled  "a"  were  spawned  within  boundary  1.  The  heavy  line  is 
triangle  shell  #1.  The  triangles  labeled  "b"  were  spawned  from  shell  #1  and 
produced  shell  #2.  It  is  obvious  that  the  procedure  may  be  repeated  to 
include  all  nodes. 
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Figure  10.  Spawning  with  successive  shells 


26.  McCON  is  programmed  to  use  the  procedure  of  successive  shells 
whenever  the  total  number  of  node  points  is  in  excess  of  180.  (A  mesh  formed 
with  180  nodes  requires  approximately  45  sec  for  triangle  generation.) 
Concentric  circles  are  used  as  the  (imagined)  convex  shape  to  build  successive 
shells.  The  centers  of  these  concentric  circles  are  located  at  the  mid-point 
of  the  minimum  and  maximum  coordinates  of  the  nodes.  The  radius  of  each 
successive  shell  is  determined  by  requiring  that  at  least  a  certain  number 
(currently  set  to  60)  of  nodes  are  contained  within  each  shell.  If  a  non- 
contoured  region  is  also  specified,  the  center  of  the  concentric  circles  is 
located  at  the  mid-point  of  the  nodes  specifying  the  noncontoured  region;  the 
distance  from  the  center  to  the  furthest  specified  node  on  the  noncontoured 
boundary  is  computed,  and  the  first  successive  shell  contains  all  of  the  nodes 
lying  within  that  computed  distance  from  the  center. 

27.  The  procedure  of  successive  shells  gave  a  tremendous  reduction  in 
execution  time  for  two  reasons.  First,  the  number  of  nodes  to  be  examined  as 
candidate  nodes  was  reduced  to  approximately  60.  Second,  it  was  no  longer 
necessary  to  examine  every  previously  spawned  triangle  to  check  for  shared 
edges.  Since  each  successive  shell  is  convex,  it  is  assured  that  as  triangles 
are  spawned  within  any  given  shell,  that  edges  can  only  be  shared  with 
triangles  within  or  outside  of  the  previous  shell.  The  successive  shell 
procedure  resulted  in  an  execution  time  that  is  proportional  to  the  number  of 
nodes  rather  than  to  the  square  of  the  number  of  nodes. 

28.  The  successive  shell  procedure  does  have  the  disadvantage  that 
"peculiar"  triangles  may  be  generated  at  the  shell  boundaries.  This  may 
result  in  some  unwarranted  "waviness"  in  contour  line  drawing.  Therefore,  in 
the  event  that  the  contour  plots  are  being  directed  to  a  plotter,  the  number 
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of  nodes  per  spawning  cell  is  set  such  that  it  never  exceeds  one -half  of  the 
number  of  nodes  available  for  triangle  generation.  This  increase  in  the 
number  of  nodes  per  cell  will  lead  to  an  increase  in  execution  time.  However, 
since  most  plotters  draw  so  slowly  (compared  with  drawing  lines  on  the 
screen),  the  execution  time  increase  will  not  be  particularly  noticeable. 
Contour  plots  sent  to  a  plotter  will  (for  large  numbers  of  nodes)  generally  be 
of  a  higher  quality  than  screen  plots. 

Contouring  Nonconvex  Areas  Bv  Composite  Fragments 

29.  The  triangle  mesh  generation  criteria  will  always  result  in  the 
data  points  being  surrounded  by  a  unique  convex  boundary.  This  boundary  is 
generally  adequate  for  topographical  applications.  However,  if  contour  plots 
are  desired  for  applications  such  as  stress  analysis,  the  regions  in  which 
meaningful  contours  may  be  drawn  are  often  limited  by  the  physical  shape  of 
the  object  being  stressed.  These  physical  regions  often  do  not  yield  a  convex 
shape  which  conforms  to  the  physical  boundaries.  Any  nonconvex  shape  may  be 
subdivided  into  two  or  more  shapes  which  are  convex.  McCON  provides  an  option 
to  restrict  contour  drawing  to  an  internal  convex  region.  This  option  is 
similar,  but  opposite,  to  the  option  to  prevent  drawing  in  a  noncontoured 
area . 

30.  Figure  11a  shows  the  contours  of  the  shearing  stresses  within  an 
embankment  and  its  foundation.  (The  small  crosses  or  tics  on  the  drawing 
represent  the  location  of  the  data  points.  The  outline  of  the  embankment  and 
foundation  was  drawn  using  the  "template"  option).  Contour  lines  are  drawn 
outside  the  embankment  boundaries  because  the  generated  triangle  mesh  (shown 
in  Figure  lib) ,  which  must  result  in  a  convex  shape  which  encompasses  all  of 
the  data  points,  contains  triangles  beyond  the  physical  boundaries.  It  is 
obvious  that  the  contour  lines  drawn  through  "space"  are  not  desirable. 

Figure  11c  shows  the  contours  within  the  foundation.  These  contours  were 
created  after  restricting  contour  drawing  to  the  convex  region  formed  only  by 
the  data  points  within  the  foundation.  (The  use  of  the  option  to  restrict 
contouring  to  a  specific  region  required  that  the  user,  by  use  of  the  movable 
screen  "cursor",  identify  the  four  data  points  at  the  corners  of  the  founda¬ 
tion.)  The  triangle  mesh  is  shown  in  Figure  lid.  Figure  lie  shows  the 
contours  which  resulted  from  restricting  contour  drawing  to  be  entirely  within 
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g.  Combined  composites 


h.  Composite  mesh 


Figure  11.  Preparing  a  "composite"  contour  plot 


the  upper  portion.  The  triangle  mesh  for  the  embankment  portion  is  shown  in 
Figure  Ilf.  A  composite  drawing,  as  shown  in  Figure  llg,  may  be  produced  by 
first  contouring  the  lower  portion  and  then,  without  removing  the  paper  from 
the  plotter,  contouring  the  upper  portion. 

Handling  Discontinuity  Or  "Fault"  Boundaries 

31.  The  preceding  discussion  demonstrates  that  the  parent  region,  which 
is  defined  by  all  of  the  data  points,  may  be  broken  into  any  number  of  convex 
regions  simply  by  identifying  which  internal  data  points  are  along  convex 
boundaries  of  the  subregions.  A  composite  contour  drawing  may  then  be 
produced.  McCON  also  provides  a  variation  of  the  above  option  for  specifying 
the  location  of  connected  lines  which  indicate  where  a  contour  discontinuity 
is  expected.  For  example,  if  the  trend  of  a  geologic  "fault"  is  known,  then 
it  would  be  expected  that  the  elevation  of  the  top  of  a  given  formational  unit 
(e.g.  top  of  limestone)  would  be  discontinuous  across  the  "fault"  (i.e., 
contour  elevations  should  not  be  interpolated  between  data  points  that  are  on 
opposite  sides  of  the  "fault").  The  use  of  this  option  requires  that  a  data 
file  be  previously  prepared.  This  data  file  must  contain  x,y  pairs  (one  pair 
per  data  file  line)  which  define  the  "fault"  location.  If  there  are  more  than 
two  x,y  pairs  (more  than  one  line),  then  the  lines  which  connect  the  x,y  pairs 
(in  the  order  given)  must  form  a  convex  shape  always  turning  in  the  right-hand 
direction.  (The  "fault"  line  may  not  have  "switchbacks.") 

32.  Additional  nodes  (which  are  displayed  on  the  screen)  are  created  at 
regular  intervals  along  each  "fault"  line  segment.  The  contour  plot  may  then 
be  created  as  a  "composite"  plot.  First,  invoke  the  option  that  contours  are 
restricted  to  be  within  a  convex  area.  Only  the  nodes  situated  on  the  right- 
hand  side  of  the  x,y  pairs  (from  the  "fault"  data  file)  will  be  used  to  define 
the  contour  plot.  The  elevations  of  the  additionally  created  nodes  along  the 
"fault"  line(s)  are  computed  using  Equation  2  (with  the  Zj  restricted  to 
points  within  the  contoured  area) .  Then,  a  second  contour  plot  is  prepared  by 
invoking  the  "noncontoured  region"  option.  The  non-contoured  area  will 
consist  of  the  nodes  on  the  right-hand  side  of  the  "fault."  The  elevations  of 
the  additional  nodes  are  computed  with  the  Zj  restricted  to  points  outside  the 
non-contoured  area. 
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33.  Figure  12  contains  a  series  of  drawings  which  demonstrate  contour¬ 
ing  in  the  vicinity  of  a  "discontinuity"  boundary.  The  data  represent  water 
table  elevations  on  either  side  of  a  vertical  impermeable  boundary.  Since 
water  cannot  flow  across  an  impermeable  boundary,  interpolation  of  elevations 
across  the  boundary  should  be  avoided.  Figure  12a  shows  "all"  of  the  nodes 
used  to  form  the  triangle  mesh  (original  plus  "created"  nodes).  The  solid, 
bent  line  (which  was  drawn  via  the  "template"  option)  is  located  along  the 
impermeable  boundary.  Figure  12b  shows  the  contour  plot  which  would  result  if 
no  "fault"  boundary  was  prescribed.  The  contour  lines  are  continuous  across 
the  "fault"  boundary.  Figure  12c  shows  the  triangle  mesh  which  resulted  from 
not  including  the  "fault."  The  contour  lines  are  discontinuous  at  the  "fault" 
boundary,  but  are  continuous  beyond  the  ends  of  the  "fault"  line  as  shown  in 
Figure  12d.  The  mesh  for  the  "fault"  boundary  case  is  shown  in  Figure  12e. 

Profiles 


34.  After  a  triangular  mesh  has  been  formed,  it  is  a  relatively  simple 
process  to  create  an  "elevation"  profile  along  any  section  of  the  mesh.  The 
x,y  location  coordinates  and  the  z  value  (elevation)  for  the  vertices  of  each 
triangle  is  then  known.  After  the  user  specifies  the  beginning  and  ending 
locations  of  the  ends  of  the  profile  section,  McCON  computes  the  x,y  location 
of  numerous  (one  location  every  5  screen  units)  profile  points  along  the 
section.  Equation  1  is  applied  to  locate  the  triangle  in  which  the  profile 
point  lies.  The  elevation,  z,  to  assign  to  each  profile  point  is  computed  by 
assuming  that  the  elevation  of  any  point  within  any  given  triangle  is  given 
by: 

z  =  ax  +  by  +  c  (3) 

where  the  constants  a,b,  and  c  are  found  by  substitution  of  the  known  values 
of  x,y,  and  z  at  the  three  vertices  of  the  triangle  in  which  the  profile  point 
lies.  The  profile  drawings  may  either  be  superimposed  directly  on  the  contour 
plots  as  shown  in  Figure  13a,  or  the  (same)  profile  may  be  produced  as  a 
separate  drawing  as  shown  in  Figure  13b.  Several  profiles  may  be  produced  at 
one  time  via  the  "GRID"  option  of  McCON.  Examples  of  multiple  profiles  are 
shown  in  Figure  14a  and  Figure  14b. 


22 


a.  Location  of  "barriers"  and 
elevations  of  nodes 


Figure  12.  Contours  across  a  "fault"  boundary 
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APPENDIX  A  -  USER’S  GUIDE  FOR  MCCON 

A  floppy  diskette  containing  the  executable  (MICR01.EXE  and  MCGRID.EXE)  and 
Fortran  source  language  (MICR01.F0R,  MICR02.F0R,  MICR03.F0R,  MICR04.F0R,  and 
MCGRID. FOR)  programs  is  attached  to  the  inside  back  cover  of  this  report.  A 
batch  file  (MCCON.BAT)  which  calls  the  programs  into  execution  and  a 
demonstration  data  file  (MCDEMO.DAT)  is  also  included.  To  utilize  the  dis¬ 
kette,  follow  these  steps: 

1.  At  the  C:  prompt  create  a  new  directory  (e.g.,  enter  MD\CONTOUR) . 

2.  Get  into  the  new  directory  (i.e.,  enter  CDNCONTOUR) . 

3.  Insert  the  diskette  into  the  floppy  drive  (presumably  the  A  drive) 

4.  Enter  COPY  A:*.* 

5.  Enter  INSTALL  (the  programs  will  be  installed  on  the  hard  drive). 

6.  Remove  the  diskette  and  return  it  to  the  storage  envelope. 

7.  Enter  MCCON  to  execute  the  programs. 

Note:  The  demonstration  data  file,  MCDEMO.DAT,  may  be  executed  by 

simply  entering  a  carriage  return  at  all  program  prompts  (with 
the  exception  of  entering  a  "Y"  when  requested  at  the  outset 
of  the  actual  contouring  activity) . 

The  contouring  programs  are  brought  into  execution  by  entering  "MCCON".  The 
MCCON  command  first  issues  the  command,  MCGRID,  which  requests  the  name  of  the 
(separately  prepared)  data  file  which  contains  the  pertinent  information  and 
performs  certain  enhancements,  if  requested,  on  the  data.  A  file  named 
MAP. DAT  is  always  created  upon  the  exit  from  MCGRID.  The  MCCON  command  then 
issues  the  command  MICR01  which  accesses  the  MAP. DAT  file  and  performs  the  ac¬ 
tual  contouring  functions.  You  may  (at  the  C:  prompt)  enter  MCGRID  (to 
prepare  a  new  MAP. DAT  file)  or  MICR01  (to  operate  on  the  current  MAP. DAT  file) 
directly  from  the  keyboard  if  desired. 

All  of  the  program  queries  listed  below  may  be  answered  by  a  default  (DEF) 
response,  (i.e.,  simply  by  a  carriage  return).  The  default  for  ALL  Y/N 
queries  is  N  (no).  The  following  instructions  give  the  default  response, 
where  applicable,  for  other  than  Y/N  queries. 

Audible  beeps  issued  by  the  program  indicate  that  a  carriage  return  must  be 
entered  to  cause  the  computations  to  proceed. 

The  numbers  m  parentheses  (e.g.  (1.01))  refer  to  the  questions  asked  during 
execution  of  the  contouring  programs. 

o  DATA  INPUT  ACTIVITY. . (1 .01  THRU  1.16).. 

(1.01)  THRU  (1.06)  RELATE  TO  THE  MANNER  IN  WHICH  THE  INPUT 

DATA  FILE  IS  TO  BE  INTERPRETED.  THE  MCGRID  PROGRAM  ISSUES  THESE  FROMPTS . 

(1.01a)  ENTER  (ORIGINAL)  DATA  FILENAME  (DEF=MCDEMO.  DAT) :  Execution  of  the 
programs  requires  that  an  ASCII  data  file  containing  the  x  and  v  coordinates 
and  the  z  (elevation)  values  of  each  of  the  points  be  separately  prepared. 
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This  file  may  be  created  via  an  editor  or  by  another  program.  Each  line  of 
this  file  must  contain  (as  a  minimum)  one  set  of  x,y,  and  z.  The  x  and  y 
coordinates  may  be  given  in  a  random  order.  Multiple  z  values  may  be  included 
on  a  line.  The  order  of  the  x,y,  and  z  values  is  arbitrary;  however  they  must 
be  in  the  same  order  on  each  line  within  a  given  file.  The  default  response 
will  access  an  example  data  file  (MCDEMO.DAT) .  This  file  may  be  treated 
either  as  a  "free-field"  file  or  as  a  "columnar"  file  (see  (1.01b)). 

(1.01b)  IS  FILE  " FREE-FIELD "  OR  "COLUMNAR"  (F(DEF)  OR  C  ?:  The  data  file 
structure  may  either  be  "free-field"  or  "columnar".  If  the  file  is 
"columnar",  the  data  columns  must  be  "right-justif ied" .  The  program  will  ex¬ 
amine  the  input  data  file  and  produce  a  series  of  "X's"  to  denote  the  columns 
as  shown  in  the  following  example. 

(1)  (2)  (3)  (4)  (5) 

XXXXXX  XXXXXX  XXXXX  XXXXXX  XXXXXXX  (computer  made  heading) 

1043.2  653.4  198.4  156.43  234.54  (1st  line  in  data  file) 

986.3  1785.3  67.4  1245.52  (2nd  line  in  data  file) 

etc.  ..  ...  ...  ...  ..  (etc.  ...  ..  ) 

Presume  that  Column  (1)  is  x,  Column  (2)  is  y  and  Columns  (3),  (4),  and  (5)  are 
z  values.  If  a  value  is  missing  (e.g  2nd  line  Column  (4)),  then  any  contour¬ 
ing  operation  (see  (1.07)  thru  (1.13))  which  involves  such  a  missing  value 
will  cause  the  entire  line  to  be  ignored  (i.e.  the  information  in  the  line 
will  be  used  provided  that  the  operations  do  not  involve  a  missing  field). 

"Free-field"  files  are  characterized  by  separating  the  numbers  by  one  or  more 
spaces  or  by  commas  or  by  /' s.  The  sample  data  file  above,  could,  in  free- 
field  format  be  represented  as: 

1043.2,653.2,  198.4,  156.43,234.54 

906.3,1785.3,67.4, ,1245.52 

or  by 

1043.2  653.2  198.4  156.43  234.54 

906.3  1785.3  67.4, , 1245.52 

If  any  data  element  is  missing  in  a  free-field  file  (signified  by  (,,))  then 
the  element  is  interpreted  as  "0.0";  -  not  as  a  "missing"  bit  of  information. 
ALL  data  in  a  "free-field"  data  file  MUST  BE  NUMERIC. 

A  "columnar"  data  file  is  recommended  since  this  form  provides  for  data 
manipulations  that  are  not  possible  with  a  "free-field"  file. 

(1.02)  WANT  INVERSE  POWER  GRIDDING  ?  (Y/N):  This  option  is  provided  so  that 
the  (random)  input  data  points  may  be  manipulated  to  create  a  regular  rectan¬ 
gular  grid  of  data  points.  A  13X13  (15X15  if  the  border  described  under 
(1.05)  is  retained)  point  grid  is  created.  The  z  values  (or  elevations)  of 
the  grid  points  are  computed  by  the  formula: 
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)  J 

Zi  =£  [Zj  /  (di  J )»]  /  £  [l/(du)«]  for  all  dij  ?  0 

i  i 

Zi  =  Zj  for  any  dij  =  0 

where  Zi  is  elevation  at  the  grid  point  i, 

Zj  are  the  elevations  at  the  given  original  points,  j, 
dij  is  the  distance  between  point  i  and  j, 
and  Q  is  an  exponent  chosen  by  you  (usually  2) . 

The  effect  of  imposing  inverse  power  gridding  is  to  smooth  the  data.  This  op¬ 
tion  should  probably  be  used  for  cases  where  the  original  data  points  tend  to 
cluster.  If  good  spatial  coverage  exists,  this  option  is  not  recommended. 

Also  see  (2.32)  below  for  an  alternate  method  for  data  smoothing. 

(1.03)  WANT  ORIGINAL  DATA  POINTS  "AS  IS"  ?  (NO  FILLING 

OUT  TO  EDGE  OF  PLOT)  (Y/N):  The  default  (N)  will  cause  a  rectangular 
border  (which  just  surrounds  the  original  data  points)  of  points  to  be 
created.  The  elevations  (z)  of  the  points  that  are  created  and  lie  on  the 
border  are  computed  as  discussed  in  (1.02)  above.  Generally,  the  addition  of 
these  fictitious  points  will  cause  a  more  pleasant  appearance  to  the  contour 
plot.  A  "Y"  response  will  result  in  leaving  the  data  points  "as  is",  without 
any  border  enhancements.  The  query  for  this  option  is  not  in  effect  when  in¬ 
verse  power  gridding  is  invoked. 

(1.04)  ENTER  VALUE  FOR  Q  (POWER  FOR  INVERSE  FITTING )  DEF=2:  The  Q  value 
defined  under  (1.02)  is  to  be  supplied  to  this  query.  This  query  is  invoked 
only  when  inverse  power  gridding  is  selected.  Higher  Q  values  will  cause  the 
elevation  at  a  grid  point  to  be  more  closely  related  to  its  nearest  (original 
data  point)  neighbor.  The  closer  neighbors  should  indeed  have  more  weight  in 
determining  the  assignment  of  a  grid  point  elevation,  but  within  reasonable 
limits;  therefore  a  good  value  to  select  is  "2"  (the  default  response). 

(1.05)  REMOVE  APRON  FROM  AROUND  DATA  POINTS  ?  (Y/N):  The  default  (N)  will 
cause  a  rectangular  apron  to  be  placed  around  the  original  (or,  when  invoked, 
inverse  power  gridded)  border  (1.03).  No  original  data  points  will  then  lie 
on  the  extreme  edge  of  the  contour  plot.  The  elevations  of  the  fictitious 
border  points  are  computed  as  described  under  (1.02).  A  "Y"  response  will 
cause  the  border  to  NOT  be  generated. 

(1.06)  DO  NOT  HONOR  ORIGINAL  DATA  POINTS  ?  (Y/N):  The  default  (N) ,  will,  when 
inverse  power  gridding  is  selected,  retain  the  original  given  data  points  in 
addition  to  the  gridded  points.  The  default  response  is  usually  recommended 
since  inverse  power  gridding  tends  to  "chop-off"  the  peaks  and  valleys.  A  "Y" 
response  will  not  add  the  original  data  points  to  the  gridded  points  (i.e., 
the  original  data  points  are  deleted) . 
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(1.07)  THRU  (1.14)  RELATE  TO  THE  IDENTIFICATION  OF  THE  DATA  ELEMENTS. 

The  following  discussion  refers  to  the  sample"columnar "  data  file  below. 

1  234  56  7  8  9 

XXXXX  XXXX  XXXX  XXXXXXX  X  xxxxx  xxxxx  xxxxx  xxxxx 

24-A1  6226  3756  5150.23  D  14.54  13.50  23.22 

24-A2  6444  3890  5160.34  D  15.43  23.45  16.45  15.34 

24-A3  6489  3967  5170.33  A  17.34  24.44  13.23 

24- A1  6226  3756  5150.23  D  23.22  12.80 

25- A1  6422  4201  5180.33  A  14.54  12.32  15.44  13.34 

. . .  etc . 

Column  1  may  be  thought  of  as  the  borehole  label,  Column  2  as  the 
y  coordinate.  Column  3  as  the  x  coordinate.  Column  4  as  the  ground 
surface  elevation.  Column  5  as  the  aquifer  identifier,  and 
Columns  6-9  as  depths  to  the  water  surface  on  different  dates. 

(1.07)  WHICH  COL  IS  THE  " IDENTIFIER "  (0  IF  NONE,  DEF-0 ):  The  computer  screen 
will  show  the  first  few  lines  (as  above)  of  the  input  data  file.  You  may 
select  the  column,  if  any,  which  is  to  be  used  as  an  identifier  (i.e.  the 
borehole  label  column  1  above) .  If  there  is  no  column  containing  such  infor¬ 
mation,  enter  0.  If  you  do  select  a  column  as  being  the  identifier  column, 
then  in  the  event  that  the  same  identifier  number  is  repeated  within  the  file, 
only  the  applicable  information  contained  in  the  last  occurrence  is  retained. 
For  example,  borehole  24-A1  is  repeated.  Therefor,  the  information  associated 
with  the  second  entry  will  be  used  for  any  operations  involving  Columns  8  or  9 
and  the  first  entry  will  used  for  any  contouring  operations  restricted  to 
Columns  1  thru  7.  This  option  applies  only  to  "columnar"  data  files.  Al¬ 
phabetic  data  may  be  contained  in  the  "identifier"  column. 

(1.08)  WHICH  COL  IS  THE  ", SELECTOR "  ?  (0  IF  NONE,  DEF=0)i  A  data  file  may  be 
prepared  that  contains  multiple  sets  of  data  to  be  separately  contoured 
provided  that  the  file  contains  a  column  that  identifies  the  different  sets. 

In  the  above  sample  data  set.  Column  5  is  the  "selector"  column.  If  one  were 
contouring  the  ground  surface  elevation,  then  a  choice  of  "0"  to  this  query 
would  be  proper.  If  one  were  wishing  to  contour  the  depth  to  water  within  a 
single  aquifer,  then  the  proper  response  would  be  "5".  This  option  applies 
only  to  "columnar"  data  files. 

(1.09)  ENTER  SELECTOR  " KEEP  VALUE "  ( DEF=0 ):  This  query  appears  only  if  you 
have  responded  with  a  non-zero  answer  to  (1.08)  above.  Enter  the  characters 
which  identify  which  data  set  you  wish  to  contour  (i.e.,  for  the  example,  en¬ 
ter  a  "D"  or  an  "A") . 

NOTE:  The  next  four  queries  (1.10  thru  1.13)  have  default 

responses  that  anticipate  that  the  original  data  file 
has  only  one  x,  y,  and  z  triplet  (in  that  order) 
entered  as  the  first  items  on  each  line. 
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(1.10)  WHICH  COL  IS  "X"  ?  ( DEF=1 ):  Enter  the  number  of  the  column  which  con¬ 
tains  the  x  coordinate  (i.e.  "3"  in  above  example). 

(1.11)  WHICH  COL  IS  "Y”  ?  ( DEF=2 ):  Enter  the  number  of  the  column  which  con¬ 
tains  the  y  coordinate  (i.e.,  "2"  in  above  example). 

(1.12)  WHICH  COL  IS  Z1  ?  (-)  SIGN  IF  DESIRED  (DEF=3):  You  may  contour  the 
value  of  a  single  column  in  the  data  file  or  the  sum  or  difference  of  any 
columns.  A  (-)  sign  before  the  column  number  indicates  a  subtraction  opera¬ 
tion.  Suppose  that  it  is  desired  to  contour  the  elevation  of  the  water  table 
surface  on  the  date  associated  with  Column  7.  Then  a  response  of  ”4"  is 
proper  at  this  time  as  this  will  initially  set  z  to  the  ground  surface  eleva¬ 
tion. 

(1.13)  WHICH  COL  IS  Z2  ?  (FOR  Z1+  OR  -  Z2)  OR  O(DEF)  WHEN  DONE :  Enter  the 
number  of  the  column  (preceded  by  a  sign  if  necessary)  which  is  to  be 
added  or  subtracted  to  Zl.  A  "0"  response  signifies  that  no  further  data 
manipulations  are  wanted.  To  contour  the  water  table  surface  on  the  date  as¬ 
sociated  with  Column  7,  the  proper  response  would  be  "-7"  (i.e.  contour  (4)- 
(7)).  This  query  will  keep  repeating  (up  to  5  times)  until  a  "0"  response  is 
entered. 

(1.14)  ENTER  NEW  CHECKING  DISTANCE  IF  DESIRED  (DEF=7):  The  data  are  checked 
to  see  if  any  points  are  at  or  near  the  same  coordinate  locations.  If  any 
points  are  closer  than  (what  will  ultimately  be)  the  specified  (default=7 
screen  units)  checking  distance,  the  points  are  combined  into  a  single  data 
point  at  the  average  elevation  of  the  affected  points.  Points  which  have 
identical  coordinates  are  combined  regardless  of  any  specified  checking  dis¬ 
tance.  The  default  value  of  7  is  recommended. 

(1.15a)  ENTER  MULTIPLIER  FOR  Z  (RETURN  IF  HAPPY):  Following  all  data  calcula¬ 
tions,  the  minimum  and  maximum  z  values  are  displayed  on  the  screen.  You  may 
multiply  all  of  the  z  values  by  the  constant  you  supply.  The  default  is  1.0. 

(1.15b)  ENTER  MULTIPLIERS  FOR  X  AND  Y  (RETURN  IF  HAPPY):  The  minimum  and  max¬ 
imum  x  and  y  values  (and  their  differences)  are  displayed  on  the  screen.  To 
change  these  values,  enter  a  pair  of  constants  (positive  or  negative, 
separated  by  a  comma  or  space) .  All  x  and  y  coordinates  will  be  multiplied  by 
these  constants.  A  negative  multiplier  will  invert  the  (indicated)  axis.  The 
default  response  (1.0, 1.0)  will  cause  a  file  named  MAP. DAT  (which  contains  the 
x,y,  and  z  values  of  the  points  to  be  contoured)  to  be  created  (or  replaced) 
and  execution  to  be  transferred  to  the  contouring  activity. 

(1.16)  WANT  NORMAL  (N)  OR  FINE(F)  INVERSE  POWER  DISTANCE 

GRIDDING  (DEF=N)  ?  (F  OR  N) :  The  default  (N)  response  will  generate  a 
13X13  point  rectangular  grid  of  points  (see  1.02).  The  elevations  of  those 
points  are  computed  by  inverse  power  distance  averaging.  A  response  of  "F" 
will  cause  a  finer  mesh  to  be  created.  The  square  of  number  of  points  on  an 
edge  of  the  grid  plus  the  number  of  original  points  (if  retained)  will  be  set 
up  to,  but  not  exceeding,  999  points. 


o  CONTOURING  ACTIVITY. . (2.01  THRU  2.44)..  PROGRAM  MICROl  ISSUES  THESE  PROMPTS. 


(2.01)  THRU  (2.20)  REQUEST  INFO  TO  PREPARE  THE  CONTOUR  PLOT 


(2.01)  ENTER  DATA  FILE  NAME. 

OR  ENTER  ' TRIPLE ' 

OR  ENTER  '*D (raw)  ' 
OR  ENTER  '*M(ake)  ' 

OR  ENTER  '*G(rid)  ' 

OR  ENTER  ’*R  (econ)  ' 


FIRST  DEFAULTS  TO  ’MAP. DAT' 

THEN  TO  LAST  FILE  NAMED. 

TO  RE-CONTOUR  A  TRIPLE  TRIANGLE  MESH 
TO  DRAW  A  PREVIOUSLY  SAVED  PROFILE 
TO  CREATE  A  PROFILE  FROM  A  PREVIOUSLY 
GENERATED  AND  SAVED  MESH 
TO  DRAW  GRID  BOX  PROFILES  (A  MESH  MUST 
HAVE  BEEN  PREVIOUSLY  SAVED) 

TO  REDRAW  THE  CONTOURS  FROM  A 
PREVIOUSLY  SAVED  MESH. 


OR  ENTER  'S'  TO  STOP :  Enter  a  file  name  or  one  of  the  following. 

(1)  A  carriage  return  will  default  to  the  file  MAP. DAT  which  is  normally 
created  during  the  previous  phase  of  program  execution.  (2)  If  you  are  re¬ 
contouring  a  "triple  triangle  mesh"  (see  2.32)  enter  "TRIPLE".  (3)  A  response 
of  *D  (see  2.23)  will  enable  you  to  re-draw  a  previously  generated  and  saved 
profile  (see  2.26).  (4)  A  response  of  *M  will  enable  you  to  create  new 

profiles  from  meshes  that  have  been  previously  saved  (see  2.21)  without  having 
to  re-generate  the  contours.  (5)  A  response  of  *G  (see  2.41)  will  provide  a 
means  to  generate  multiple  profiles  (horizontally  and  vertically)  on  the  same 
drawing  provided  you  have  previously  saved  a  mesh.  (6)  A  response  of  *R  (see 
2.24a)  will  redraw  the  contours  provided  you  have  previously  saved  a  mesh  (see 
2.21).  The  contour  intervals  may  be  changed.  (7)  An  "S"  response  exits 
McCON. 


(2.02)  ENTER  COLOR  FOR  LINE  WORK  ( DEF=YELLOW) :  The  default  color  for  line 
drawing  on  the  screen  is  yellow.  Permissible  responses  are  ’RED’, 

’BLUE' , 'WHITE' , 'CYAN' , 'GREEN' , 'YELLOW , ’BLACK’ ,  or  'MAGENTA'. 

(2.02)  ENTER  BACKGROUND  COLOR  (DEF=BLACK) :  Same  as  above  for  background  color. 

(2.03)  ENTER  1  OR  2  WHERE  1  IS  FOR  SMOOTH  CONTOURS  (DEF=1) 

AND  2  IS  FOR  SHARP  CONTOURS:  The  default  is  1.  A 
response  of  2  prevents  smoothing  of  the  contours.  This  response  would  be  ap¬ 
plicable  for  contouring  graded  terraces,  or  embankments  and  the  like. 

(2.04)  ENTER  1  FOR  SCREEN.  2  FOR  PLOTTER  (DEF=1) .:  The  output  may  be  directed 
either  to  the  screen  or  to  a  plotter. 

(2.05)  ENTER  SIZE  REDUCTION  FACTOR  (DEF=100) 

(NUMBER  BETWEEN  1-100):  Drawings  directed  to  the  plotter  may  be 
reduced  in  size  as  a  proportion  to  the  number  selected.  A  response  of  100 
will  give  full  sized  page  drawings. 

(2.06)  PAPER  SIZE  ?  ('A'=SMALL,  'B '-LARGE)  (DEF-A):  This  program  is  designed 
to  drive  small  flatbed  plotters.  An  "A"  response  means  to  plot  on  8.5  inch  by 
11  inch  size  paper.  "B"  indicates  11  inch  by  17  inch  paper.  Depending  on  the 
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configuration  of  your  plotter,  you  may  also  need  to  set  plotter  "switches"  for 
different  paper  sizes. 

(2.07)  WANT  PLOT  IN  SINGLE  COLOR  ONLY  (Y/N) : Multiple-pen  plotters  may  be  set 
to  draw  using  different  pen  colors.  However,  if  you  are  preparing  plots  to  be 
reproduced  in  monochrome,  you  may  wish  to  respond  "Y". 

(2.08)  THE  ABOVE  ARE  USED  TO  CHOOSE  A  SCALE  TO  FILL  THE  SCREEN. 

WANT  TO  CHANGE  THESE  VALUES  ?  (Y/N):  The  program  computes  (and 
displays)  the  minimum  and  maximum  x  and  y  coordinates  of  the  data  points.  A 
scale  is  then  automatically  computed  such  that  the  resulting  contour  plot  will 
"fill"  the  screen.  A  response  of  "Y"  to  this  query  will  enable  you  to  change 
the  scale.  This  option  may  be  important  if  you  are  planning  to  produce 
several  drawings  which  you  wish  to  be  to  the  same  scale.  If  you  respond  "Y", 
you  will  be  queried  to  supply  new  minimum  and  maximum  x  and  y  limits. 

(2. 09)  DO  YOU  WANT  TO  DESCRIBE  AN  INTERNAL  CONVEX  AREA 

TO  CONTAIN  ALL  CONTOURING  ?  (Y/N):  A  "Y"  response  indicates  that  you 
want  to  restrict  the  area  in  which  contours  will  be  drawn.  The  location  of 
all  of  the  original  plus  grid  plus  border  points  will  appear  on  the  screen. 

The  following  question  (2.10)  will  then  be  asked. 

(2.10)  ENTER  t  OF  POINTS  TO  DESCRIBE  AREA  (DEF=0) 

OR  ENTER  NAME  OF  A  "FAULT"  FILE:  Three  or  more  points  are  required  to 
describe  the  area.  The  default  "0"  will  skip  to  query  (2.12).  A  numeric  non¬ 
zero  response  will  cause  query  (2.11)  to  be  asked.  If  you  respond  with  a  file 
name  (of  a  previously  prepared  "fault"  file)  that  file  must  consist  of  a 
series  of  x,y  pairs  (in  the  same  units  as  the  original  data)  which  define  a 
boundary  to  "cut"  through  the  original  (plus  gridded)  points.  The  order  of 
the  x,y  pairs  (one  pair  per  data  file  line)  must  be  such  that  the  original 
(plus  gridded)  points  one  wishes  to  retain  are  on  or  are  to  the  right-hand 
side  of  all  the  lines  connecting  (in  order  given)  the  x,y  pairs.  The  bound¬ 
aries  defined  by  the  x,y  pairs  are  not  required  to  "close".  New  "fictitious" 
points  are  generated  (and  displayed  on  the  screen)  at  regular  spacing  between 
the  x,y  pairs.  The  elevation  (z)  of  these  "fictitious"  points  are  computed 
via  the  inverse  power  equation  (1.02),  but  only  the  "retained  "  points  are 
used  to  make  the  estimate.  The  normal  use  for  specifying  a  "fault"  file  is 
to  break  the  contour  plot  into  two  areas.  The  first  area  is  contoured  by 
specifying  that  there  is  an  "area  to  contain  all  contouring"  (2.09)  and  then 
using  a  "fault"  file  to  define  this  area.  The  contour  plot  will  then  consist 
of  contour  lines  drawn  up  to  the  boundary  defined  by  the  x,y  pairs.  The 
elevations  at  the  "fictitious"  points  along  the  boundary  are  estimated  using 
only  the  elevations  of  the  points  within  the  "area  to  contain  all  contouring". 
The  output  should  be  sent  to  the  plotter.  Then  the  second  area  is  contoured 
by  specifying  a  "non-contoured  area"  as  described  by  (2.12)  below.  Specify¬ 
ing  the  same  "fault"  file  will  result  in  a  another  contour  plot  in  which  the 
contour  lines  are  only  drawn  outside  of  the  "non-contoured  area"  (i.e.,  within 
the  region  excluded  from  the  first  area) .  The  elevations  along  the  (same) 
boundary  points  are  estimated  using  only  the  elevations  of  the  points  outside 
of  the  "non-contoured  area".  Again,  the  output  should  be  sent  to  the  plotter, 
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but  do  not  replace  the  paper.  Generally,  the  composite  contour  plots  will  be 
discontinuous  across  the  "fault"  boundaries. 

(2.11)  DESCRIBE  AREA  BY  POSITIONING  CURSOR  AT  LOWER  LEFT  OF 

POINTS  (IN  A  COUNTERCLOCKWISE  BANNER)  AND  PRESSING  "A". 

(CURSOR  IS  MOVED  WITH  ARROW  KEYS  OR  -,6,  <,  AND  >  KEYS. 

- THE  -,6,<.  AND  >  KEYS  GIVE  FINE  MOVEMENT):  Proceed  around  the 

required  area  in  a  counterclockwise  fashion  by  placing  the  cursor  just  to  the 
lower  left  of  each  of  the  points  which  will  define  the  area.  You  will  see  a 
cross  appear  at  the  selected  points  after  striking  "A"  at  each  point. 

(2.12)  IS  THERE  A  NON-CONTOURED  AREA  (Y/N):  A  "Y"  response  will  let  you 
proceed  to  describe  an  internal  convex  area  in  which  no  contours  will  be 
drawn.  Queries  (2.10)  and  (2.11)  will  follow  a  "Y"  answer. 

(2.13)  ENTER  BEGINNING  CONTOUR,  CONTOUR  INTERVAL,  HEAVY  LINE 

INTERVAL  ft  NO.  CONTOURS.  (0  MEANS  SPAN  RANGE  OF  ELEVATIONS ): 

OR  'N'  FOR  NO  CONTOURING  (DEF=AUTO)  (H  FOR  HELP )  —>:  The  minimum  and 
maximum  elevations  are  shown  on  the  screen.  You  choose  the  beginning  contour 
elevation,  the  contour  interval  (either  as  a  positive  or  negative  value),  the 
"heavy  line"  interval,  and  the  number  of  contour  lines  to  be  drawn.  The 
"heavy  line"  interval  must  be  a  multiple  of  the  contour  interval.  These 
"heavy"  lines  will  plot  in  a  different  color  on  the  video  screen  or  be  drawn 
as  thick  lines  on  the  plotter.  A  zero  entry  for  the  number  of  lines  to  be 
drawn  is  interpreted  to  mean  that  contour  lines  will  continue  to  be  drawn 
starting  at  the  beginning  contour  value,  incrementing  by  the  contour  interval 
until  the  range  of  elevations  is  spanned.  For  example,  presume  that  the 
screen  shows  that  the  minimum  elevation  is  138  ft.  and  the  maximum  elevation 
is  243  ft.  A  response  of  140,10,50,0  would  cause  the  140,  150,  160,  ..230, 
and  240  (ft.)  contours  to  be  drawn.  The  150  and  200  ft.  contours  would  be 
"heavy."  A  response  of  154,2,10,10  would  result  in  the  154,  156,  158,  ..170, 
and  172  (ft.)  contours  to  be  drawn  with  the  160  and  170  ft  contours  to  be 
"heavy."  A  response  of  200,2,2,1  would  result  in  only  the  200  (ft.)  contour 
to  be  drawn.  The  default  response  will  result  in  automatic  scaling.  The  con¬ 
tour  interval  will  be  set  to  be  evenly  divisible  by  2,4,5,  or  10  and  the  num¬ 
ber  of  contours  not  to  exceed  13.  The  default  response  would  cause  the  140, 
150,  ..  240  (ft.)  contours  to  be  drawn. 

(2.14)  LEAVE  BOUNDARY  OF  CONTOURED  REGION  UNDRAWN  ?  (Y/N):  The  default  (N) 
will  cause  the  lines  which  define  the  convex  outer  boundary  of  the  data  points 
to  be  drawn.  All  contour  lines  lie  within  this  boundary.  A  "Y"  response  will 
result  in  no  drawn  outer  boundary  and  the  contour  lines  will  be  left 
"dangling"  at  edges. 

(2.15)  DO  YOU  WISH  TO  DRAW  A  TEMPLATE  ?  (Y/N):  This  option  provides  a  way  to 
enhance  the  resulting  contour  plot  with  some  "art-work".  If  a  "Y"  response  is 
given,  the  following  query  will  appear. 

(2.16)  ENTER  TEMPLATE  DATA  FILE  NAME  (DEF=TMPLT. DA T) :  Enter  the  name  of  the 
file  which  contains  the  information  to  draw  the  "art-work".  This  file  (which 
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must  be  previously  prepared)  consists  of  a  set  of  x  and  y  coordinates  (in  the 
same  units  as  the  original  contour  data)  along  with  an  indicator  to  raise  and 
lower  the  pen.  If  the  file  looked  like: 

100  100  1 
100  150  1 

200  150  -1 

300  200  1 

300  500  -1 

a  line  would  be  drawn  from  (100,100)  to  (100,150).  Another  line  would  be  drawn 
from  (100,150)  to  (200,150).  The  pen  would  then  be  raised  (the  -1  means  go  to 
the  next  coordinate  with  pen  up;  a  1  means  pen  down)  and  positioned  at 
(300,200).  A  line  would  then  be  drawn  from  (300,200)  to  (300,500).  This  op¬ 
tion  provides  a  way  to  add  physical  features  (such  as  building  locations)  to 
the  drawing. 

(2.17a)  WANT  TO  SEE  ELEVATIONS  PLOTTED  ?  (Y/N):  A  "Y"  response  will  cause  the 
z  values  (elevations)  to  be  printed  on  the  drawing.  If  the  data  points  are 
close  together,  the  drawing  may  become  cluttered  when  this  option  is  used. 

(2.17b)  WANT  LABELS  PLOTTED?  (Y/N):  A  "Y"  response  will  cause  the  labels  (or 
identifiers) ,  if  any,  of  the  original  data  points  to  be  printed  on  the  draw¬ 
ing.  Again,  the  drawing  may  become  cluttered  if  the  data  points  are  close 
together. 

(2.17c)  LEAVE  X-Y  AXES  ANNOTATIONS  OFF  ?  (Y/N):  The  X  and  Y  axes  may  be 
labeled  with  tic  marks  and  coordinate  values  if  desired.  The  coordinate 
values  which  are  printed  are  automatically  calculated  as  "nice”  numbers. 

(2.18)  LEAVE  ELEVATION  TIC'S  AT  DATA  PTS  OFF  TOO  ?  (Y/N) 

(ENTER  'A'  TO  INCLUDE  GRID  POINTS  ALSO):  The  default  (N)  causes  small 
crosses  (tic's)  to  be  drawn  at  the  location  of  the  ORIGINAL  data  points.  A 
"Y"  response  will  cause  the  tic's  to  remain  undrawn.  A  response  of  "A"  will 
cause  the  ORIGINAL  PLUS  ADDED  DATA  POINTS  tic's  to  be  drawn. 

(2.19)  ENTER  TITLE  OF  PLOT—  (IF  NONE,  HIT  RETURN):  The  title  of  the  plot  (if 
any)  will  be  printed  on  the  drawing. 

(2.2*1  SKIP  LEGEND  ?  (Y/N):  The  contours  lines  ?“e  labeled  "A",  "B",  etc.. 

The  'egend  (which  will  be  printed  on  the  right  margin  of  the  contour  plot) 
lists  the  numerical  value  of  the  labels. 

(2.21)  THRU  (2.44)  RELATE  TO  OPERATIONS  FOLLOWING  THE 
COMPLETION  OF  THE  TRIANGLE  MESHES  AND  THE  PRODUCTION  OF 
THE  CONTOUR  PLOTS. 

(2.21)  DO  YOU  WANT  TO  SAVE  THIS  MESH  (SO  THAT  YOU  MAY 

LATER  DRAW  OTHER  PROFILES)  ?  (Y/N):  This  query  occurs  after  the  con¬ 
touring  is  completed.  A  response  of  "Y"  will  let  you  save  certain  results  of 
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the  mesh  generation  which  will  enable  you  to  later  rapidly  create  profiles 
along  any  section  of  the  mesh.  If  "Y"  is  the  response  the  following  question 
appears. 

(2.22)  ENTER  A  FILE  NAME  IN  WHICH  TO  STORE  THE  MESH. 

(DEF=M.MSH):  Choose  a  file  name  in  which  the  mesh  generation  data 
will  be  stored.  The  default  file  name  is  M.MSH. 

(2.23  WANT  TO  DRAW  A  PROFILE  ?  (Y/N):  A  "Y"  response  indicates  you  wish  to 
create  a  profile  of  the  elevation  data  from  some  given  point  to  another  point. 
This  query  will  occur  following  the  completion  of  a  contour  plot  or  as  a 
result  of  the  *M(ake)  option  discussed  under  (2.01).  Further  questions  about 
the  nature  of  the  profile  will  ensue. 

(2.24a)  WANT  TO  REDRAW  THE  CONTOURS  ?  (Y/N):  A  "Y"  response  will  permit  the 
contours  to  be  redrawn  (provided  a  "mesh"  was  previously  saved) .  The  contour 
interval  may  be  changed  if  desired. 

(2.24b)  WANT  THE  VOLUME  UNDER  THE  MESH  ?  (Y/N):  The  volume  of  the  space  con¬ 
tained  between  some  reference  elevation  (next  query)  and  the  contoured  eleva¬ 
tion  for  the  entire  mesh  area  may  be  computed.  The  units  of  the  volume  are  in 
terms  of  the  units  implied  in  the  original  data  file. 

(2.24c)  ENTER  REFERENCE  ' BASE '  ELEVATION  (DEF=0.0):  The  minimum  elevation  of 
all  data  points  is  shown  on  the  screen.  Enter  the  elevation  you  wish  to  be 
used  as  the  reference  for  volume  calculations.  The  default  response  will 
cause  an  elevation  of  0.0  to  be  used  as  the  reference.  After  the  profile  is 
drawn,  the  volume  will  be  printed  at  the  bottom  of  the  drawing  as  will  the 
area  of  the  mesh  and  the  average  elevation  (the  average  elevation  is  the 
elevation  which,  if  used  as  the  reference  elevation,  would  give  a  volume  of 
zero) . 

(2.25a)  ARE  PROFILE  LINE  ENDPOINTS  TO  BE  LOCATED  WITH  THE  CURSOR  (C) 

OR  FROM  A  NAMED  FILE  (F)  ?  ( DEF=C ):  The  endpoints  of  the  profile  line 
may  be  located  either  by  positioning  the  crosshair  cursor  or  by  supplying  the 
(x,y)  coordinates  of  the  two  ends.  If  the  response  is  "C"  (for  cursor), 
queries  (2.25b),  (2.25c)  and  (2.25d)  will  be  issued.  If  the  response  is  "F" 
(for  file) ,  you  will  be  asked  the  name  of  a  one  line  file  containing  the  x-y 
coordinates  of  the  endpoints  (i.e.,  Xi ,yi  -  X2,yz)  or  you  may  respond  that  the 
file  name  is  CON  (for  console)  and  directly  type  in  the  endpoint  coordinates. 

(2.25b)  WANT  CURSOR  TO  INDICATE  NEAREST  (TO  UPPER  RIGHT) 

DATA  NODE  POINT  ?  (Y/N):  The  data  node  points  will  be  displayed  on 
the  screen.  A  "Y"  response  will  cause  the  location  of  profile  endpoints  to  be 
positioned  exactly  on  the  node  point  located  closest  (to  the  upper  right)  of 
the  crosshair  cursor.  A  "N"  response  will  position  the  profile  endpoints  at 
the  position  of  the  crosshair  cursor. 
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(2.25c)  POSITION  CURSOR  AT  1ST  END  OF  PROFILE  LINE  AND  BIT  RETURN 

..USE  <,>,  ‘ ,  AND  6  KEYS  TO  'FINE'  POSITION  CURSOR :  Position  the  cursor 
at  one  end  of  the  line  defining  the  section  along  which  you  want  a  profile. 
Move  the  cursor  with  the  "arrow"  keys  or,  for  fine  movement,  with  the  <,>,', 
or  6  keys.  Strike  the  "return"  key  when  cursor  is  properly  located. 

(2. 25d)  POSITION  CURSOR  AT  2ND  END  OF  PROF.  LINE  AND  BIT  RETURN :  Same  as 
above  except  to  locate  other  end  of  profile  line. 

(2.26)  Want  profile  (D)rann  on  this  plot,  sent  to 

(F)ile  or  (B)oth  ?  (Enter  D,F  or  B)  DEF=D:  When  this  question  is 
asked,  a  plot  of  the  data  points  will  be  visible  on  the  screen.  A  response  of 
"D"  will  cause  the  profile  to  be  drawn  (after  a  question  regarding  the  scale 
of  the  drawing)  directly  on  the  current  plot  skewed  along  the  positions  indi¬ 
cated  via  query  (2.25a).  A  response  of  "F"  will  send  the  profile  data  to  a 
file  (which  is  named  by  (2.27).  Profiles  saved  via  the  "F"  response  may  be 
drawn  later  by  responding  *D(raw)  to  query  (2.01).  A  response  of  "B"  will 
cause  both  actions  ("D"  and  "F")  to  occur. 

(2.27)  ENTER  A  FILE  NAME  (DEF-D. PRO) :  Enter  a  file  name  in  which  to  store  the 
profile  data.  The  file  name  entered  here  will  be  required  by  the  *D(raw) 
response  to  query  (2.01)  . 

(2.28)  ENTER  A  ONE  LTNE  DESCRIPTION:  The  description  you  enter  will  be 
retained  in  the  file  named  above  (2.27)  and  will,  when  the  profile  is  redrawn 
via  the  *D(raw)  response  to  (2.01)  be  printed  at  the  top  of  the  drawing.  If 
no  description  is  desired  press  the  "return"  key. 

(2.29)  ENTER  SECTION  LABEL  (LIKE  A  FOR  A-A):  The  section  label  you  enter  will 
be  retained  in  the  file  named  by  (2.27).  This  label  will  also  be  placed  at 
the  ends  of  the  profile  section  on  the  drawing  showing  the  data  point  loca¬ 
tions  . 

(2.30)  ENTER  Y-AXIS  LABEL  ( DEF=ELEVA TION) :  When  the  profile  is  redrawn  via 
the  *D(raw)  response  to  (2.01)  this  description  will  be  printed  as  the  label 
of  the  "elevation"  axis. 

(2.31)  ENTER  YOUR  MIN  AND  MAX  PROFILE  ELEVATIONS 

(DEF=MIN  AND  MAX  ELEVATION  OF  ALL  POINTS  ON  PLOT)  :  The  minimum  and 
maximum  elevations  for  the  all  data  points  within  the  entire  mesh  and  along 
the  profile  line  are  shown  on  the  screen.  The  values  you  chose  will  influence 
the  "scale1  of  the  resulting  profile. 

(2.32)  WANT  TO  GENERA TE  A  TRIPLE  TRIANGLE  MESB  ?  (Y/N):  This  option  provides 
yet  another  method  to  "massage"  the  original  input  data  (i.e.  an  alternate  to 
the  inverse  power  scheme  discussed  under  (1.02)).  All  contouring  is  ac¬ 
complished  by  creating  a  mesh  of  triangles  which  connect  every  data  point 
(original  plus  grid  points  if  any).  The  effect  of  a  positive  response  to  this 
query  ' s  to  add  a  (fictitious)  data  point  at  the  centroid  of  each  of  the  "old" 
triangles.  A  file  named  "TRIPLE"  which  consists  of  the  old  plus  the 


All 


"fictitious"  points  is  created.  When  the  program  cycles  back  to  (2.01)  you 
must  respond  with  "TRIPLE"  to  create  a  new  contour  plot  that  considers  the  new 
points.  This  option  may  also  be  used  to  refine  "inverse  power"  grids.  If  the 
response  is  "Y"  the  following  question  is  asked. 

(2.33)  LINEAR  OR  INVERSE  SQUARE  RULE  ?  (L  OR  I)  DEF=L:  The  elevation  of  the 
centroid  of  the  "triple  triangles"  may  be  computed  in  two  ways.  A  "L" 
response  will  compute  the  elevation,  Zc ,  of  the  centroid  via  the  equation, 

Zc  =  a  +  bx  +  cy,  where  the  constants  a,  b,  and  c  are  determined  by  substitu¬ 
tion  of  the  x,  y,  and  z  values  at  the  nodes  of  the  triangle  containing  the 
centroid.  A  response  of  "I"  will  cause  the  centroid  elevation  to  be  computed 
via  the  "inverse  power"  equation  given  under  (1.02)  using  only  the  three 
points  which  define  the  triangle.  The  default  is  "L". 

(2.34)  DO  YOU  WANT  TO  DRAW  TEE  MESH  OF  TRIANGLES  ?  (Y/N):  A  response  of  "Y" 
will  cause  the  triangle  mesh  to  be  drawn. 

(2.35a)  WANT  TO  DRAW  A  PREVIOUSLY  SAVED  PROFILE  ?  (Y/N):  This  is  the  query 
you  will  get  by  using  the  *D(raw)  response  to  (2.01).  Answer  "Y"  to  create 
page  size  profiles  of  "saved"  profiles. 

(2.35b)  1ST  PROFILE  OF  A  SERIES  ON  SAME  PLOT  ?  (Y/N):  More  than  one  profile 
may  be  drawn  on  the  same  plot.  This  option  will  provide  a  means  to  retain  the 
same  scaling  parameters  for  multiple  profiles.  Answer  "Y"  for  the  first 
profile  to  be  drawn;  "N"  responses  from  that  time  on  will  cause  subsequent 
profiles  to  be  drawn  on  the  same  screen  (or  page). 

(2.36)  ENTER  SAVED  DATA  FILE'S  NAME  (OR  "Q"  TO  QUIT)  DEF=D.PRO:  Enter  the 
name  of  a  previously  saved  "profile"  data  file  or  enter  "Q"  to  return  to  query 
(2.01)  . 

(2.37)  WANT  THE  PROFILE  STRETCHED  TO  FILL  THE  PAGE  ?  (Y/N):  The  default 
response  (N)  will  let  you  draw  this  and  subsequent  profiles  to  the  same 
horizontal  scale,  i.e.,  if  one  profile  is  longer  than  another,  the  drawings 
will  also  be  of  different  widths.  A  response  of  "Y"  will  cause  each  profile 
to  fill  the  page  width. 

(2.38)  ENTER  YOUR  ZMIN  AND  ZMAX  (DEF-AUTO) :  The  minimum  and  maximum  eleva¬ 
tions  encountered  along  the  profile  section  are  shown  on  the  screen.  Enter 
the  values  you  wish  to  see  as  the  minimum  and  maximum  values  for  plotting  pur¬ 
poses.  The  default  response  will  cause  automatic  scaling. 

(2.39)  HOW  MANY  Z  UNITS  PER  PLOT  INCH  (DEF=AUTO) :  Enter  the  number  of  verti¬ 
cal  (elevation)  units  per  (approximate)  plot  inch.  The  resulting  plot  may  not 
be  more  than  5  inches  tall.  An  entry  which  would  result  in  a  plot  in  excess 
of  5  inches  will  be  rejected  and  this  query  will  again  appear.  The  default 
response  will  cause  automatic  scaling. 

(2.40)  ENTER  NAME  OF  "MESH”  FILE  ( DEF-M.  NSH  OR  "Q"  TO  QUIT):  This  query  will 
result  when  responding  either  *M(ake)  or  *G(rid)  to  (2.01).  Enter  the  name  of 
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the  file  in  which  mesh  data  was  previously  stored  (2.22).  The  default  will 
assign  "M.MSH"  as  the  file  name. 

(2.41)  WANT  TO  DRAW  GRID  BOX  PROFILES  ?  (Y/N)‘.  This  option  permits  the  simul¬ 
taneous  creation  and  drawing  of  multiple  profiles  located  within  a  rectangular 
grid  box  of  your  choosing.  A  response  of  "Y"  will  cause  a  request  for  the 
"MESH"  (see  2.40)  file.  You  will  also  be  asked  (see  2.41a)  to  describe  the 
box  dimensions  in  which  the  profiles  are  to  reside. 

(2.41a)  ENTER  PROFILE  GRID  BOX  ORIGIN  OF  COORDINATES,  ANGLE  OF  BOX  TO  BORZ, 
LENGTH  OF  BOX ,  i  HEIGHT.  (X,  Y,ANG,L,H) 

OR  ENTER  'S'  FOR  SAME  AS  LAST  TIME  (DEF=AUTO):  Enter  the  coordinates  of  the 
lower  left  corner  of  the  box  in  which  the  grid  profiles  are  to  be  drawn  (in 
terms  of  original  data  coordinates)  ,  the  angle  the  box  makes  with  the 
horizontal  (usually  0) ,  the  length  of  the  box  and  the  height  of  the  box.  En¬ 
ter  'S'  to  use  same  grid  box  as  previous.  The  default  response  will  generate 
a  box  which  fits  inside  the  span  of  the  data. 

(2.42)  WANT  PARALLEL  (TO  BASE)  OR  NORMAL  GRIDS  ? 

INPUT  P  OR  N  (DEF=P):  Respond  appropriately  for  the  required  set  of 
profiles.  A  response  of  P  will  cause  the  profiles  to  be  drawn  parallel  to  the 
base  of  the  grid  box.  A  response  of  N  causes  the  profiles  to  be  drawn  normal 
to  the  base  of  the  grid  box. 

(2.43)  HOW  MANY  PROFILES  ON  THE  PAGE  ?  (DEF=11):  Enter  the  number,  n,  of 
profiles  to  be  created.  The  profiles  are  evenly  spaced  with  the  1st  and  last 
profiles  being  drawn  just  slightly  inside  the  appropriate  edge  boundaries  of 
the  mesh.  The  default  response  of  "11"  will  result  in  11  profiles  drawn 
(horizontally  or  vertically)  centered  on  the  middle  of  the  region.  A  series 
of  'n'  files  will  automatically  be  created  (for  further  processing  under  the 
*D  main  menu  (see  2.01)  option).  These  files  are  named  with  the  same  root  as 
the  ??.MSH  file  (see  2.41)  with  the  extension  .OP,  .IP,  etc.,  (for  profiles 
drawn  parallel  to  the  base  of  the  grid  and  .ON,  .IN,  etc.  for  profiles  drawn 
normal  to  the  base  of  the  grid) .  That  is,  if  the  mesh  file  is  named  M.MSH  and 
5  profiles  parallel  to  the  base  of  the  grid  box  are  to  be  created,  the  profile 
files  will  be  named  M.OP,  M.1P,  M.2P,  M.3P,  and  M.4P.  These  files  may  then  be 
selected  (see  2.36)  to  prepare  individual  profiles  under  the  *D  option  (see 
2.01)  . 

(2.44)  ENTER  SCALE  EXAGERRATTON  FACTOR  (DEF=1.0):  The  default  response  (1.0) 
will  result  in  the  'height'  of  the  profiles  to  be  about  1/2  inch  (maximum). 

The  value  of  your  response  is  used  to  cause  an  exaggeration  to  the  height  of 
the  profiles. 
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APPENDIX  B.  FORTRAN  LISTING  FOR  McCON 


. A  CONTOURING  PROGRAM . 

WRITTEN  BY 
JOHN  B.  PALMERTON 
GEOTECHNICAL  LABORATORY 
UATERWAYS  EXPERIMENT  STATION 
U.S.  ARMY  CORPS  OF  ENGINEERS 
VICKSBURG,  MISSISSIPPI 
AUG  1901 


C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 

C. .THIS  PROGRAM  REQUIRES  A  SET  OF  (N)  X-Y-ELEV  DATA  POINTS 
C..A  GRID  IS  CREATED  WHICH  CONNECT  THE  DATA  POINTS  INTO  TRIANGLES 
C.  THESE  TRIANGLES  WILL  COMPLETELY  OCCUPY  THE  CONVEX  REGION  DEFIN 
C.  BY  THE  EXTERIOR  POINTS.  THE  GRID  IS  CREATED  BY: 

C.  1.  CREATING  A  SEED  TRIANGLE  WITH  EDGES  1-2-3  LABELED  COUNTERCLKWSE 
C.  2.  IN  ORDER  OF  EDGES,  FORMING  NEW  TRIANGLES  BY  SELECTING  THE  POINT 
WHICH  LIES  OUTSM  OF  THE  EDGE  OF  THE  GENERATING  TRIANGLE 
AND  FORMS  THE  LARGEST  ANGLE  WITH  THAT  EDGE. 

REPEATING  THIS  PROCESS  FOR  ALL  EDGES  OF  ALL  TRIANGLES  UNTIL 
EDGES  ARE  SHARED  BY  TWO  TRIANGLES  OR  A  CONVEX  BOUNDARY  IS 
FORMED. 

C.. CONTOURING  IS  DONE  BY  ASSUMING  A  LINEAR  INTERPOLATION  OF  ELEV 
C.  BETWEEN  TRIANGLE  VERTICES.  AS  EACH  NEW  TRIANGLE  IS  FORMED  THE 
C.  CONTOUR  LINES  WHICH  CROSS  THE  SHARED  EDGE  ARE  ARE  DRAWN  AS  EITHER 
STRAIGHT  LINES  OR  CIRCULAR  SEGMENTS.  (IF  THE  USER  WISHES  TO 
INSURE  THAT  ALL  CONTOURS  ARE  STRAIGHT  LINES  BETWEEN  INTERPOLATION 

^Bi?to^Ti^LmygE[^DFiHRSi.is  useful  for  *epresenting 


INED 


C. 

C. 

C.  3. 
C. 

C 


C. 

C. 

1: 


C. 

1: 

c. 

c. 

c. 


.DEFINITION  OF  SOME  PARAMETERS 


Mfl’f11 

VERTEXtM, I) 


POINTS 

ERED 


NODE(I)  . LINKED  LIST  OF  THE  NUMBERS  OF  THE  ENDPOINTS 

WHICH  HAVE  THE  OTHER  END  AT  POINT  I. 

^NC^CLtSOOMANNS.FACT.PX.PY.PZT.CLS 

IPOW’  0  ,  , 

/NcON , LCON ( 999 ) , MED , HCON 

/D^V.PCX, PCY.PSFT. COLOR, CLRBCK,PSIZE,LPEN(8) 

vv.„,v.i/,nvi.v.>/COLl.COL21C0L3,C0U 
COHMON/WINDO^XWMlflj  XWmAx.  YWHIN,  YWMAX 
. . . .  . DNAXxTyDMIN.YUDMAX.FSOFT 


COMMON/DAREA/XyDNIN, XUDMAX. YUDnlN. 
DIMENSION  X(99?],Yj$?9),Z($99] 
CHARACTER  IANSH.OEVH.CLSH,  JNC‘1 
CHARACTER  FNME  li,BNME*12,FMflE  14 
CHARACTER  C0L0R*4.CLRBCK*J,P$IZEJ1 

mm  wio 

CHARACTER  C0L1*4.C0L2*4,CQL3’4 
?  ^BLAC’/RED  ’, 

DATA  IFIRST/O/ 


5933 


5938 


88 


VTA  IF 
JEL=25, 

FNME=,  ’ 

BNME=  MAP.OAT 

feku«$ 

IFflANS.EQ.  YM 
IF  1  IANS . EO .  y  j 
GO  TO  5933 

mm 

DO  88  1=1.72 
DO  88  J=1.54 
IBNll, J’=0 

m\to. 

YWMIN=0, 

YUflAX =780 . 


A . C0L4  *4 

’ CYAN’ . ’ GREE ’ , ’ WHI T ’ . ’ MAGE’ , '  BLUE ’  / 


ENTER  A  Y  )  ’ 

GO  TO  5938 
PRINT* /SET  CAPSLOCK 


B1 


T.Nfi- ' 

hc'on=: 


PR 

PR 


NT'. . PLEASE  NOTE  ! . 

NT'., FOR  QUESTIONS  THAT  ARE  TO  8E  ANSUERED  Y  OR  N 
NT’., ANT. RESPONSE, | INCLUDING  A  CARRIAGE  RETURN)  OTHER  THAN  . 
NT’..  Y  MEANS  NO  .A  BEEP  MEANS  HIT  RETURN  TO  CONTINUE. 
NT’,'  2.01 r 


PRINT’.  ENTER  DATA  FILE  NAME.  FIRST  DEFAULTS  TO  ’MAP. DAT’ 
.  ,,  THEN  TO  LAST  FILE  NAMED.  , 

RIPLE  ,  TO  USE  TRIPLE  TRIANGLE  MESH 
Diray)  TO  DRAW  A  PREVOIUSLY  SAVED  '. 


INT’ 

INT’.  OR  ENTER 
INT’. OR  ENTER 
+  PROFILE 
PRINT*.  OR  ENTER 
+’  PREVIOUSLY 
PRINT1,’ 

PRINT’, 'OR  ENTER 
♦  MUST 
PRjNT* 


*H(ake)”  TO  CREATE  A  PROFILE  FROM  A’. 


*G ( rid ) ’ 


PRIN 
PRIN 
PRIN.,, 
READ 
IFfFNME.EQ 
REWIND  0^ 


’.OR  ENTER  ”*R(econ)’ 


OR  ENTER 
HE 

THEN 


GENERATED  AND  SAVED  MESH*  , 
v  TO  DRAW  GRID  BOX  PROFILES  (A  MESH’, 

HAVE  BEEN  PREVIOUSLY  SAVED)’  . 

'  TO  REDRAW  THE  CONTOURS, FROM  A' 

,  PREVJOUSLjf  SfVED  MESH.’ 


W 


1881 


PEN 
CLOS 

PN 

1882  OPEN 
CLOS 

1883  OPEN 

1884  O^N 

1885 


09 


ATUS=’OLD’ .ERR-1881) 
/ii  ilk'  iuu  »vimyy=  OLD  jERR-1882) 


>7ATUS= ’ OLD’ ,ERR=1885) 


U-  ^T\VH  J 

w/.viATDS  =  VDE(.uii.  ,  ,  , 

9. PlLE=SKUNKv , STAtUS= ’ OLD ’ , ERR=1884 
09 ,  STATyS=rDELI^Tcr^ 1 
9.PlLE=vSAyC0N 
09.STATUS-’DELEfE 

.losei^Mtus^deleiE 

1886  mm 

1887  CONTINUE 
STOP 
ENDIF 

IF(NCOLER.EO.O  THEN 
NCOLER=} 

PRINT*, [COLOR  CHOICES  ARE: 

^(■ilTtR  COLOR  FOR  LINE  WORK  (DEF=YELL0W)  > 

-to 


reoa  ... 

BLACK, G: 


MAGENTf.CYAN’ 
AND  WHITE 


8904 


col6r.eo.,red  COLOR^’RED  , 

COLOR. EQ.  GREN, )  COlOR=  GREE’, 

COLOR. NE. 'RED  ’ .AND. COLOR. NE. ’BLUE’ .AND. COLOR. NE.  WHIT’ .AND. 
AA_OR.NE. IMAGE; . AND. COLOR. NE., YELL;  .AND.  COLOR.  NE,  CYAN’ .AND. 
+C0LOR.NE. ’GfEf  .AND. COLOR. NE.’BLAC’]  C0L0R=YELL’  ,  , 

PRINT,, (A\JY  12.02]  ENTER  BACKGROUND  COLOR  (OEf =BLACK )  >  ’ 
READ  ’(A4r.tL.RkK 
IF  CLRBCK.Et. ,REO’ )  CLRBCK=  RED  ’, 

IF  CLRBCK.EO .GREN  )  CLRBCK=GREE’ , 

IF  CLRBCK.NE.  RED  , AND. CLRBCK.NE.  BLUE  , AND. CLRBCK.NE.  WHIT  , 
♦AND. CLRBCK.NE. [MAGE  .AND. CLRBCK.NE.; YELL  , AND. CLRBCK.NE,  CYAN'. 
♦AND.CLRciCK.NE. ’GREE’. AND. CLRBCK.NE. ’8LAC’  CLRBCK=BLAC’ 

DO, 8904  1=1.5  ,  ,, 

IF  COLOR. EO. COL  I 

IFICLRBCK.EQ.COL(I) 
p8N7^o8EI=JJ, JJ+8 

‘Mu 


JJ=I 
KK= I 


^^l^)L=l-8 
IFTcCLI. NE. COLOR. 
CONTINUE 

DO  7702  I=KK , KK+8 


.COLOR. AND. C0L1.NE.CLRBCK)  GO  TO  7701 
KK+8 
)  L=L-8 


^{l.GT., 

,  m  6l.2ENE ! COLOR. AND. C0L2 . NE . CLRBCK . AND. C0L2.NE. C0L1)  GO  TO  7703 

7702  CONTINUE 

7703  00  7704  1=1.8 
COL3=COLI I ) 

^ IF' COp. NE jCg^OR.  AND. ML^ . NE^C! '  RBCK.  AND. COL3.NE.COL1. 

77O4*C0NTINUE°L3'N  'C°L2  G  °  77° 

7705  COL 4 =COL OR 


B2 


If  .  . 

,  :NnE.FQ. 
FNHE.tG.  *0, 
FNHE.EQ.  *H 
FNHC  =n  ,c' 


FNME=BNf1E 

GO  TO  7777 
GO  TO  7777 


9. FI 


pi  n  m 

LE=FNhE,STATU$=,OLO’ , ERR=567 ) 


567  INO-O 


568 


189 


CONT 
READ 
READ 
READ 

Ife, 

NRilnsI^ 

CLOSE  1 89 

saw? 

RIB: 

PRINT* 

PRINT* 

PRINT* 

PRINT* 

HHP:: 

mu  , 

k  iNQ.&o.g)  then 

PR  NT*!’  —DATA  FILE  ’  ,FNHE,  ’  NOT  FOUND- 
TQ  5 

^;^:81tahFile  is  \fnhe,’  .. 


HIS  PlIoGRAH  EXPECTS  UP  TO  999  NODES  IN  FILE  HAP. DAT 
HAP. DAT  IS  A  FREE  FIELD  FILE  OF  X.Y.Z  VALUES- -1  SET 
PER  LINE.  THE  POINTS  MILL  BE  SCALED  TO  FIT  SCREEN.’  , 
TEMPLATES  HAY  BE  DRAUN  ON  THE  CONTOUR  HAP.  THE  TEMPLATE 
QATA  HUST  BE  IN  A  SEPARATE  FILE;--  A, FREE  FIELD  FILE9, 

F  -X.Y.PEN-  VALUES.  A  NEGATIVE  PEN  VALUE  INDICATES 
HE  END  OF  A  TEHPLAtE  (I.E.  PEN  IS  UP  TO  NEXT  POINT).’ 
-  COORDINATE  SCALE  AS  SCALE  FOR  HAP. DAT)' 


TH _ 

(USE  SAHE 


NC-1 

NC=1 

NC=2 


_  vFNHE.EO. ’TRIPLE’ )  IP0W=1 

876  CO)lT{NyE 
JNC=  1 

NC=0  ,  , 

IF}  JNC.EO. 

IF  JNC.EO.  1 
IF  JNC.EO. ’2’  NC=2  , 
IF(nC.EQ.1.0R.NC.EQ.2)  GO  TO  877 

877  CONTINUE 
REMIND  09 

Sir-1  °del=i- 

7777  CONTINUE 

PRINT*,, FPR, THE  OUTPUT . 

fRINT  'A 

.  , 

PRINT*,  . WAIT 

HH=1 


TER  1  FOR  SCREEN,  2  FOR  PLOTTER  (DEF=1 )  >’ 
,  JNC 


,  JC.EO.  1 
JNC.EO.  2 
JNC.EO.  3 
HH.ES. i 
HH.EO  • 

Remind  o§ 

PC Y=  1 . 0 
PCX-PCY  , 

5008  PRINT’ 

PR  NT 
READ 

lik'd; 


HH=1 
HH=1 
HM=2 
HH=3  , 

v=, ibhh! 

.V=  HP4 
EV=  EPS  ’ 


ENTER  SIZE  REDUCTION  FACTOR  (DEFflOO)’ 
[NUH.  BETMN  1-100)  FOR  PLOTTER  >’ 

PC Y= 100 . 

m-i 

PERt 


B3 


rcEAjliJ.  . FRR-500S !  PCY 
CL0E,E(u3. ST ATUS=' DELETE'  i 


ifipcy.eq.o. 

pcy=pcy/i6o. 

PCX--PCY, 

ODTklT  1  1  \  ! 


PC  Y= 1 00 . 


mkk 


|pER  SIZE  ?  ( ”  A”=SHALL.  ” B”  =LARGE )  ( DEF =A )  )  ’ 

E!,i;f  PSIZE=’A’ 

0,  9  PCX=. 825  PCY 

J\,j2.07)  “ANT  PLOT  IN  SINGLE  COLOR  ONLY  ?  (Y/N)  >  ’ 


IF|PSI|E.EQ,  p  1  PCX=. 825’PCY 

REAd[’  222)  iLi2' 07  WANT  PL°T  IN  $IN6LE  C0L0R  °NLY  ?  lY/Nl  }  ’ 

IFIUNOE  ’f]  THEN 
00  8004  1=1.8 
8004  LPEN[I)=I 

PRINT*  PEN1  $  USED  JO  DRAW  CONTOURS  AND  LEGEND.’ 

PgtNT*, ,PEN2  IS  .  .  .  DATA  TIC  POINTS  S  ELEVATIONS.’ 

PRINT',  PEN3  5  .  TEMPLATES. 

PRINT',  PEN4  S  .  .  .  PROFILES.’ 

PRINT*., PEN4  IS  .  .  .  HORIZONTAL  (PARALLEL)  GRID  PROFILES’ 

PRINT;,, PE^2  IS  VERTICAL  (NORMAL)  GRID  PROFILES* 

Print*’ *so,make  sure  you  know  what  is  in  the  pen  holder!’ 

CALL  UfeELL 
CALL  UPAUSE 

IF(IANS.EQ.'Y’)  THEN 
DO  8002  1=1,8 
8002  LPENlIM 

Bi 

XUDMIN=1Q. 'PCX 
XUDflAX=75.*PCX 
YUD(1IN=10 .  PCY 


YyDMAX=75.  *PCY 

IHFNHE.EO.  ty.OR.FNNE.EQ.  *H’ 


+  IF(FNHE^Q()’*H’]  CALL  INPR0?1ha1, ZHIN,HN, XHINX, XHAXX, YNINY. 

+IF?MAXY  SM'O  oi '  CALL  FENCE(ZnAX-ZnlN’n,1’X,1INX-!<,1AXX'Y,1INY- 

+IF(FNNE!eo:^R  )  CALL  FENCE ( ZMAX , ZMIN, MM, XMINX, XMAXX, YMINY. 
^  YHAXY, 5N, 0. 1 J 
CALL  PEftD 
GO  TO  5 


GO  TO  5 
f PC^LL=60 

222  i^hnteo|  ipce 

CALL  ORAWCONf IPC 

$litaPicLJ)  PRIN 
IF (MM .£0 . 3 J  STOP 
CLOSE  09 
GO  TO  5 


mb* 


’ENTER  SPRINT  TO  SEND  TO  PLOTTER’ 


NE  ,CONTER(ZINT) 


mm 


COHHON/CELLS/NUNTIoOO 

!M  Imm 

COMNON/NODEL/NODEI170 


(9?9),Y(999i.Z(999),cX(3).CY(3),CZ(3), 
.?RNriOOjl ,  IRD(IOO)  ,KRO,AVX,AVY 


ZUINDO/JfWMIN.  xwmAx,  ywAin.  ywmax 


CHARACTER  COLri.COLC’A.rn:  '  4 .  ( 01  4 ’ A 
CHARACTER  0EVVC0L0R,4.uR6CI'‘4.PRIZE'! 
INTEGER  SIDE. VERTEX. FI. F2.P3.F 
0PEN151.FILE=  SAVCONS  I 
“  "  '  '  2000 


=0 

=2,KRD 


isr1"- 


AVV!”2 


RD{ JJ-1 ) . ANO. IO.LE. IRD(JJ) )  THEN 


..  , _ 2  CALL  SEED 

continue 

00,200  L-1,3 
"  3.1)  THEN 


THEN 


9916  NCRIC=0 
8810  CONTINUE 

m  ‘ 


IS(ll).LT.O)  60  TO  200 


C  SET 


ASH=30. 

CX1=X 


AN6L?  EQUAL  TO  15  DEGREES 
V^?EX( 


707 


708 


GO  TO  708 


CONTINUE 
LL  IN=0 


'A=CXi^CY2-Ytlt)tCX2*lY(i)-CYl)+X(I),(CYl-CY2) 


A'. 

(A. 


^2=(ixi- 


„  „„„  lOiA.LE. i'.Ol'GO  TO  250 

iMW0.10.250... 


.GE.rLO^i 


X[I 


SD3=(CX2-X[ I 
IF (IDASH.E 


**2+ 
**2+ 
D  AD 


tIcYMilli 

=8024503 


”2 

2 


B5 


USE  STATEMENT  ABOVE  TO  MINIMIZE  PERIMETER  OF  TRIANGLE 
USE  STATEMENT  BELOW  TO  MAXIMIZE  VERTEX  ANGLE  Of  TRIAN 


II  limviUkVi  VI  rw  j 

IF ( AO . GE . ASM]  GO  TO  250 


VfRTEX  ANGLE  OF 
7S0RT!S02’S0?) 


TRIANGLE 


CONTINU 

IG-VFRT 

kn=n5de 
IF  KN.E 


0,0) ' GO  TO  88i 5 


ilJtM 

,0  GO  TO  4 


& . 

LCC=LL-LP*100000 


00  40  K= 1 . KN  , , 

L INO=NOOE  t LCC ) / 100000 

ll=nodeTlccI 

LP^LL/IOOQOO 
LCC=LL-LP^  100000, 

xi=xTvertexTm, IS  1 
X2=X  VERTEX  M,I$  2 
X3=X l LINDT 


X4=X  iP  , 

SS.BHIWill 


I 

CALL  AR 


CALL  AREA(X1,: 
IFlA.LLO.Ogl, 


CALL  AREAIX2, X3, 1 

{Fific^tio  }!  t! 

NCR  C=3 
GO  TO  8810 
END  F 

tar,  , 

ip 

LCC=LL-Lr*100Q00 
00  45  K=1,KN  , , 

L IND=NODE  t  LCC 1/100000 
LL=NODEILCCI 
LP=LL7100Q06 
LCC=LL-LP*100000 
X1=X  LIND]  , 

!i;i  v|H:  1 

X4-X  I]  , 

Y^Y  VER?lx(n,IS 
Y3= Y  VERTEX  M, IS 
Y4=Y  I)  f 
CALL  AREA  XI. X2. 

ifTa.le.o.ooIH 

CALL  AREA(X1.X2, 

IfInc^S£7e6:i)  ti 

NCR IC=3 
GO  TO  8810 


wwu'k) 

THEN 


iiiraim 

feiiiiiiii 


K2.X3. Y1.Y2.Y3.A) 

So  to  is 

X2,X4,Y1,Y2,Y4,A 

W  45 


C  IF  T 


h[  M  if ?H^N  150  DEGI 
MIOASn.EO.O. ANO.AD.LT. -1.73)  TI 
IOASft-1 
GO  TO  9916 
ENOIF 
ASM= AD 
NU  =  I 

IFJIP.EO.NP.AND.NCRIC.EO.O)  THEN 
IFfNW.EO.O)  GO  TO  260 
NCRl C = 1 

j-jljQ 

GO  TO  8816 


IREES  FIND  MINIMUM  PERIMETER 
HEN 


B6 


ENPIF 

>n  contini.-; 

IFiNU.EQ.O)  GO  10  200 
15-HB+ i 

P1N0U  IS  THE  PARENT  TRIANGLE 
MNOW=n 

VERTEX  MS.l  =NU 

VERTEX  (IB. 2  =VERTEXIN.  ISll) 

VERTEX (MB, 3  =VERTEX(H.  IS  2 


m\m 


mki 


CALL  UPENIX (VERTEX fnB , 2  ,Y  VERTEX  MB, 

tt  m  Sltil  ;!  IP  i 

CALL  UHOVEf  X(Il)*X(l61+Xll3)]/3..  y(1 


xxnB=MB 

CALL  UPRNTllXXnB,  INTE  ) 

CALL AREA(X|I1),X(I2),X|I3),  , 
♦  ,  Y  II  ,Y  121  Y  13  ,A 

IFlA.LT.O. I  GO  to  296 
IQ=VERTEX  MB, 2) 

VERTEX  HB,2J=VERTEX(HB, 3 


tli*Y(I2)+Y(I3l)/3.  ] 


VERTEX  MB  3)=I0 

00.295.11:1,3  _ 


CALL.0R6eR(NG,P(2), 


(NG,P(i),P(2)) 

=NG 

(MG,P(2),P(3J) 

(NG,P(3),P(ll) 


1F0N=0  , 

MM ' 

ll=nodeTnw 

l  pr 

tell! 

LINO=NOQE  LC 
klANGLEIST 

lp=ll/iooqoq 


:N?f PAREN)=- 


IABS(SI0E(I1PAREN,  IPAREN) ) 


m  m 

i5B!&Io£fnB’ 


00000 
0  TO  311 


100000 
,KN  , 

to1 

I  WHICH  SHARES  NODE  NU  IS  f12 
P‘100000 


TO  310 


ftl) ) .EQ. I ABS ( SIOE ( f!2 , H3 ) ) )  GO  TO  350 


SIDE (MB, HI  = - I ABS 

pwriAg‘ 

IF(H2.NE.H  THEjl 
CALL  USETr  INTEV 


ANO.SIDEjflB,Hl).LT.O.)  GO  TO  320 


$(SIDE(fl2,N 


H23-H2  , 
Kl=-$IDE(f12,H3) 

K2=K1/106o 

K3=K1-K2M000 
GO  TO  320 
END  IF 

Iiz=n2  , 

J1=-SI0E(H2.H3) 


POT‘.’IFON=3’ 

THEN 


GO  TO  320 
320  CONTINUE 
300  CONTINUE 

310  CONTINUE 

311  CONTINUE 
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,IC) 


541 

542 


B , I 2 , 13, NU , NOLO 


1  +O.Oi*ZINT 

- NT 

JNT 
INI 


♦u.ui 

Mm 

-0.01*21* 


543 

544 


GO  TO  542 
CONTINUE 

JffMw 

00  543  Ic=l,fcc 

mk'v  ,| 

If  1 

IF  22.E0.ZT  4 
ZN=AHIN1  ZT  1 

zh=ahax|i2t 

+cSi^N|^?t  ftbf  2^1 2 , 1 3 .  NU ,  M°L0 ,  K  J ) 

CONTINUE 
IF jlFON.EQ.O) 

IFlIFQN. EQ. 21 

Mf1 
M 

60  TO  548 

IF?I^0N.EQ. 1)  THEN 
IF0N=IFON-l 
HN0U=nZ 


-ZT 

=ZT 

=ZT 

=ZT 


GO  TO  200 
THEN 


548 


60 

»  ?" 

200  ‘ 


fc=l  .3 

""  rtf 


n„. 

no, 


400 


_N0  IF 
CONTINU 
00.546  _ 

m 

l-VERTEX 

NW-I, 

547 
NUE 
305 
NUE 


(H.GT.HB) 

'  ^0 

-  A _ 7GT.0) 

HC0N=NCON 

NC0N=0 


TO  546 
TO  546 


GO  TO  400 


THEN 


B8 


00  45A  ]P=1, NP 
!  =  !-:' flflP) 


.  ORAU  EDGE 
KCONlNUflCO 
CALL  USE T( 
OPEN (6 1 . F I 
WRITE  f 6i . * 
DO  llil  U 
CALL  UMOVE 

mm 


. f'T.ZMAXj  ZHAX=Z[Z] 

.lt.zhin!  z«in=z! I ! 


N+l  =KCON( 1 ! 
COL} ) 

L£=  EDGECON’ 

]  NUMCON 
l.NUflCON  , 
(iCKCONfl)  0 

nfflitl  I; 


{ KCON ( I  + 1 ) ) , Y ( KCON ( I + 1 ) ) 


m!:^IABS(SI°E(J!K),)  SID£(J,K)=IABS(SID£ (J,K) ) 


«ia: 

11; , 
CALL  USET iCOLl 
CALL  UHOVEIXtl 
CALL  yPENlill? 

CALLEUslt(C0L6i 

S&iCrl.L 


GT.O)  GO  TO  660 


),X(I2),Y(I2) 


SRTEX 
RTEX 
RTEX 


I:B)  gS  II 


GO  TO  668 
662  CONTINUE 

668  MWh 


3  I=liSTART,HB 

35=1,3  ..  .. 


•GT.O)  THEN 


ii=k7iooo 
I2=K-I 1*1000 


. EO.  XfllN.  ANO.  X  12 
,E0.  XflAX . AND .  X  12 


.EO.YHIN 
.EO.YflAX 
.EQ.XHIN 
.  EO. XMAX 


b'J  IU  i«4 

END  IF 
CONTINUE 
MST  ART  =MB 
60  TO  385 


GO  TO  383 
GO  TO  383 


B9 


35i  (1?T  AR  T  - 1 
335  CONTINUE 
r=nsTART 
CALL  PUAIN 
«50  CONTINUE 
NUt1BR=lBL 
NP=NPT 


^51) 


CLOSE.. 

RETURN 
END 

SUBROUTINE .NODES] 


/ 


HEN 


NAV.NPMAX 


40 


P|pNT*.  MORE  THAN  2000  TRIANGLES  --  PROBLEM  TOO  BIG 

Ml  IF 
H=I  ,  , 

LL=NOQE ( I) 

LP=LL/ lOQyOO 
LCC=LL-LP*100000 

I2=LP*100000+L{c 
node(iT=Tlp+i7*iooooo+lcc 

KCC-LCl 

IFIKCC.EO.r  -- 
NODE  [ft  )=N0! 

NODE  NAV)= 

NODE ( NAV+1 
NAV=NAV+2 
RETURN 
END  IF 

h=lcc  , 

LL=NODE(KCC 


99),CX(3),CY(3),CZ(3), 
PX,PY,PZT,CL$ 
luu^IRO'/lOO)  ,KRD,  AVX,  AVY 


**2 

l),Y(P2).Y(I),Al 


(D.GT.DS)  GO  TO  200 


BIO 


DsrD 

P3-T 

CONTINUE 

CALL  ARE  A  ( X ( PI ! ,X ( F2 ) . X ( P3 ) . Y ( PI ) , Y ( F2 ) . Y ( P3 ) . A ) 
jPlA.LT.O.)  GO  TO  300 

P2=^3 

P3=0 

VERTEX  1,1  =P1 

VERTEX  1,2  =P2 

VERTEX  1.3  =P3 

N6=P1*1060+P2 

IPlP2.LT, PI)  N6=P2* 1000+P1 

SIDeT 1 , 1 1 =NG 

NG=P2*lQ0Q+P3 

|f(P3.lT  P2i  NG=P3  1000+P2 

S1NG=P1M000+P3 

l:R:i 

m  MB 


CALL  NODES  M  M  I 
END 


SUBROU  INE  DLAREAlAl. A2,A3,B1, 
0=S0RT  MAX  IA1-A2 |s*5+t01-6?  * 
♦  .  Al-Aol*  2+ f B1-B3J * *27 

A=A1  B2+A2*B3+A3*B1 -A1 *B3-A2^8 
IFJABS  A/D).LT,0.04*DJ  A-O. 
iFj^GT.-O.OOODOOl.AND.A.LT.G. 

END 


*2+ ( 82-B3 ) *  *2, 


B1-A3*B2 
.0000001)  A=0. 


SUBRJUTINE  ARE 

fuUftN 

END 


1001 . AND. A. LT. 0. 0000001 )  A=0. 


INTEGER  SIDE, VERTEX, h, fa, P3 

HiIMWeximii 


imm 

CYlJ  =Y  VERT 
CZtJ  4  VERT 
CONT  NUE 
RETURN 
END 


SUBROUTINE  ORDER  N.N.I 
N=M*10Q0+I 

IFII.LT.M  N=I *  1 000+n 


SI"* 

"8- « ■ » . '  2 . 1 = . » J 1 

common/stuff/zt  u]  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  , 

COfinON/COOR/NCCjt 999) , X ( 999  ,Y{999 ) ,Z(999  ,CX  3  ,CY  3  ,CZ  3  , 
+  VERTEXi^OOO . 3 } . SIDE  t  ^000 , 3 J , NP 

COMMON/, CRAP/GDEL,Mfl  , 

COHMON/BOX/IBNuJ,  54 )  ,  , 

8totos}8Kifat?S^fo£St”’tL*B“-PSI“-LP“<'1 

COMMON/WINDO/XUMlft ,  XUHAx  ,  YWfllN,  YHMAX 
CONNON/DAREA/XUDMI<l,  XUDMAX ,  YUDMIN,  YUDMAX ,  FSOFT 
C0MM0N/HVY/HEAVY,NUMH(500) 

COMMON/HHVY/HLET  52) 

CHARACTER  HLET*2 

CHARACTER  COCl*4.COL2'4,COL3,4.COL4’4 
CHARACTER  DEV*4 , fcOLOR *4 , CLRBCK k4 , PSI ZE * 1 
CHARACTER  2  LET152) 

CHARACTER  BLET *2  , 


Bll 


ioj! . !d\! . ;a\i 

;r\; 

M' 

:t\'. 

. u\, . . x\. . ,y 

z\ 

A\ ' 

s\  . 

,E  :.;F  \  .  6 

,H\ 

• '  \  • 

•j\’. 

;n  ...N  .  0  , 

P\’ 

'  Q\' 

R\  . 

u\  . ’v\’  ’u\’ 

TA  IFLOP  1/ 

’XV 

’  Y\ ’ 

’IVI 

Itfifiiie- 

IF?ABS(HH].LT.Q. 1 "HEAVY )  IB0LD=2 


mu 

IFIISCC. 


IFIISCC.ES.OT  REtCIRN,  , 

nwM  if  i  t  '  a 

£ALL_L0CTZZ,ZTT2)  ,2T(3) 


,X12,Y12,I1,I2,ISCC) 
,XHB,YHB,I1,I3,ISCC) 
,  XMB,  YMB, 12, 13, ISCC) 


100  ^ij|jgpZ,2T(2) ,21(31, XMB, YH6, 12, 13, ISCC 

200  Siitfh'i,hU,f‘ll"’W',!’U’,StCI 

IFIISCC. EO.O]  RETURN 
do  TO  450 


HH.EQ.1)  KB0LD=1 
763  JBOLD-l.KBOLO 


CALL  USETICOL  ) 
iniBOLO.GT.il  CALL  USETIC0L2) 

Y2=Y12 

ifl  ( XH-X2)  f*2+{Y(1-Y2}**2,LT. 

s4  ,o  ” 


FBL-O. 75*FS0FT 


Yl=  YtIB+Y 
A-  XM+X2 
Y3=  YM+Y2 
do  0  350 
300  Xl=  XH+X2 
Yl=  YM+Y2 
X3=  XM8+X 
Y3=  YfIB+Y; 
350  CONt  NUE 
XN=X 

c 

02=SQRT  ; 


mm  v 


IK 

XN0U=X2 

YN0W=Y2 

pili! 

CALL  UflOVEt X1+ 
CALL  UPENIXNOU 
GO  TO  400 


l;fo3)....Al  IS  ANGLE  1-2-3 

!  ttj 

>1**2 

SL-03) )/SL) 


0.01)  THEN 


. AND. XN0U-X2. 
.0R.XNQU-X1.N 
‘SINIANGJ.Y1- 
L*SIN(ANd) , YN 


, YNOW-FBL 


ANG=AUN2(YN0U-Y1,XN0W-X1) 

'rails)] 


TANA=R 

A1=ATA 


SUCH  THAT  THE  OISTANCE  BETWEEN  THE  PI  &  THE  CURVE 

S  L02-(X2-X1?**2+Vy2-Y1  )^2L  SET  °°EL:1  F0R  SHARP  c0NT0URS'- 

CH0RD=02*C0sfl . 5707963-A1 )**2 

OD2=SoIt  oIS1 il-5707963'A1^ i/(l.-COS(1.5707963-Al)) 

..  IFlOCH.LT  002-0.05]  THEN 


C  NEW  COORDINATES  FOR  POINT  1  (IF  NECESSARY) 


B12 


60  TO  378 
END  IF 

C  . . .SIMILAR  TRIANGLES. . 

iMqktTchokd’dz/toJ-choroi  ! 

AA’ATAN2( Y2  Y1  a2  XI) 

C  AND  YC'COOROINATES’OF "center  OF  CURVE  OF  RADIUS  RAO 
XC=X1 -R* SIN  AAj'SIG 
YC=YlfR  COS  AArSIG  ,,,  , 

RA,D=t i X 1-XC  **2+l Yl-YC )  2) 

C  HAKE  SURE  THAT  THE  CURVE  CROSSES  THE  COMMON  EDGE  OF  THE  TWO 
C  TRIANGLES.  .IF. NOT.  HOVE  THE  PT  POINT  1]  CLOSER  TO  POINT  2. 

Km 


imjai  aw*  11  closer  10 


837 


tmm 

FB=-2. *Yt 

FC-SXj^Xl+^’XC+YC'YC-RAD^.'SXrxC 


-X(Il)) 


.  IF(.F2..LT.Q.l  GO  JO.  834 


*GK*GB 


itS, 

iALLCQNSlS 

CALL  UARC(XC+FBL*SIN(ANG).YC-FBL*COS(ANG) .DEL  1 
CALL  UWHEREIXNOW.YNOWT 
ANG=ATAN2(Y3-Y2,X3-X2) 

isaaim;'051*11 

.XN6U.6T.XWMAX)  GO  TO  890 
(YNOW.LLYWHIN  OR.YNOW.GT.YWMAX)  GO  TO  890 

.  _  BANG=1Y3-Y2)/(X3-X2) 

IGN=5. 


!  Y1-6?*FBL*C0S 
,YU6.  *FBL  COS 


ftUPEN(X3+FBL*SIN 
AVY.EQ.O.OR.IBO 
OW.LT.XWMIN.OR. 


B13 


IF  BANG, GT ,f.  I  VSK'N--3Q. 

f  ALL  UH(iVE  I  f  XNOU+K  J + V3 ICN* F SOF T  >  / ^ .  I  YNna*  T'  1  / 2 .  ■ 
!BX=( ! IXNOU-XUniN] 7FSOFTJ/60. 

I B Y  = f ( IYNOU-YUHIN  /FSOFTi/60. '+2. 

IF  I6X.LT. 2 )  IBX=2 


I. 

00  893 


BX. GT.  7 
BV.LT.2 
~Y,GT.  5 


893 


!b2= 


IBX=71 

IBY=2 


Y®) 


XK J=K J  ,, 

CALL  USETf ’ INTE’ 
CALL  US£T(C0L2) 
KF=-26 


8M;pi 
=1 


GO  TO  890 


890 


ONTINUE 
ENOIF 
763  CONTINUE  , 

CALL  USET (COLl) 
RETURN 
CONTINUE 


450 


ANG=ATAN2(Y12-YH, X12 
KBOLD=IBOLD 
IF (MM. EO,  1 )  KB0LD=1 
00  784  J60L0=1.KB0LD 
CAI^L  USET  j  COL  J 1 


-XM) 


.GT.l)  CALL  USET ( C0L2 ) 


«Lg 


-  .  JOL.. 
CALL  UMOV 
CALL  UPEN 
764  CONTINUE 


Wcoul 

END 


COS(ANG) 1 /2. ) 


,CX(3),CT(3),CZ(3), 


fa>mRET,,R" 


pS|lll!iRi&E5;t5§IJ:SSMS'CLR!CK'RSI!£'LPE",!l 

COHNON^DAREA/xyon  ifl ,  XUONh ,  YUDfUN.  YUOflAX ,  FSOF  T 
CHARACTER  DEV*4^ 


ifTS*(S£S.cofE8R;‘«:^TH|;‘'RSIZ£'1 

—  -llT- 


.. OUTPUT  TO  SCREEN  OR  PLOTTER... 
DEV=OEVC  ,  , 

me® 

IF  OEV.EO.  HP4 
IF  DEV.EO.  HP4 
CALL  USTART 
CALL  UERASE 

CALL  UBAr; _ 

CALL  USET  MEOIU 
IFlDEV.EO ,  IBMH  , 

CALL  USET  ’PERCENT 


T_  USE 
IFIN.EO.I, 
IF1H.EQ.2 
CALL  -■ 
CALL  USETT 


ALL  USET(COLOR) 

CALL'UDAREAC 10. 
|ALLUOAREAiiq 


.  PCX, 75. *PCX, 10. *PCY, 75.  PCY 
.  "PCX, 75 . *  PCX .10. ‘PCY.75. *  PC  Y ) 
. , YUMAX) 


X , YUMIN 
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=5FT=i . 

F!n. EG. 21  P$FT=  67 
CALL  UPSET!  VERT  ,40’PSFT 
CAi  L  UPSET !  HORI  ,20.*PSFT j 
|||yQ.2!  CALL  uSET!  HARO  ] 

END  IF  , 

IFTh.E0.3)  THEN 
..OUTPUT  TO  PRINTER 
DEV=DEVC  ,  , 

fflwmjoM.) 

CALL  HIe^PERCEMT’) 

CALL  UD ARE A | 10 . .75. . 10. .50. ] 

CALL  UWINDOtXWrilN.XHMAX, YMhIN, YUMAX) 
CALL  U$ET( ,SOFT* )' 

CALL  UPSET  f VERT  ,30. 

CALL  UPSET!  HORI  ,24. 

RETURN 
END  IF 
RETURN 
END 

SUBROUTINE  PEND 
CALL  UF LUSH 
CALL  UENO 
RETURN 
END 


SUBROUTINE  INF f NC ) 

CHARACTER  FLNE*14.B0RD*2, HONORS, TYP*4 
CALL  UHOVE  810.,5i0.) 

X-775. 

Y=750. 


NC.6T.26 


_ ,  X=81Q. 

NC.LT.26)  Y=NC*25. +100. 

.  Y.LT.550. )  Y=550. 

CALL  UPRN^lf’CONTOURSV , 'TEXT' ■ 


RETURN 
END 

SUBROUTIN' 
IN  _  ' 
cotinoN 


BROUTINE  PR0FI 
iTEGER  SIDE ,  VER 

nnoN/cooR/kci 
♦  ,  yERT 
COHHON/DVICE/DEV. 
COflflON/TRI/flBL.HB 


conriON/wiNDO/xwniN.  xumax  ,  ywmin.  y 


Jznx,  zmn,  nn,  xhinx,  xiiaxx,  YfiiNY,  ypiaxy?  sn,  rdask) 
aCX(3),CY(3),CZ(3), 

, PSIZE , LPEN ( 8 ) 


CHARACT 
CHARACT 

CHARACT..  . . 

CHARACTER  RDASK‘1, ILX’l 
I^P*’  i  SH  ' 

IP:T  ’’ 

CALL  CLEAR 
PRINT’ 

PRINT,,.. 

READ  Ul 
IFlIASK.f 


VV'-l  4  i  VvLt  H  I  VULv  4lVVL*  H 

DEV'A.tOLOR  i.CLRBCK1*. PSIZE  1 
IP* 1 , 5 ASK* 1 . F ILEN* 12 , DESC  *80, S 
MNMix*40,PhLE*20 


SLASH*1. IPP*2.YLABEL*40 


25a)  ARE  PROFILE  LINE  ENDPOINTS  TO  BE  LOCATED  WITH’ 

,  THE  CURSORIC)  OR  FROM  A  NAHED  FILEIF)  ?tDEF=C I  ! 


$K 


JASK=’C’ 
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:r'  JASF.EO.  'c  !  GO  TO  7744 
PRINT',  ENTER  NAME  OF  FI!F 
PRINT.  I a\  OF  ENTEP 
READ  •  A20)  ,  PF 1 1  E 
IFlPFILE.E'j.  CON T  THEN 


INT’.'ENfRR.NAHE  Or  FM  WITH  Y1.Y1.X2.Y2  ' 

INI.  IA\  .  OF  ENTEP  CON  IF  FROM  KEYBOARD  > 


PRINT  TA^J  ,  TYPE  in  X1.Y1.X2.Y2  ) 
READ* .Xl.Xz.Yl.Y2 
GO  TO  7 1U 

ENDIF  ,  , 

OPEN | 33. FIL£=PFILE,ERR-3334, STATUS^  OLD 
READ133.  I  Xl.Yl,xi,Y2 

xi=  xi-$u6mAci 

xi=  Txi-xniNxj/SfD+ioo. 

8=  p2-XOX]^H)  +  100. 

Yl=  Y1-SU8YJ*FACY 
Yl  =  (Yl-YMNY]/StthlOO. 

Y2=  Y2-SUBY I *FACY, 

Y2=  IY2-YMNY]/SM)  +  1Q0. 

CALL  UCRCLE  Xl,Yl,5.J 
CALL  UCRCLE  X2.Y2.5. 

GO  TO  7755 


7744  CALL  CLEAR 

print*,; 

PRINT,  A\ 


SSl'IMa 

^RINT^t2.25c) 

PBT*IT«  ’  pntlTt 


I NY  /SH  +100. 

if "iiMW  H  J®‘?£i  W 1,0  ura ' 


PRINT*,’  POSITION  CURSOR  AT  1ST  END, OF  PROF.  LINE  S  HIT  RETURN’ 
callTuhovB(io§  Soo  )AND  6  KEYS  10  FINE  PO»nION  CURSOR. 

ifIilleo^I’caiIl^loseskxi.yi) 

^llT^,[2.25d)’ 

PRINT*.’  POSITION  CURSOR  AT  2ND  END  OF  PROF.  LINE  &  HIT  RETURN’ 
CALL  U§ET 
CALL  UCRCL 


m  n 

CALL  UB 


cleTxI ! Y1 , 5. ] 


CALL  KURSI 
IFTILX.EO, 


2  Y2  IP) 

’CAlL  CLOSEST (X2.Y2) 


ALL  CLEAR 

^X|ltEQ.X2.A^.YJ.E0.Y2)  GO  TO  7902 


CALL  UCRCLE I 
CALL  UFLUSH 


F ( A§S ( X2-X 1 } 
FIX2.LT. Xl! 


W ,H£" 


GO  TO  123 
ENPIE 


IF! ABS( X2-X1J.LT.0. 1. AND. Y2.LT.Y1)  THEN 

X2=TX 

Y2=EY 

X 1  =  TX 

Y1=TY 


123  CONTINUE 
XB1=X1 
YB1 =Y1 
XE1=X2 

Will  1=1,10*1,2 


)  GO  TO  335 


TFjXBl.LE.U.O.AND.XEl.LE.O.O)  GO  T 

IH  8:51  m\m\\ 

IF piBl.GT.O.I  THEN 
X6i=Xl+F*  X2-X1J 
YB1=Y1+F’(Y2-Y1I 
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IFfXEl.GT.n.O'  THEN 

sFuxi-P’Ir'y-Ml 

YE1-Y2-F’<Y2-Y1) 

ENDIF 

1.M5L  WILL  INCLUDE  POINTS 


%  m,  iw 

CALL  CORCORl J) 


IN  NON-CONTOURED  AREA 


ifVxbi.gt.o. 

CALL  AREA  CX 
CALL  AREA  CX 


iftki-. 


ll%2 

2|.CX(3 


, YB1 , A 1 

-y.gi.A2 


, XBl.CY 

nm 

m 

IFlXEl.GT,  .  ,  . 

'  "  . . .  XEl.CY 

XEl.CY 


II  \  A 

CALL 

CALL 


HEN 

,CX 

,CX 


5  m  i  ;(!  Ml:!ij 

1  iXEl.CY!3j,CY  li.YEl,A3) 


vnww  ni.iiniv'i.ki  \  j.  /  •  v  i  iw;  ■  ■  w  X 

m  ^  3 M  IHB 

xe1=-xe i 


YE1--YE1 
END  F 
ENOIF 

334  CONTINUE 
333  CONTINUE 

335  CONTINUE 
Xls-XBl 
Y1=-YB1 
X2=-XE1 
Y2=- YEl 

CALL  USET  C0L1) 

CALL  UCRCLEIXl, Yl,5. 


5. ) ,ABS( (Y2-Y1 )/5. ) ) 


IN. 


100  CONTINUE 
ZMIN=99 . E30 
ZnAX=-ZnTN 

10  Jsl.NJ 
.  J=-999 

■IND  TRIANGLE  POINT  IS 
CALL  UALPHA 
DO  200  IG=II,(1B+II-1 
KK=0 

IFTI6.GT.HB)  KK=flB 
I~ IG-KK 

IF  I.LT.HBL)  GO  TO  200 
CALL  CORCORl II 
CALL  AREA  CX 
CALL  AREA  CX 
CALL  AREA  CX.„, 

IFIAl.LE.O. AND. A2.LE.6.AN0. A3.LE 
CALL. UALPHA 


END 

200  CON...... 

210  CONTINUE 

88  m*m 

PRINT,  Af 
READ  !A1 
IFCJASK. E 


orof i  1 
or 


e  !D)rawn,on  this  plot,  sent  to  ’ 
(Bjoth  ?  (Enter  D,F  or  B)  DEF=D  > 


<  iii  n  vn  ,k«i  .  .  1  JASK= ' 0 ’  .  . 

F 1 J ASK .  NE .  ’ 0’ .  AND . JASK . NE . ’ F ’ . AND. JASK .NE.  B 

IF( JASK. EQ.’F’. OR. JASK. EQ.’B’l  THEN 

^NTCl^)' 

12.271  ENT^R  A  FILE  NAME  f DEF=D . PRO )  >  ’ 


GO  TO  88 
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REAP  IA!2) \FTiFN 
IP  I FILEN.  EO.  1  FILEN--  P.PRO 
OPEN!  14.FILE=FILEN. ERR;8  / 1 
60  TP  R9 

37  PRINT’ fFILEN.  ’  ALREADY  USED-  HIT  RETURN,  THEN  TRY  A  NEW  NAME’ 
CALL  U|ELL 
CALL  UPAUSE 
‘0  TO  91 

'  ‘LEA.. 

NTER  A  ONE  LINE  DESCRIPTION  (OR  RETURN) 

C 


89  CALL  CL? 
PRINT, , 
READ  ’A 
CALL  CL* 

DDTfelT  T 


(2.29)  ENTER  SECTION  LABEL  (LIKE  A  FOR  A-A)  > 


iFfYLAfiE^EO 
+  YLABEL-  ELE 

vx=$n*Y  * 

VY=Sh* 

URITEt 
VX=SM* 

vy =sn* 

WRITE! 

CALL  UT10VE 
CALL  UPRNT 
CALL  UHOVE 
CALL  UPRNT 
CALL  UM 
CALL  UP 


IP, SLASH 

TER  Y-AXIS  LABEL  (DEFAULT=ELEVATION) 
1  YLABEL 


_  VATION’ 
X1-1Q0. ) +XHINX 
Yl-iOO. j+YhINY 
4*1  VX.VY 
xi-lOO.  J+XfllNX 
10.  J+YPIINY 
VX.VY 
1-iO, ,Y1  , 
IPP ,'tEX  ’ 
2+26, . Y2 
IPP,  tpxt 


') 


40) ’ )  YLABEL 
XT(I),YT(I),2T(I) 


1**  •  )  “  • 

14 

K.EQ.’F’)  GO  TO  7902 
LEAR 

PRINT  12,  ZflN.ZflX.ZMN.ZHAX 

34  PR  NT*,  2  Jl)  ENTFR  YOUR  MIN  AND  MAX  PROFILE  ELEVATIONS. 


».  'ilF' 

DEF=MIN  AflD  MAX  ELEV  OF  ALL  POINTS  ON  PLOTS) 

rs  M  i/ni  ' 

ZZf1IN=ZHN 
ZZHAX=ZHX 
THE' 


...RAX, 

INMAX. 

JNMAX.F 


K¥a;OT&ax 


[|h,^Afu^’OELETET 


,  ZZMAX 


7902  CONTINUE, 

'*  E0.2)  THEN 

.If' 

PST  ART ( 2 ) 

UERASE 

CAL  75.  *PCX,  10.  *PCY, 75.  *PCY) 

THEN 


ON  ( ZHX,  ZMN,  0,  XMINX,  XMAXX ,  YMINY,  YflAXY,  SR) 


&Lxlofi^-*N! 

CALL  UPRNT 
CALL  UftQV* 

CALL  UPRN 
IFllP.NE. 


Y1.E0.Y2)  RETURN 
PPi’tEXr) 


CALL  UPRNTlt ’ *\’ . ’TEXT’ ) 
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<JJ  CONTINUE 
XP= XNOW 
YP= YNOU 
COLOR-CLR 
CALL  USET(COLOR) 
RETURN 
ENDIF 

callRus^ticolor) 

RETURN 

END 


REJURN 


OLDfISH’12 


(2.23)  WANT  TO  DRAW  A  PROFILE  ?  (Y/N)  > 


CALL  CLEAR 

print, 

read  Mai  ..... 

CONTINUE  ,  ,, 

IFUUjE^lh  ThIn  '  ,  . 

102  PRINT1,  M2. 40)  ENTER  NAME  OF  'MESH  FILE  ( OEF  =  .OLDHSH.  ’  I’ 


1456 


THEN 
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PRINT  I  A'\ 
READ  l A12 

Umil 


OR  0  TO  our 


.ntlESH 
.  .0  T  RETURN 


FIflESH=OLD«SH 
OPEN  1 81  .FILE=F  IflESH,  STATUS:  OLD' .  ERR= 1  QOl 
QL0t!$H:F|f1ESH 

*  7 HE  NOT  FOUND.  TRY  AGAIN’ 

SUBX.SUBY 
FACX.FACY 

Xf1IN,XNAX,YNIN,YHAX 
ZN^piSi.XtlAXX 

B,nBL 

Xjl),Y(I),Z(I).NCCm 
?VERTEX(I,J), J=l,3) 


ANT  TO  REDRAW  THE  CONTOURS  ?  (Y/N)  > 
.ROASK 


m: 

CALL  CLEAR 

PDTUT*  tMT 


uaAt  the  volume  under  the  MESH  ?  (Y/N)  > 
.'y^then 

PRINT*. "HINUinun  ELEVATION:  \ZI1IN 

_ 

m;jt!gi:;R!fiiEFELE:o.o 

AlOrl  REFEL 


NtER  REFERENCE  ’’BASE”  ELEVATION  (0EF=0.0)> 


1931 


CX(2) ,CX(3),CY( 1) ,CY(2) ,CY(3) ,A) 
)/3. 

vbLl* (ZMAX-ZHIN)/ (V0L2-V0L1 ) J+ZMIN 
CALL  UDAREAl 10. *PCX, 75. ‘PCX, 10. *PCY, 75. ‘PCYl 

mM' 

CALL  ?s|aRI(1) 

O0A^a!^R^cI,75.  *PCX,  10.  *PCY,  75.  *PCY) 

CALL  U 

mi., 

CALL  UlIOV! 

CALL  UPEN 
CALL  UMOV 
CALL  UP£N 
CONTINI 
CALL  U 
CALL  UP 
CALL  UP  _ 

CALL  XYAX 

MW 


COLA 

C0L2 


XHINX.XHAXX. YniNY.YHAXY.SH.il 
THEN 


CALL  RDCON 1 ZMAX . ZMIN , 1, XHINX. XMAXX, YHINY, YHAXY.SM! 
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HU 

READ 
READ 
READ 
READ 

READ,ul,  , 

1928  WfllJi 
DO  1'"'*  - 

1929  READ 


a  . '  FFF 
§  .  FFF 
81,’  FFF 
§  FFF 
81,’!  FFF 
l.NF 
, . /  FFF 

'  T  i  Ai.  —  _  > 


THEN 

NIN’ XMAXX, YMINY, YfiAXY,  SH,  RDASK ) 


£ALL  Vjuit 

ofSff 

fill  Bf, 

CALL  UPRNT1 

Fall  B!i 
SfM! 
CALL 


,VERT,, 30. ‘PSFTl 

,hoRi;.i[«psft) 


[vi«'  , 

'WL  M 


JHBMS&i-C 

*,'|V»»V^«£fa£,v,„L,.  .,1I0I 


cal' 

M'JIWm”’  T0  "m  mnLC !  <»"' 


CmK  1W,  XfliNX,  XdAXX,  YHINY,  YMAXY,  $fl) 


CynHON/TRj7n5L  jj8 

““T'  ^  ’  •",m-as 
iszisy^'M *®  iMEMC0" 


^E**8,TITTLEU8)«i.PASK’1 
%IP*1,THPLT*20 


.■WBCL-I  I.8L«f  11 1  - 1 
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+  !  I  • :  s>\ ! . ;  o\ ; . ;  r\ ; . ;  s\ ; . ;  t  \ ; ,  ;  v\ ; . 

*  ,y\. v\  .  z\  .  a\  .  S\  .  c\  .  D\  , 


SET^A^ktob  6  ’ lPL^^P^AN’ ' ’ »AGE’ * ’ 


IFlABSIXUDMIN-lG ,’PCX) ,GT.O. 001 . AND. AB$(XUDMAX-75. ’PCX] 
♦  .GT.O.pOl!  I FLAG= 1 


^ALLFupsef?; 

M 


^RTf^Q^Psfl’FSOFT 
HORI  ,20.  PSFT  FSOFT 


CLEAR 

=9.E30 


znAX=-2niN^ 

AND.X(I).LT.XUHAX.AND. YED.GT. YUHIN.AND. 
+  Y  Ij.LT.YUMAX)  THEN 
IF  2  I  .GT.ZMAX]  ZNAX=ZlI) 


IFIZ  I  .LT.ZMIN 
ENDIF 

2307  CONTINUE 
CALL  CLEAR 
HEAVY=0, 

CONVAL-  . 


ZMIN=Z(I 


B^GliliiNG  toUioUR.^CON^pUR  INTERVAL.  HEAVY  LJNE’ 
PRINT*,  UPTIPLE,  AND  NO.  OF  CONTOURS (0  MEANS  SPAN  RANGE)/ 

/OR.’IT  ,F6R  NO  CONTOURING  (DEF=AUTO)  (H  FOR  HELP)  >  ’ 

DC  Ah  MAAftl*  /'ftAlUAI 


a  i  ivviunkiuv,  ,ny 

IF ( CONVAL. EQ.’HV 
CALL  CLEAR 


.OR. CONVAL . EO .  ’ N’ )  CONV AL= ’ 990999 ,1,0,1’ 
THEN 


FOUR  NUMBERS  ARE  TO  BE  ENTEREOtOR  PRESS  RETURN  FOR  AUTO)’ 


THE  2ND  IS  A  +  OR 
ETURN  FOR  MORE  INFO)  > 


W’ W  f  KaSK  Rl fpRKnETUl?N‘F5R*H0RE /NFO)  >  ' 

PRINT’,  ’fp**  (SCREEN).  THE  4TH  IS  THE  NO.  OF  LINES  TO  BEV 

*[F  DRfUN  m  flEANS  SPAN  RANGE  OF  DATA).  PRESS  RETURN  FOR  MORE)  ’ 
READ  ’(Air,  A  ASK 

PRINT’.fl.E.  SO^.^.iO.S  MEANS  START  AT  ELEV  5020  AND  DRAW  THE’ 
PRINT  ,  5020,5018. 5016, 50U. 5012,  AND  5010  CONTOURS . UITH  THEV 


THE  3RD  CAUS 
OF  THIS  VALU 
COpOR  (SCREE 


H  ARE  MULTIPLES. OF’ 
TER)  OR  ANOTHER’ 
f  LINES  TO  BEV 


ONVAL.EO. ’ 


ZDMZMAX- 


NC«=0 

GO  TO  7503 


L0G?L0G-1. 


’ip;0 

. 0 * *LOG ) 

nt  %rinm,zm 


IFIHEAVY.LT. zintT  the 


^»:EPoR;fp!Toz*ZINT'HEAVY’NCn 

IF (HEAVY. ^T. ZINT)  THEN 

Hf^VY, LINE  MULTIPLE  MUST  BE  GREATER  THAN  INTERVAL’ 
PRINT  MA\)  ,’  PRESS  RETURN  TO  CONTINUE  )  ' 
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READ  (All  . AASK 
GO  TO  2302 
END  IF 

GO  TO  2303 

2302  CLOSE  14?. STATUS:  DELETE 
GO  TO  2301 

2303  CONTINUE 

CL0SEU3, STATUS:  DELETE  ) 
7503  CONTINUE 
CALL  CLEAR  , 

PRINT  V(A\  ’ . 


,STATUS:'DELETE’) 


PRINT  t(A\)  . 

'  12. J7aT  UANT  TO  SEE  ELEVATIONS  PLOTTED  ?  (Y/N|  )  ’ 
READ  (All  , HELV 
CALL  CLEAR 

DDT  kiT  “Till* 


5 

READ  ’  A1 

pr^t^L£a 
♦^2.18)  LI 
PRINT/  A 
READ  ’(Al 
CALL  CLEAI 
PRINT*/"; 

DOTkIT  »  i  *' 


PgljJTS  f?/9)  E 


.LfcAVE  X-Y  AXES  ANNOTATIONS  OFF  (Y/N)  ?  >  ’ 

,  AASK 

AyE  LOCATION  JIC’ ’ $  AT  DATA  PJS.  OfF  TOO  ?  (Y/N)’ 
J.  (ENTER  Av’  TO  INCLUDE  GRID  POINTS  ALSO)  ) 


}9)  ENTE?  TITLE  OF  PLOT--! IF  NONE,  HIT  RETURN)’ 


CALL  UCHP 
DO  1923  I 

!F  NCC(l) 
FlNCCtl] 
NDIF  , 
FINCC  I) 
FINCCllj 
ALL  UHOV 
ALL  UPEN 
CALL  unovi 
CALL  UPEN 
IFTHELV.Ei 
CALL  UHOVI 
CALL  UPRN 
PNOIF 


000)  GO  To 


.GE. 1000 

eftffil 


ALL  USE!  COLA) 

r1 


I  §  5 SR 


SilSIllI' 


6^FS0q,Y(I)+6.*FS0FT) 


-EO.U  THEN 


\  LH'J  .  V*  v  ■  X  /  IIIL.I1 

CALL  UUINDO  1..1OA0..0 

ca^.udareaIio/pcx/5 


,780.) 

*PCX,10. 


*PC Y , 75 . *PCY) 


CALL  USET(COLl) 

FjNCH.GT.S2. AND. HEAVY. EQ.O. 
FINCH. GT. 52. AND. HEAVY.NE. 0,1 

F  NCH  GT.500)  NCH=500 
Nc~2 

777  ClTnC)=CL(1)*2INT*(NC-1] 

IE  ZINT.GT.fi. Q. AND. CL  NC.GT. 


IF  ZINT.GT.fi. 0. AND. CL  NC 
IF  ZINT.LT. 0.0, AND. CLINC 
IFlNC-l.EO.NCH  GO  TO  77 


NCH=52 

NCH=52*HEAVY/ZINT 

NCH=52 

NCH=52*HEAVY/ZINT 


GO  TO  778 


LINCi.LT. ZHIN)  60  TO  778 


GO  TO  777 
778  NC:NC-1 


ra« 

IsSliBxft 


HEAVY)  THEN 


NUHH  f I ]  =  J J 

EonIinue 

DO  760/  1=1. NC  , 
IFlNUHHlIl.NE.O!  GO  TO  159 
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7607  CONTINUE 
HEAVY =0 . 
ENOIF 

159  CALL  CLEA 
PRINT, Ha 


(2^15)  ENTER  TEMPLATE  FILE  NAME  (DEF-NONE!  ) 
THEN 


CALL  UCHPE _ _ 

open! 12, file-tmplt, Status-  old’ ,err=i59) 


JPN=rl 
49  REAP (12 
X7s 
Y7= 


X7-§U8X( 
UBYj* 


=59 3  X7,Y7,IPN 
FACX 
‘  Y 


CALL  UUHEREfX77, Y77) 
AN6=ATAN2(Y7-Y7).X7-X77 


CALL  UMOVE. 

isM-W® 


iH 


N.GE.OlCALL, 

UM0VEu7,Y7) 


FSOFT'SINl 


*FSOFT‘CO: 


CALL  UMOVE  X7-1‘F§0FI*§ 
"  UPENIX7-1.* - 


IF  f IFLAG . EQ , 1 )  THEN 

CALL  IJUINOO  XllflIN,  XUMAX .  YUMIN.  YUMAXj 

CA|A  UDAREAfXUONIN, XUDMAX , YUDAlN, YUOMAX 


*,J.AND.IU. 
SEtlCOLlj 
,(81,*  NUMCON 
DO  7744  1=1, NUMCON 

8Wl«L*.Ul.M|X2,Y2 

>.! 


NINY,YHAXY,SH,l) 


CALL  UM&VE. .... 

CALL  UPENIX2.1 
7744  CONTINUE 
CALL^EAR 

+  ' JflurL^E  BOUNDARY  OF  CONTOURED  REGION  UNDRAWN  ?  (Y/N)  > 
READ  (All  ,BASK 


.  MN1 ,  MN2 ,  MN3 ,  MN4 .  MN5 ,  HN6 
,ZTU) 


543  CONTINU 
GO  TO  1548 


MN3,MN4,MN5,MN6,K3l 
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1550  CONTINUE 
;54<)  CONTINUE 
CLOSE  1 81 1 
CALL  CLE A' 
PRINT, ;(A 


FO!ii:W 


UANT  TO  ORAU  THE  TRIANGLE  flESH  ?  (Y/N!  ) 


PRiNT^A^)’  ’  ERASE  SCREEN  ?  (Y\N)  ) 

MM  'HKthen 

CALL  UE&ASE, 

CAj^UBACKGlCLRBdO 

DO  510  I=nBL,HB 
II=I 

CALL  CORCORdI  ,  ,  , 

CALL  UPlOVEfCXf . J.CYC 
CALL  OPEN  CX?2,  ,tY  2. 

CALL  UPEN  CX  3  ,CY  3 

WP  bllj’CY1 

ENDIF 

. . DRAWING  TITLE  AND  LEGEND... 

IF  (in.EQ.l.QR.nn.EO.21  THEN  , 

CALL  UUlNDOU.,lO?Q.,0.,780.  _ 


.DRAWING  TITLE  AND  LEGEND. 
.QR.nnjQ.2l  THEN  . 


[Ifl.’^rCXj  lu6!  vPx.O.  *PCY,  100.  *PCY) 

,V^|ffSo.;PSFTj 

CALL  f)QAftEA[lQ.’PCXvlQO.*PCX,0.*PC 


IF  IW.EQ.I  CALL  OQAREATIu.  PCX. 100. *PCX,0.  PCY. 100. * PCY ) 
IF  IW.EQ.l  CALL  UPSET?  VERT,, 3ft. 

!F  IW.EQ.I  CALL  UPSET  ’HORI’, 15.  , 

Ftnp.E0.2.AND.IW.E0.0)  CALL  tlSEt  ’HARD  ) 

mtir 

CALL  UPSET (  HORI  ! 15. j 

!|K«.) . . 


0. ,0. ,67. 1 


BIft  f 

CALL  UPRN 
CALL  CLEA 
PRINT  y[A 


,;H8Al)’)  (TITTLEtl ) ,  1=1, 48) 


‘kIH  if ?  ,m ' ' 

IF^ASKJQ.’Y’)  GO  TO  162 

IPOW=NC 

KPOU=Q 

DO, 8900, 1=1.  NC  , 
IF(NUnH(I).CT.O.)  KPOW=KPOW+1 
8900  CONTINUE 

IF(KPOW.GT.O)  I POW=KPOW 
CALL  INF(IPOW) 

PX=775 . 

YF=700. 

IFIIPOW.lt. 26  YF=IP0W*25. +50. 
IFIYF.LT. 500 . J  YF=500. 


call  unovE 

CALL  UPRNT1 
CALL  UPRNT1 


Jojio 


C?1,8C,INC 


f  ffE 

YF= YNG 
PX=XNG 
ENDIF  ,  , 

CALL^UliJvl(PX,YF) 


G.O. AND. HEAVY. NE.G. )  GO  TO  500 
THEN 


WilffiKIklWixT’ 


IfVaCL jG?T^O. OOOi . AhSIaCL . LT . 0.0001)  fCL=0 
IF(A8SlACLl.LT.100.)  call  UPRNTlfACL,  real 
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TFlAR^fAH  I.GE.1QQ. I  CALL  UPRNT1 ! ACL.  IN7E  1 

iff rc.eo. n  then  , 

CALL  UPRNT1  V  .  TEXT  > 

CALL  UUHERE ( XN6, YNG7 
ENOIF 


,CY(3),CZ(3), 


J=0 

DO  100  1=1, Nf 
IF  XX.GT.X 
IF  YY. 

D=  XX- 
(D._ 

IN=D 


100  '  NUE  , 
NE.O) 


I  GO  T 
GO  T 
>+(YY- 
THEN 


100 

1C 


THEN 


RAUC0N  (IPCELL) 

it  - 


COHtlON/CL 

IBiiirfpiMiiir"'" 

*P  /NUtillOOO) 


Yj999|,Z(999),CX(3),CY(3),CZ(3), 
?fACt^PX,PY,PZT,CLS 


CELLS 


COMMON 

connr 

comnu., - 

common/nodel/.n 


„ _  ,  i  nviu  xxjuw  ;  j  \  iuu 

mmlnimhm.m'njm.m, 

COMNON/DRAUS/MBLNlim  100  J ,  NUfjBR 


,  IRD( 100) ,KR0, AVX,AVY 


ivJOO)  NAV.NPMAX  .  , 

■  t  PCc6[f  c6l  aL0R  1 CLR8CK  - PSI ZE  ’ LPEN  f  8 } 

Au,  SUBMUBY 


_ _ ELL  ( 1 

DIMENSION  HARK1999 
CHARACTER  HLET k2 

TUADATTCD  AOr*  T  AA 


0) ,XTM(50),YTM(50) ,ZR(2) , I AKSHN ( 50 ) 


1!  ®»,C0L3,,C0U. 


|  Ifc 

character 

CHARACTER 

CHARACTER  .^S^CjEMi ALABEK 99 


1,PASK*1 

a , tnPLT*20iPQT* 1 
XY*12,Fr1ME*14,F0RCFILE  20 
LSH* 1 , 8L ABEL* i 1 , BLABl 1 1 J *  1 

”’H  ”’I  ”’L  ” 

x\  !  Y\  .  z\  ? 


;  'fv”fv 
*  ;r  ;:;n  ; 

+  ’U\’.’V\’ 


u\’ 
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HAT  A  . 

Format 
OPEN!  8 
00  434  1=1 


COLRMELL  ,’RED  . ’GREE  .  BLUE  ,  C> 
TI3F16.fr.Ili.Al.AlO) 

8,FIL£=  MCINFO. OAT  1 


OAN'.'flAGE’.'UHIT',' 


DO  534  1= 
MARK! I )=0 


,.50 
.=6 

=1.999 


300 


250 


200 


7781 


7779 


READ  8, 
READ'* 

READ 
CIOS... 
NPMAX=0 
JNPLT=1 
j^ONTINUE 

Zt1IN=99.  E30 
Zf1AX=-99.£30 
XHIN=99.E30 
XflAX=-99.E30 
YPI1N=99.E30 
YMAX;-9$,E30 

or  ir»lA  " 


Air )  ip 
FACX.FACY 
!  SUBX , SUB Y 


READ(9  ' 

DCADfO  * 


READ 

R"‘ 


rIao 

IP?.. 

IP, 

IP 

IF 


xX(l 


9!* 

9,1 


Ai4 )’) 
FACX, 


fiwe 

FACY 
" BY 
X,_. 
xniN= 


,E«Mf 
.LT.xniN  : 
.GT.XHAX  XHAX=X 
.LT.YHIN  YMIN=Y 
.GT.YNAX  YI1AX=Y 

•\m  szz( 

a  1 

Mf2°° 

Mi250 


I ,  Z  ( I ) ,  NCC  ( I ) ,  ASPACE ,  ALA6EL  ( I ) 


999  NODES..  IGNORING  EXTRA  NODES. .HIT  RETURN’ 


7789 


PR 

cal 

cal 


*u 


PCELL=999 
Ell =999 


i  (IODES=  ’  ,NP, 1  . PRESS  RETURN’ 


ACX  +SUgX  ,  X«AX= 

'  . .  YHAX=  , 

TO  CHOOSfc 


T°f7779. 


Xtli.  , 
YflAX/FA1 


THE  ABOVE  ARE  USED  TO  CHOOSE  A  SCALE 
.WANT  TO  CHANGE  THESE  VALUES  ?  (Y/NT 
ER  NA{IE  OF  FILE  GIVING  XMIN, XPIAX  YftlN, 


♦  SUBX 

Wl 


THE’ 


YflAX  ’ 


,  f  ILEXY 


.  W:- 

ENOIP 


7779 


GO  TO 
THEN 
N  AND  XMAX 


ENTER  NEU  XU 
ENTER  NEU  YfllN  AND  YflAX 


OPEN  44, F 
READU4  * 
CLOSE 
XMIN= 
XMAX= 

YfllN 


K 


_=FILEXY,STATUS=’OLD’ ,ERR=7789) 
XMIN,  XPIAX,  YfllN,  YMAX 


U\ 

XHIN-SUBX  ‘FACX 

mXjfYHAX-SUBYj‘FACY 

SX=IXHAX-XHIN)/900. 

SY=  1 YMAX-YMIN ) /640 . 
SM=SX 

ji  w  sm:sy 

. 

Hfl.EO.l  ft® 


F 

F  I 


JNPL 

Finn. 

Finn. 


IF 
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IFi.JNPLT.EO.lt  CALL  UEPASc 
CALL  UBACKG! CLP8CK  i 
Xf1INX  =  Xf1IN 
V(1INY=YM1N 
XflAX  X=Xf1AX 

yhaxy=yHa 

..+100. 

snl+ioo. 

xQi=ioo: 

II 


1077  CoiJlNUE 

"  J‘  ' 

EN^F 


. JNPLT.EO. 1) 
ALL_CLEAR 


THEN 


00  1914  1= 
CALL  UtlOV' 
CALL  OPEN 
CALL  UflOV 
CALL  UPEN 
1914  CONTINUE 
CALL  CLEAR 
NFALT=0 


jT^T.EQ.O)  THEN 

-  «p 


PRINT *,  (2,P9  DO  YOU  WANT  TO  DESCRIBE  AN  INTERNAL  CONVEX  AREA’ 
PRINT, '(A\  Y  T0  CONI  AIN  ALL  CONTOURING  ?  (Y  OR  N  )  1 
READ(’,222),)aSK 

CALL* UERASE  Y  ™EN 

at  sj«(cLRBcK) 


Hi  Ms * cW,L 


OF  POINTS  TO  DESCRIBE  AREA — ( DEF =0 )  ) 
'FAULT'  FILE  )  ’ 

C0N=0 


m 


43,  ,ERR= 


....  JCON 
=  DELETES ) 


8047 


GO  TO  8047 
8046  JCgN=-l 

?iic’ST*TUSs  DELETE  ) 

INUE 

teH]°yrE0-21 60  T°  8044 

u,  F  jL^{ONVA^TATUS=  ’OLD’,  ERR=8044 ) 


1ST  AND  LAST  FAULT  POINTS  ELEV  UITH  ALL  GIVEN  POINTS 

102  1=1.2 


XRxXtHfi 
YR:YTM  1 


THEN 
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EH,  DIF 

IF  1 1.  E-j.il  THEN 

IH=~-ORT 1 1  sTNi  !L-1  l-KTf1(  III!' 

IF  I IH.EG.O!  I H= 1 

XIH=IH,  .  , 


’2  +  1  ,'Tmi  IL-1  '-YTft!  IL ! !  "2 :  /20. 


XR=xth{IL-l  J  +  (IH-1]MXTM{  IL)-Xin(IL-l  j  )/xih 
YR=YTf1( IL-1  +  IH-lJMYTri  IL  -YTH  IL-l)  /XIH 
ENOIF 

Dp|(XR-xljj  j^2+^YR-Y(  J))'*2 

ZR( lj-ifij 
60  TO  1102 
ENOIF  ,  ,, 

ST=ST+Z  Jj/OO 

sb=sb+i./d6 

1103  CONTINUE, 

ZRI IJ-ST/SB 
1102  CONTINUE 

OPEN! 19.FILE=  SKUNK’] 

00  520  |=1 ^NP 

CALL  AREAixtM(iK).XTM(IK+l),X(I),YTH(IK),YTn(IK*l),Y(I),Al) 
IF(Al.GE.O.O)  GO  TO  1521 


-  j. - ...luiinj  inmiinTiM  A\iM  i  ...» _ _  .  .... 

IFTAI.GE.O.Oi  GO  to  1521 

519  CONTINUE . .  , 

WRITE C 19, * )  X(I),Y(I),Z(I),NCC(I) 

GO  TO  520 

1521  WRITE (28, * )  X ( I ) , Y C I ) , Z ( I ) 

520  CONTINUE 
ENDFILE  19 
REWIND  19 
NP-i 

1817  READ(19,',END=1818)  X(NP),Y(NP),Z(NP),NCC(NP} 
NP-NP+1 

GO  TO  1817 

1818  NP=NP-i 
NNPP=NP+ 1 


CALL  USET1C0L21 
DO  1701 ,IK=l,ILrl 


IH=S9RT(  XTM (IK  -X 

S#1*-"1  1h.| 

»8  ™  i;M» 

X  NP]=XTf1  IK  +  -1 
Y  t  NP ) =YTMi IK  j  +  -1 
CALL  UnOVElX(NP)-3 
CALL  UPENlXfNP  +3.,V(NP] 
CALL  UMOVElxfN  1  . YtNP)-3 
CALL  UPENlXtNP  .tlNPhJ. 
SO=T  SORT  ft  XHAX-if1INJ**2+ 
NCCjNP  =0 


XJM|IK|-XTH(IK+1))**2+(YTH(IK)-YTM(IK+1))**2)/20. 


m 


XTfl  K+ll-XTM 
|TH  K+IJ-YTM 


YflAX-YMIN ) * *2 ) /20 .  )**2 


GO  t6  1704 

ST=It+Z(31/D0 
SB=SB+ 1 . 7o6 
1706  CONTINUE 

IF?iAKSHN(IK) 


X?NP|-Xtj|pl+(Y(NP)-Y(j))*'2 
O.EO.O. )  THEN 


F  IAKSHN  IK 
F  IAKSHNUK 
F(ipO.EO.l) 
EUINO  28 
£AD28  *  EN 


.EQ.O  100=1 
.EO. 1. AND. I.EO. 1 
.EQ.  LAND.  IAKSHN 
THEN 


100=1 

K+l) .  EQ.O. AND. I.EO. IH)  100=1 


(NP1-X 

.EQ.O. 


ZINP]=Z8 
GO  TO  1704 
ENOIF  , 
$T=$T+Z8/D0 
S8=SB+ 1 . /DO 
GO  TO  6071 
6072  CONTINUE 


ZfNP  rST/SB 
1704  CONTINUE 
1702  CONTINUE 


CONTINUE 

rii  L  !ic.FTi(n|_nfi  i 
Z  •  NNPP  T  =  7R  !  i  I 
Z'NPI=ZR!2> 

CALL  UFLUSH 
CALL  U6ELL 
CALL  UpAUSE 

CL0SE128. STATUS-' DELETE’ 
GO  TO  1523 
ENDIF 

CALL  CLEAR 
PRINT’,,  1 

t” 


: 1  M|"M  m 


PRIN. 

PRINT’,’  --—THE  ,6 
CALL  UMOVE  f 300. , 300 ! T 

call8ubell’JC0n 

W»*p’Wfl 


T  Of’ 
•N6  CR' 
>  KEYS.’ 


THEN 

THEN 

2+TyP-Y(IC))”2 

THEN 


LC0N1 1 1= _ 
CALL  USET 
CALL  UMOV 
CL  UPEN 
LL  UMOV 
LL  UPEN 
LL  USET 
"  unov! 

LCON 


, X (LCON (J+1)),X(I),Y( LCON ( J ) ) , Y ( LCON ( J+ 1 ) ) , 
10 ,5‘ 

«!  SI  ii  *(I),vm.Z(It.NCCUl 

932  CONTINUE 

523  CONTINUE 
REWIND  19 
NP:1 

252  REAO(19,’,END=254)  X(NP),Y(NP),Z(NP),NCC(NP) 

NP=NP+1 
GO  TO  252 

254  CONTINUE 
NP=NP-1 
1523  CONTINUE 

CL0SE(19. STATUS^  DELETE  ) 

CALL  UERiSE, 

CALL  UBACKG(CLRBCK) 

GO  TO  1077 
ENDIF 
ENDIF 
NC0N=0 


22 


CALL  CLEAR  ,  , 
IFIJASK.NE.  Y 


.  NT’ 
THEN 


ASK. EO.’  Y 
CALL  UE&ASE, 

CALL  UBACKGf CLR8CK 1 
CALL  UOUiLN 
DO  3344  IM,NP 
CALL  UMOV- 
CALL  UPEN 
CALL  UHOVT 

Call  ijpentx 

3344  CONTINUE 


IS  THERE  A  NON-CONTOURED  AREA  (Y/N)-? 
,’ . WAIT’ 


B31 


AAflJ  .* AM  o  F Am 
PRINT’, 


+  OR  ENTER  NAME  OF  A  FAULT 
READ  U40)  ,  CONVAL 
IF(CONVAL.EQ.'  :  NC0N=0 

" . V-wOT 

,A40r  ‘ 


OF  POINTS  TO  DESCRIBE  APEA--ID£F=0)  ) 
FILE  )  ! 


IFjCQNVA! 

uritea 
remind 


R 


CONVAL 


43.  , ERR=< 


h&HP) 


GO  TO  6605 
6603  NCON'-l 

CL0S£(43,$TATUS=’DELETE’ 

6605  CONTINUE 

IFJNCON.EQ. 1 . OR . NCON. EO. 2 )  GO  TO  6601 
IFINCON.LT. Ol  THEN  ,  , 

0PENU4, FILE=CONVAL. $TATUS=  OLD  .ERR-6601 
NFALT=1 


IL=1 


71  CLOSE ( 44 ) 

Tl  zTI  -1 

[19,FILE=’ SKUNK ’ 


NNPP=NP+ 

DO  1404 
IH=SORT , ... 

IF  IH.Ed.O 
XIH=IH  , 

S8lMiEIW 

NP=NP+1 


Xfnllk)-XTH(IK+1)  )**2+(YTM(IK)-YTf1(IK4l)  )**2)/20. 


XlNP}:XTH(lK)+jl-l}MXTH|]K+l)-XTH(IK))/XIH 

NcSpi-o K  +  _i  YTn  K+i  *yth  K  /xih 


1405  COI 
1404  COnV 


AKSHN 
AKSHN 
INUE  , 

USET  COLOR) 
...  INUE 
CALL  UFLUSH 
CALL  UBELL 
CALL  UPAUSE 
0  52|  1=1 jNNPP-1 


IK+1J .EQ.O. AND. I.EQ. IH)HARK (NP ) -1 


DO  522 


522  CONTINUE  , 
URITE1 19, * 
521  CONTINUE 


31] 


XTM(IK+1 ) ,  X(I ) ,  YTM(IK) , YTHl IK+ 1 ) , Y ( I ) , A1 ) 
GO  To  521 


SO 


= I SORTt (XWAX-XMIN ) * 
I:NNPP-NP 

ST=0 . 0 

DO  1408  J= 1 , NNP P — 1 
A2=0 . 0 

.Up  AL^  POINTS  TO  COUP 


Tia 


ID0=O 


.NNPP.OR.I.EQ.NP 


*2+ ( YMAX- YMIN) *  *  2 ) / 20 . 1**2 


EV  OF  1ST  AND  LAST  FAULT  POINTS.. 
TO  1411 


N(IK).EQ.O)  ID0=1 
IlfHNPP^P- 


IF 1 1 AKSHN 
DO  6302  I. 

IF  ( MARK (I  I ) . EO 
6302  CONTINUE 
IF  1  IDO. EO. 1 
DO  1409  IK= 

CALL  DL ARE A  , 
IF(Al.GT.O.O)  GO  T 


ID0=1 


GO  TO  1411 

IL,1 


XTM(IK+l).X(J).YTM(IK),YTM(IK+l),Y(J).All 

1411 
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lino  row T t M if 
■-O' I?  140: 

Hit  Ci'iNTTNIlf 

'  06*lxUI-xiJ)]”2*f>fH-Y(j))*-; 

IF ( DD. tO. 0, )  THEN 

mum) 

Ed  tO  1407 
ENDIF  ,  ,, 

ST-ST+Z(J]/00 
SB-SB+l./dO 
HOB  CONTINUE 

z{i]=st7sb 

1407  CONTINUE 
GO  TO  958 
ENDIF 

NUftCON=NCON 
.  CALL  CLfAR 
PR  NT*,,  (2.11) 

PR  NT*.  DESCRIBE  AREA  BY  POSITI 
PR  NT*.  POINTS  (IN  A  COUNTERCLOi 
PR  NT*,  CURSOR  IS  flOVED  WITH  Al 
PR  NT*’  ---THE  6.U  KEYS  G 
CALL  Ufl0VE(300.,300.) 

DO  2086  I=1,NC6N 

CALL  KUR^N(XP,YP,IP) 
DCL=999.E30 

?P(X?fcT?GT^pf  THEN 
IF(y(ICI.GT.YPI  THEN  ,  .... 

{(MM  ®"-,IICI1  2 

Mi 

1 

2087  CONTINUE 


AREA  BY  POSITIONING  CURSOR  AT  LOWER  LEFT  Of’  , 
N  A  COUNTERCLOCKWISE  HANNER)  AND  PRESSING  C?  . 

5  m  m  IrayUt.'r ' K£,s- 


il 

CALL  UttOVE 
CALL  UPEN 
CALL  USETl 
CALL  UflOVE 


DO  889?  Hl.NUHCON 
8899  KCONilHcONU)  ,  , 

LCON  NtON+lHLCONl) 

OPEN! 18. FILE-  SKUNi 
REWIND  18 
DO  952  1=1, NP 

call5dlaremxIlcon(j)),x(lcon(j+i)),x(i),y(lcon(j)),y(lcon(j+i)), 

IFfA.LT.-0.60i  00  TO  952 
954  CONTINUE 
WRITE! 18,  )  I 
952  CONTINUE 


REWIND  18 
0PENI19.FILE=  SKUNK  ) 

REWIND  i? 

976  READI18  ‘  END=958)  I 

call  dlarea|x|lcon(J)),x(lcon(j+i)),x(i),y(lcon(j)),y(lcon!j+i)), 
+  IF(ABS(A)At!oAi)  THEN 


REWIND  1? 
READ(18.‘.E 
DO  957  M. 
CALL  DLAREA 
t  Y 

IF ( ABS ( A ) . L 
WRITE (19 . *  T 
GO  TO  m 


ENDIF 

957  CONTINUE 
GO  TO  976 

958  CONTINUE 
REWINO  19 
NCON=  1 

1952  READU9  ’.ENDM954!  I 

lcon(nc6n)=i 

mw 

1954  NC0N=NC0N-1 
REWIND  19 
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CLOSE* 1; ! 

END  If 

OFENI li . FILE=  TAP. OAT  .STATUS:  UNKNOWN  » 
REWIND  11 


1$ITe!i  :1)NX!I), 
XPL:XUH  N 


Y ! I ) . 2 1 1 ! , NCC ! 1 ) . ASPACE , ALABEL ( I ) 


XPR:XWMAX 
YPL = YUtl  IN 
YPR:YUf1AX 

IF (JNPLT.EQ. 1  THEN 
CALL  CLEAR 

fir  v  uant  10  (cilip  oR  ,z)oon  area  (def=n)  ?(c  or  zi  >5 

c  C°?  eIr' ' ^  R- Z00M- E0-  ’ 2 ’ 1  THEN 

■L  ■  1  “ .’  UINOOU  FROM  (F)ILE  OR  (C)URSOR  (DEF=C)  >  ' 


IF ( ZOOM  Efi. ' C’ . OR. ZOOM. EO. ’ 2 ’ )  THEN 

pfilNT„TA\)Y  UINOOW  FROM  (F)ILE  OR  (C)URSOR  (DEF=C)  >  ' 

ffflo  AIWgO  TO  8346 
CALL  CLEAR 

PRINT  ,, ENTER  NAME  OF  FILE  WHICH  CONTAINS  COORDINATES  OF’ 
PRINT,  AY;,’  LOWER  LEFT  CORNER  AND  WIDTH  OF  WINDOW  )  ' 
READ,  (A2Q)  FORCFILE  ,  , 

Anrn  (it  ni  rifAnrcti  r  triTiir.1  ai  r\’  roo-rti/n 


OPENI43. FILE^FORCFILE. STATUS:’ OLD’ ,ERR=9346 
GO  Tu  9347 

9346  CALL  CLEAR,,  , 

PRINT,  (AX  ’  ’ ’  FILE  ’, FORCFILE,’  NOT  FOUND. .PRESS  RETURN  >  ’ 
READ  ’  Al]\pQT 
CALL  CLEAR 
GO  TO  9367 

9347  CONTINUE  . 


READU3  .*) 

mil 

YPL= {  YPL- 


XPL , YPL , WIDT 


»-$ubx-xhinx  /sn Moo. 

-SUBY-YMINY  /$H  +100. 

/sn 

XPR=XPL+UIDT 
YPR=YPL 
GO  TO  8345 
CONTINUE 

PRINT'  VUSE  CURSOR  TO  LOCATE  LOWER  LEFT  CORNER  OF  AREA’ 

i  t  ft  ft:  ( P| 

CALL  UMOVE{XPL-io.,*PL) 

CALL  UPENfXPL+10, ,iPL) 

CALL  UHOVEIXPL.yRL-IO, 

CALL  UPENTxPL, VPL+10. ) 

PR^NT^^USE  CURSOR, TO  INDICATE  WIDTH  (OR  HEIGHT)  OF  AREA’ 
CALL  (jflQVEuPL.YPLl  ... 


CA 

8345  CO 


XPR, YPR, IP) 


ABS  XPR-XPL  .GT. 
(ABS  YPR-YPLI.GT. 


XUnltUXPL 

XWHAX=XPR 

YWflIN:YPL 

YWHAX=YPR 

ENDIF 

CALL  UERASE 
CALL  CLEAR 
CALL  PEND 
END  IF 

2301  CONTINUE 
ZflIN=9.E30 
ZMAX:-MN 

M.?!i 


:!81l 


1=1 ,NP 
.GT.XWHIN.AI 


+  Y  .LT.YWMAX  THEN  ,  , 

|  J1  :MI  »H) 


.  LT.XWMAX.AND.YlD.GT.YUMIN.  AND. 


2307  CON 
PR  I 


(2. 13)’ 


PR  NT’.  MIN  EL. =  ’.ZMIN  ’  MAX  EL.=  . ZMAX 

PR  NT'.  CNTER  BEGINNING  CONTOUR.  CONTOUR  INTERVAL,  HEAVY, LINE’ 

PR  NT',’  MULTIPLE.  AND  NO.  OF  CONTOURS! 0  MEANS  SPAfl  RANGE). 

+^RORT;’r  F flR  NO  CONTOURING  (DEF=AUTO)  (H  FOR  HELP)  >  ’ 

I^CONVALEO'.^NONEIt.OR.CONVAL.EO.’N’)  CONVAL:’ 999999,1.0,1’ 

rr  i  A/Miuii  rr\  ’  u  ’  »  rucki 


IrINT'^FOUR’ NUHBERSEARE  TO  BE  ENTERED ( OR  ENTER  RETURN  FOR  AUTO)' 
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♦'THF  j$T  IS  THE  BEGINNING  CONTOUR  VALUE.  THE  2ND  Is  A  *  0* 
Pr{nT\[A^E'3RDCCAU5ES  ALI^CONTOUR  LINES  WHICH  ARE  MULTIPLES,’ 


yraffii'0  MWffi  « 
sW-feW  ?«K  l\  ® 


(PLOTTER]  OR  ANO„.. 
NO.  OF  LINES  TO  6E 

Ml  i' 


»*  ITT 5020.-5  Io:6  HEAM5  START  ST  f[EV  5020  AND  ORAUJHE’ 
*;  5656, 5018, 50i6,50U. 5012.  AND  5010  CONTOURS . UITH  THE 
PRINT*.  5020, AND  5010  tflULTtPLES  OF  10)  DRA^N  BOLD. 

W’  A F  ^T"pfi^ss  *ETURN  T0  C0NI1NUE 

I Mti/U.™" 

8L0G=LO6l0(ZD j 

ImMU®:1- 

IF[ ^OEL. GT . 5 . )  JDEL=10 


L0G=L0G-1. 
LOG) I+l. 

JOEL-10 

0EL=4. 


ll=!zO?ZINt  +0.0001*ZI1IN/ZINT 

ZZHIN=IZ*ZINT+ZINT 

NCH=Q 

GO  TO  7503 

ENDIF 

OPEN (43, FILE=  PI345ER  ) 

URITEUJ.MAAOr)  CONVAL 

B43^,ERR=2302)  ZZtIIN.ZINT, HEAVY, NCfl 
AVf.fO.Q J  60  TO  230$ 

^^AEA^Y  lI2eTHuC!tIPLE  MJST  BE  GREA 
"  TO  5302 


MJST  BE  GREATER  )  OR  =  TO  INTERVAL’ 


GO  TO  2303 

ot(imtatuss  DELETn 

CL0SE?43, STATUS=*  DELETE ’ ) 

SWBB.EQ.1J  NED=0 

IF  NFA(.T.EO,lJ  GO  TO  7564 

’2,14  LEAVE  BOUNDARY  OF  CONTOURED  REGION  UNDRAWN  ?  (Y/N)  )  ’ 

BFAM*  SSoT  TACK 


»• 

IFIJASK 


22)  JASK 


7564  PRINT,’  A\  ,’,’  (2. 
READ  ’(A20) ’  .TfIPLT 

Jffik;' 


ENTER  TEMPLATE  FILE  NAME  (DEF=NONE )  ) 


ITHPL=1 


5it^-o’2  2  JASK 


TO  SEE  ELEVATIONS  PLOTTED  ?  (Y/N)  ) 


IFTJASK.EO.’Y’  JEL 
jlabel=6 

IFfJELV.EO.O  THEN 
JTtST =0 

DO  3200  1=1. NP  , 
IF(ALABEL(T).NE. 
3200  CONTINUE 

T 


JTEST= 1 


’  (2, 17b]  WANT  LABELS  PLOTTED?  (Y/N)  ) 

^ffiKlM^JLABEL=l 


17c)  LEAVE  X-Y  AXES  ANNOTATIONS  OFF  ?  <Y/N)  )  ’ 


♦S (2.18)  LEAVE, LOCATION, TIC’’S  AT  DATA  PJS.  OFF  TOO  ?  (Y/N)  >  ; 
PRINT  ' (A\  J ' .  (ENTER  TO  INCLUDE  GRID  POINTS  ALSO)  )  ' 

wlw22]  PASK 


IFIPASK.EO. ’ Y’ ) 
IF(JELV.EQ.l)  J 


JJELV=1 

JELV=0 
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PRINT*.  !'  i0.1  ENTEf?  TITLE  Of  FLO? — UP  NONE.  HIT  RETURN)' 
PRINT  1  AT 

READ! * ,  TA48\ ’ T  TITLE 

PRINT’.’  -PRESS  RETURN  AT  BEEP?' 

CALL  U8ELL 
CALL  UPAUSE 
CALL  PENO 

^OOhIeq!  C*  CALL  UUINDOl  1 . ,  1040.  ,0. ,  780. 1 
CALL  U^T!C0L2) 

inZoHn.^’Z’)  CALL  UUINOO ( XPL , XPR , VPL , YPR ) 


i|?z5oA:eo. ’z! )  FSOf 

IttMl 

DO  1929  1=1, SiP 
IF  PASK.NE .V  THEN 
IF  NCCfij.EO.O!  GO  TO  1929 
IFINCClll.LT. 1000)  GO  TO  1929 

CALL  UnOVE(X(I)-3.*FSOFT,y(I  ) 

CALL  UPENim  +3  ‘FSOFT.Hi!  , 

CALL  UHOVEfXf  ) , YUJ-3  ‘FSOF 

CALL  UPENlX  I  ,Y(1 )+3. *FSOFT 

CONTINUE 

CALL  UHOHE 

CALL  UBELL 

CALL  UPAUSE 

ENDIF 

CALL  USETf  SOFT  ) 

CALL  UPSET f VERT' , 30. ’PSFT’FSOFT 

CALL  UPSET!  HORI ’ , 15. ’PSFT’FSOFT 1 

IFIAASK.NE.  Y  )  CALL  XYAXISt XNINX , XOAXX , YhlNY , YNAXY, Stt, 0) 

If  Ncchf  EO*o]  GOETO  1923 
IF ! NCC IIJ.LT. lOOOJ  GO  TO  1923 

CAL1UN0VE ( X ( I ) *6, ’FSOFT , Y ( I ) +6. ’FSOFT ) 

CALL  UBELL 
CALL  UPAUSE 

IFIJLABEL.EQ.1)  THEN 

H  - 1  |^p 

I  7e6. 0)  60  TO  1223 

!M;iF  60 10  l2?l go  to  mo 


FSOFT=( { XPR-XPL ) / 1039 . ) 


THEN 

GO  TO  1929 
0)  GO  TO  1929 


If  BttBiHjIp1 ;° 1 

Myi(||IT*6.,Y(I)*6. 

^  BEL(IJ,’(10A1)')  ( 

hr!  jPfi  a|el(  i ) ,  ’  llOAl?’ ) 

GO  TO  1229 


GO  TO  1223 


(BLAB(J), J=l,10) 
(8LAB(JJ),JJ=J,10) 


nlK’jlOAl)’)  ( BLAB (J) ,J=1, 10) 
O.’’  »  60  TO  1226 


UR^TeVbCABEL. ’ (11A1 ) ’ )  (8LA8(J),J 

mfo;l?:,i8/PSFrFj»Fi) 

ULCuPSET  ’HORI’  19. *PSFT*FSOFT) 
CALL  UPRNTlfBLABEL, ’TEXT’ ) 

2) .THEN, 

CALL  UPSET  VERT  ,30. 'PSFT’FSOFT 
CALL  UPSET! ’HORI’. 15. ’PSFT’FSOFT 


END  IF 

CALL  UBELL 
CALL  UPAUSE 
ENDIF 

CALL  USET (COLOR I 
IFflTriPL. EOf It  THEN 
CALL  UCHPENlLPEN(3) 
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'.STATUS1' OLC1 


"  F  S  =  1 1 


59 


159 


12.  F !LE=TMPL 
]PN:-1 

REApii2.\END=5<>i  X7.Y7.IFN 

!?:  i?:»  W 
is  I??:8I1J1  Hi.. 

CALL  UUHERE.X77.Y77. 

ANG=ATAN2u7-Y77, X7-X77 
IFUPN.LT. o  calC.  UtlOV 
ifCjpn.ge.oI  call  upen 

CALL  UH0VE[j77+l. ’FS0F 

IFlJPN  .  “ 

FtJP(U 

JfV^PN. LT . Q] CALL  UHOVE ( X7-1 '*FSOFT*SIi 
IFlJPN. GE.Q  CALL  UPENCX7-1 . *FSOFT*SIN 
CALL  UM0VECX7. Y7 
jPN=IPN 
GO  TO  49 
CONTINUE 
CALL  UBELL 
CA^L  UPA 


.  Y77-1 ,’FSQFT’C 


*  COS ( ANG  J ) 
COSIANGJ) 


’CO 

cos 


END  IF 

IF  NCd! GT ; 52 ; and! HEAVY! NO!}  NCH=52*HEAVY/ZINT 
F  NCfl.EO.O.ANO.HEAVY.EO.6.  NCH=52. 


52. AND. HEAVY. EQ.Q.  NCH=52. 


777 


778 


if f nch ; eo ; o ; and ; heavy ; ne ; o : ) 

NCr2 
CL 


NCH=52*HEAVY/ZINT 


NC)-CL(1)+2INT  (NC-1) 
"NT.GT.Q.Q.ANP.CL  NC  .GT.ZHAX  GO  TO 
CLINCJ.LT.ZHIN)  GO  TO 
TO  778 


?78 


T.ZfllN 

.ZHAX 


NC- 


..  _ LT.O.l’HEAVY) 

KF=26 

IF(JJ.GT.26] ,KF=-26 
HLETliJ  =LET  JJtKF) 
NyHHlIjrJJ 


THEN 


EndIf 

7608  CONTINUE 
JO  7607,1=1 


7607 


IMHlIl! 

NTINUE 


.NC  , 

Ae.o) 


GO  TO  1159 


1159 


HEAVY=( 

ENDIF 

,»Nl 

IFlNCON. 

XP»IN=99 _ 

Xf1AX=-99.E30 

YHIN=99,E3Q 

YHAXr-9§.E30 


CIRCULAR  CELLS.. ABOUT  NPCELL 
THEN 


POINTS  PER  CELL. 


DO 

6 

3  1  = 

a 

ON 

h 

X 

ICON 

i 

.LT.XHIN 

XHIN=X  LCON 

I 

F 

X 

ICON 

I 

.GT.XHAX 

XHAX=X  LCON 

I 

F 

Y 

ICON 

i 

.LT.YHIN 

YPIIN=Y  LCON 

I 

F 

Y 

LCON 

1 

.GT.YHAX 

yhax=y!lcon 

I 

AVX  = 

AVY= , 
JID=0 
DO  661  1= 


XfllN+XHAX 

YHIN+YHAX 


2+(Y(LC0N(I))-AVY]”2 


D.G 
CONTINUE 
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LID-0 

DO  >62  1=1. NF 

D=i  X' I)-Ayxr’2+IY(  I)-AVY)**2 
ID=SQRTIDI +1 .0001 
IFII0.Lg.3ID)  L I D=L 1 0+ 1 

mr 


XHIN=99.E3Q 
XHAX=-99. E 30 


Bi=I,np 

■\m 

.LT.YHIN  YHIN=Y  I 
j.GT.YHAX  Yf1AX=Y  I 
Of 

IFINCON.GT.O  GO  TO  697 
AVX=|XHIN+XHAX|/2. 

AVY=f YMIN+YMAX J /2. 

697  CONTINUE 

00  650  1=1.1300 
650  NCCELL  I  =0 
00  660  1=1. NP 

D=(X(I]-AY5(]**2+(Y(I)-AVY)*,2 
ID=S0Rf(0  +1.0001,  , 

660  NCCELLl I0)=NCCELL(I0)+1 

y.1,0 

H1|s° 

00  670  1=1.1300 

no=nd+nccellT i7 

IF(K.E0.2T  NPCELL=LIO 
IF CK.NE.2j  NPCELL=IPCELL 
IF|n£>  GE.NPCELL)  THEN 

IRN  K  =ND 

K=K+1 

ND=0 

670  foNT^IUE 

l!jl?i M  ™ 6M 

k=k+i* 

680  CONTINUE 
KRO-K-1 

CALL  UCHPEN ( LPEN ( 1 ) ) 

£r^nt^a*  wait...’ 

IF(Hfl.EQ.2  .CALL  USET(’HARO’) 

ft  MA  U 
iW  fck)  66  :  fe 


mi 

.LT.YMIN  YMIN=Y 
.GT.YHAX  YHAX=Y 


20N.GT.G  GO  T 
tXflIN+XHAXl/2. 
YHIN+YHAX  /2. 


TO  697 


GO  TO  680 


YUOHIN*  YPL'(65./780 
YUOMAX={YPR* (65. 7780 


CALL  UDAREA ( XUOHIN.  XUDHAX, YUOHIN,  YUOHAX ) 

^pOOUTLN 
CALL  CONTER  ZINT) 

Efttt  ibh  f  m 

SAVE  A  MESH  FtLE  IN  FILE  ’RES&.OAr 


ONTER  ZINT) 

UINQO IXUMIN.XUMAX, YUMIN.YUMAX)  .  , 

DARtAf 10. *Pdx , 75 . ?PCX . lfi.  fCY,75. *PCY) 

C£U  ctl  C  Y tf  CT1  C  ’Mtrcti  fkAT1 


AVE  A  HESH  FILE  IN  FILE  HESH. OAT 
GO  TO  8880 
CONTINUE 

£MlT\^2,21)  00  YOU  WANT  TO  SAVE  THIS  HESH, (SO.THAJ  YOU  HAY’ 
PRINT,'  A\  V  LATER  DRAW  OTHER  PROFILES)  ?  (Y/N)  ) 

mp  'Hv« 

$INT^?2.22)  ENTER  A  FILE  NAflE  IN  WHICH  TO  STORE  THE  HESH.’ 
PRINT,  A\  1QEF=H.HSHJ  ) 

READ  ’(Alil’.F 

ifTfihesh.eo.. 


BtfiMf 

IE  t F  ihesh . eo. 


=H . HSHj 


CALL  CLEA 
PRINT*. \ 

MUi 

URITE (81 , 


FIHESH=;n.nSH’ 
SAVHSH=  H. : : : 


f=?IHESH] 
f  SUBX.SUBY 
I  FACX.FACY 
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WRITES!. 'J 
"RilEf  5l . r 


X01.XQ2.Y01.Y02 
ZflAX .  flllN.  XMINX,  XMAXX 
YMINY.YHAXY.Sf! 
NP.MB.Mei 


URI...... 

WRITE  81.' 

WRITE i 81.’ 

W^ITutl !'l ' XU ) ,  Y(  I ) ,Z! I) , NCC ( I ) 

3356  CONTINUE 

WRI?|t81 ! *  ^ ' ^VERTEX ( I , J ) , J= 1 , 3 ) 

3367  CONTINUE  , 

OPEN  65,FILE=’ED6EC0N’ . FORM= ’ FORMATTED’ 
REA0165  \EN0=1868)  NflU 

-  jf  »T  uhi 


1867  READ! 
URI  i 


1868 


y^i  tLiYVivvui  n 

*  «a«i » 


X1.Y1.X2.Y2 

Yi 


CLOSE  ILsIaTUS*’ DELETE’ 

60  TO  18^ 

ENDIF 

1823  REWIND  65 

IKN=0  ,,  ,, 

1827  READJ65,  (A72)  ,END=1828  DUM 
IKNUKNil 

60  TO  1827, 

1828  WRITE(81j  )  IKN 
REWIND  65 

1842  READ (65 , * , END=1843 )  X1.Y1.X2.Y2 
WRITE181 Xl.Yl.X2.Yi 

60  TO, 18*2 

1843  CLOSE (65, STATUS=  DELETE  ) 

1824  SPM8^FILE=’SAVC0N’) 

1767  If aJ|$5j  , END=1768 ) HN1 , MN2, HN3, HN4 , flN5,HN6 
+&R  jtE  j&lM2f1Nl ,  HN2,ftN3,MN4 ,  MN5, HN6 

1W™ 


ALL  CL 
60°fo  8882 

C  ...........  DRAWING  TITLE  AND  LE6END 

8880  CONtlNUE 

if  inn.  eo.  l .  or.  tin.  eo.  2]  then 

CALL  UWINDOll. , 1040. .0., 780. 

CALL  UDAREA ( 10. *PCX. 100.  Drv 

Stt  USpHflHjPSFT) 

SALL  UPSET t’HORI’ . 15. ‘PSFTl, 
FIHH.E0.2I  CALL  USETt^HARD’ 

ND  if 

fTiim.eo 

|[  jjJftlWitl 

CALL 
CALI 
R 
T 


PCX, 0. *PCY, 100. *PCY) 


HORI’ ,15‘PSFT 
CALL  USETTHAR 

.THEN 

jftW0.,0.,67.] 

itroT'  'in  > 


') 


I  ( TITTLE C I ) , 1=1,48) 
-1,48) 


>|lf  LE6END  ’ 
' , JASK 


_  _  (Y/N) 

READ  UAir, . 

f^ASK^Q.’Y’)  60  TO  162 
IPOW=NC 
KPOW-O 

DO, 8900, 1=1. NC  , 

IF  NUHH(I J.fiT.O)  KP0W=KP0U+1 
8900  CONTINUE 

IF(KP0W.6T.0I  I POW=KPOW 
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CALL  INF < IPOU • 

PJr775 

YP= 7nn 

IFfiP0U.LT.26) 
IF  ( YF .  LT . 500 . ! 
CALL  UflOVET 
CALL  UPRNTl 
CALL  UPRNT1 
I NC = 1 


YF  =•  IPOU  *  25 .  +50 . 
F=50‘ 


"ill* 

ZINT.’I 


f?; 

EAL 


XT’] 


00,500 
IF(NUHHilC 


0  IC=  .N 
flH  IC  .E 
. LT. 5  .) 


NC  INC 

’■fl!0. AND. HEAVY. NE.O.  ] 

th£n 


IF  ( YF 

m 

CALLCU^(PX,YF) 

KF;-26 

IF jIC.LE.26)  KF=26 
" HEAVY. EO.O.  CALL  UPRNT1 


GO  TO  500 


IF  . . 

IF (HEAVY.NE-O. 
CALL  UPRNTl 


l  vi  nn  in  \  «  TEXT 
ACL, GT. -0.0001. AND 


,C^LL  UfRNTKHLE 


ISHIklfSIxri 


500 

162 

8882 


IF  ABS  ACL 
IF  ABSfACL. 
IF ( IC. EO. 1, 
CALL  UPRNTl 
CALL  UUHEREI 
END  IF 
YF=YF-25. 
CONTINUE 
CALL  UBELL 


CALL 
CALL 


.LT. 0.0001)  ACL-Q, 

uprnti[acl,;real 

UPRNTl (ACL, ’ INTE’ ) 


C. 


UPAUSE 
[NUE81 


If 

ft 


CAD 

m 


l>|n6 


PROFILES. 

PW' 


CLEAR 


. PROFILES. 

CALL  USET(COLOR) 

CALL  UFLU$H 
956  CONTINUE 
REUIND  09 
XBI=MB 


GENERATE  A  TRIPLE  TRIANGLE  HESH  ?  (Y/N)  >’ 


THEN 


AR  OR 
ASK 


INVERSE  SQUARE  RULE  ?  (L  OR  I)  DEF=L  > 


00  516  I =MBL , MB 
II=I 


NPUT=0 

UT  DATA  FILE  NAME  IS 
WAIT.... 


'TRIPLE' 


,CX(2],CX(3),CY(1),CY(2),CY(3),A) 
516 


vl:lcZjipt!tCZ?2  WS2+ci?|i/si)  /U. /Sl*l.  /S2»l . 


/S3) 


URITEI 13.  ’  !  "'IIRY.SIIBY 
•51?  READ' i  1 . 1 . ENO=aiSl  5y.By.BZ 
gx=s(t1[Tg>-inn.  i+xniNx 
S Y=SM  f 5Y- 1 00 . 1  +  Yfl  I N Y 
! 13,1!  B  _ 


IZZ.ASFACE.ALABELf 1 1 


BX.BY.BZ, IZZ.A5PACE. ALABEL ( 1 ) 


RE,. 

222  FORM 

CALL  _  , 

IFUASK.EQ.’Y’)  THEN 
CALL  UMOVE  20..560.) 


B  ’  TRIANGLES  FORMED  ’ 

‘  YOU  WANT  TO  DRAW  THE  MESH  OF  TRIANGLES-?  (Y/N)  >’ 
JASK 


MM.EO.l 
MM.EO. ‘ 
MM.EO 
"M. 


Y’  THEN 

O. ,560. I  ,  . 

CALL  UDAREAtlO. *PCX,75. ’PCX, 10. *PCY,75.’PCY) 
CALL  UDAREA  lQ.’PCX  75,’PCX  10.*PCY  75.‘PCY 


E0.1.0R.MM.E0.2)  CALL  USEK’HARD’l 


510 

513 


RAW  T 
F(MM. 

CLR=1 
.0510  I=MBL , MB 

X8I=II 
CALL  CORCO 
CALL  UNOVr 
CALL  UPEN 
CALL  UPEN,...,. 
TALL  UPEN  CX  1 
.R=ICLR+1  , 
1CLR.GT.7) 

mp 

AUSE 


a®1 

RETURN 
END 


YflAX.SM,  IFLA6) 


K.PSIZE.LPEN18) 


o.cw.U  men  , 

XUDMIN-10. ‘PCX^ . GT . 0. 001 . AND. ABSf  XUOMAX-75. ‘PCX ) 
^A^UoIrW  k  ^PCX;  *P^x;  k ‘PCY.75. ‘PCYl 


If.UOUTLN 

Mm 

zn=fXMAX- 
"  ‘S=L0G10 


AX-XMIN) /SMl+100. 
125, 


[IQEL.GT.5 

'  M 


m)lm 

JOEL-10 


‘COG) 

)  +0. 0001 *XHIN/XINT 


20 


XX=ZZMIN-XINT 
XX=XX+XINT  , 
X=M  XX-XMIN) /« 
IFIX.LT. GO. ) 


^  SM )  +  100 . 
GO  TO  20 
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TFiK.fiE  nnAX+xci  r,n  m  ;.n 
r A!  i  Hfi'ivp  i  x . n  1 
CALL  UPENIX,  10. 1 
CALL  UMOVA | X. 34.  t  ,, 

USETT  ACENTER  1 


IfVxINT7FACX.LT. 1.  CALL  UPRNTll  XX  FACX  +SUBX, , REAL . 
IFfXINT/FACX.GE.  1 . 1, CALL  UPRNT1 ! (XX/FACXHSUBX.  INTE 
CALL  USET?  NOCENTER  ) 


30 


40 


CALL  UMOVEfX, 770, 

CALL  UPEN(X, >80. 

GO  TO  20 

CONTINUE  1(  , 

YYHAX=I[YI1AX-YNIN)/St1)+100. 


NUM 

ZD- 

BLO 


NTINUE 
MAX=([ YMAX- 
M-YYHAX/125 
-(YMAX-YOIN 


LOG^BLOG 
IF (BLOG. 
IDELMZD 
JOEL=IDE 


LOGlOfZO 


/NUM 


[ho0:?-! 

IfT ioelTgt . 5 )  ■ 

IF  1 1  DEL . EQ , 3) 
YINT=JDEL* ( 10. 
YS=YINT/SM  , 


)  LOG=LOG-1 
L0G))+1 


JDEL=10 

+0.0001*YMIN/YINT 


YY=ZZMIN-YINT 
YY7YY+YINT  ,,  , 

Y=  YY-YMIN  /SI1  +100. 
Y.LT.50. |  GO  TO  40 
“  . AX+YJ)  GO  TO  50 


ALL  UMOVE 


Y.GE.YYM 


CALL  UPENTlO. ,Y) 
IF( YINT7FACY. feE. 1 . 
- ACY.LT.l. 

urn] 


50 


IF  YINT. 

m  w . 

GO  TO  40 
CONTINUE 


CALL  UPRNT1 
CALL  UPRNT1 


({YY/FACY  +SUBY , ’ INTE’ 
( (YY/FACY  +SUBY ,  REAL 


:°. 

)0  [  AHttin «  AvnmA  «  iwiuw«  iwiihaj 

:Al  XUDHIN.  XUDHAx,  YUOhIN,  YUOflAX) 


IF? _ 

CALL  UUINI 
CALL  UDARI 
ENDIF 
RETURN 
END 


SUBROUTINE  DRUPRF(MM) 
COHMON/DVIC  " 


1)  THEN 

WHIN . XWtIAX,  YWffIN.  YIIMAX) 


YLAB2’41 


^P4^LX48KY*^i^LOR?CLRBCK’PSI2E’LPENf8) 

„SC,70.SECM1$Lll.SEC2*2.DE$C2*71.YLABEL,40,Y 
DPR^72, JASK  l.REFEL’AO.OR^FEL’AO.N^ULAB'AO 
YYO  2a|SHADE* 1 . LEVOFF* i , ZEFEL  40, 0ZEFEL*40 
OSTADE*!.  OEVOFF’l,  JK  ASK  *  1 ,  SYM*  1 , 1.EG*20 . GEL’21 . LR’ 1 


CHARACTER 
CHARACTE’ 

CHARACTE 
CHARACTE 
CHARACTE 

rHAPAT  tP 

OIMEN'1^^260] ,  Y(260) ,  Z(260) ,  IB0X(80,80) 


L  UCHPENlLPEN(ll) 
R.IIM.EO. 

ACT’  ] 

J 


CAL _ 

IFTMM.EO. 

CALL  USET 
CALL  UPSE  , , 
CALL  UPSET  ’HOR 
ENDIF 

IFiMM.EO. 2) 

CALL  USET t ’ 
ENDIF 


I  VI'  .  IIL 

Soft', 
'  VERj: 

THE(I 

SOFT') 


2)  THEN 

,:psft 

,'PSFT 


m 

REAC 

REaE 


REA.. 

WRITE 

WRIT- 

WRIT 


VX1.VY1 

Air ' 

A4C 


YL  AB> 


ZflAXr-9.  E09 
Zt1IN=-ZMAX 
YHIN=ZMIN 
XMIN=ZMIN 
YM  A  X  =  ZflA  X 
XMAX=ZMAX 

141  READ! 14. ’.END-143) 

xTu=xfi ) -100. 


DESC.SL 
,  SL 

YLABEL.SL 


X!l!,Y!l),z;i) 


BA  2 


v :  t i = v i t i - 1  no 
iFTxilJ.'LT.'xniMJ  XHIN-'XII  j 

IE !  I  I 

1 1 mwm 
IF? I2Z.EO.  -9*59)  THEN 
1  =  1  +  1 
GO  TO  HI 

Kill 

GO  ™  141 
143  CONTINUE 

closeIh) 

XSUB=X  1 1 1 
YSUB=V(1 j 
CALL  Cl  EAR  , 

w-®  v  profile  °f  a  series  °f  plots  ?  y/n  ) 

IF( JKA$K.EQ.  Y  j  THEN 
DO  101  1=1,80 
DO  101  J=1 . 80 
101  IBOXu,  J)=0 

IF(JKASK.EQ.’Y’)  THEN 
CALL  CLEAR 

PRINT*, , Verbal  axis  description:  .YLABEL 
PRINT,  A\]  /  Enter  new  description  OEF=NO  CHANGE  >  ’ 
READ  (A40)  jNEULAB 
IFINEWLA8.NE. ’  W  YLA8EL=NEULA8 


READ  (A4d)  ,NE1 
IFINEWLA8.NE. ’  1 
CALL  CLEAR 
ENDIF 

ISYt1=0  , 

IF  JKASK.EO.’Y’] 

In  JKaIk^NE.  ’Y’  ] 
HSYM=HSYH+4 
ISYH=MSYM 
ENDIF 

IFUlSYfl.GE.  16)  t 

if&askJeq/y’! 


MSYH=0 


16)  M$YH=0 


?R  NT  (A\]', 

+  12.37)  UANT  THE  PROFILE  STRETCHED  TO  FILL  PAGE  ?  Y/N  >  ’ 


CALL  CLEAR 
DO  200  1=1, NJ  , 
X(I!=(X(I  -XSUB 1*0.639 
200  Y(i!=(Y(I  -YSUB  *0.639 


XF=1 . 0  ,  , 

IFfJASK.EO.  Y 


XF=664 . 56/SORT ( (X(NJ)-X( 1) )**2+(Y(NJ)-Y(l) ) **2) 


201  Y  I  =Y( 
CALL  CL 

PPTWT*  i 


IfilwiF 

777  c8nTIN0EH  ttIN  EL=  'F8-2,11H  m  ELz  ,F8-2) 

IP'5' 


if  »i 

PR  NT  V  A 

. 7  A  701  I 


AND  znAX  ,oef:auto) 


XSET=200. 
YSET= 100. 
YGAP= 1 25 . 

CALL  CLEAR  , 

mm., 

IFIREFEL.EO.' 
ZD=  (ZflAX-Zflif 
BLGG=L0G10 ( Z( 
LQG=8L0G 

lolE=Mi8: 

JDEL--IDEL 


REFEL=0REFEL 


B4  3 


JF  !  I f)£i  Qr  f.  1  JfiPi  rjO 

if-' !de[!eo,  ?■'  ingr =4 ' 

ZIWI= JOE*  * ( 10.0**106 1 
12=1  Zf1IN7ZINlUo.  0001  ‘ZMIN/ZINT 
ZZ!1IN=IZ*ZINT 
Zf1IN=ZZf1IN 
ZZ(1AX=ZMIN 
38  ZZttAX=ZZf1AX+ZINT 

IFIZZMAX.GT.  ZflAX)  GO  TO  89 
GO  TO  38 
89  Zf1AX=ZZMAX 

ZPLOT= { ZMAX-ZMIN ) /5. 

GO  TO  94 

ENDIF  ,  ,, 

IFTrEFEL.NF  ’  ’  THEN 

0PENU3.flLE=’lH4334M 
URITEU3, MAAOV )  REFEL 
REUIND  4$ 

READU3. *.ERR=6777]  ZHIN.ZHAX 

fljj  Nj.ST«»S=Uimf) 

418  CONTINUE  .  . 


IFfJKASK.NE. 


ENOIF 

CALL  CLEAR  ,  , , 

Pm 5  2£fEL!0ZEF£L 

IFlZEFEL.EO.  ’THEN 
ZPL(jij=|ZMAX-ZHIN)/5. 

REUI  O  J  A4°  ZEF^L 

MMPi 


Y  ZEFEL=  S 

Y  1  THEN 

MANY  Z  UNITS  PER  PLOT  INCH  (OEF=AUTO)  >  ’ 
ZEFEL 


SfflHK- 


GO  TO  4 18 
ENDIF 


AX-ZITINWZPLOT  ,6T .5.05)  THEN 
LEAR 

^USING  , ZPLOT,  WILL  PLOT  TOO  TALL.  TRY  A  LARGER  VALUE.’ 


CALL  CLEAR  ,  ,  , 
IFfJKASK.NE  'Y  )  ISHA0E=’ 
IFIJKASK.EQ.’Y’j  THEN 
PRINT  MA})  ,  (JANT  PROFI 
READ  ’  A1  \lSHADE 


IFTISHADE.EO  ’s’ 

oshade=ishad£ 

CALL  CLEAR  ,  , , 

iniKA^^l’y’j  ^ 

PRINT, ’  AT]’,  AXI 


PROFILE  SHADING  REMOVED  Y/N  ?(S  FOR  SAME)  >’ 


ISHADE=OSHADE 


LABELS’ ETC 

IF(LEVOFF.EQ.’S’)  LEVOFF=OEVOFF 

OEVOFF=LEVOFF 

CALL  CLEAR 

8855  PRINT*. ySYFIBOL  ON  PROFILE  LINE 


S  LABELS,  ETC  LEFT  OFF  Y/N  ?  >  ’ 


|nT^SYM90L  ON  PROFILE  LINE  (OEF=NON£,  S=$QUARE 
ADT,(iir  ^ymt=trian^le'  °"0I^0N0'  ] 


PRINT,  A 
READ  MAl 
CALL  CLEA 

ifTsym.eq 


,  C=CIRCLE’ 


IFTSYM.EQ.  .QR.$YM,EQ.  $  .OR.SYM.EQ,  C  .OR.SYM.EQ.  T  . 

+CA  L  OR^SYM.EO.V.OR.SYM.EO/*')  GO  TO  8866 

PRINT  V  ’.SYN,’  NOT  A  VALID  SYMBOL. ...  TRY  AGAIN’ 

GO  TO  8855 
8866  CONTINUE  ,  , 

IFTSYM.NE.  ’)  THEN 

PRINT^/lf'’,  ’  LEGEND  FOR  SYMBOL  ?  (20  CHARS  MAX)  )  ’ 

READ  (A20  ’  LEG  ,,, 

URITE(GEL.MA20,Al)’l  LEG, SL 
IFl  JKASK.80.  1  THEN 
CALL  CLEAR,,  , 

PRINT  MAT  .  LEGEND  PRINTED  ON  LEFT  (DEF)  OR  RIGHT  ?  L/R  >  ’ 
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400 

412 


402 


RFAf:  ■ ! A 1  '  JR 
IFILF.EO. 

END  IF 
ENDIF 

CALL  CLEAR 
ZDIF=2PL I”  . 

IF  ZDIF.E0.0,,, 
IF  JKASK.EO.  V 
IF { JKASK . EO.  Y 
CALL  CLEAR 


!(?=■! 


ETURN 
ILEG-0 
CALL  UERASE 


smICL“CK1 

rn  *  Y  )  GO  TO  412 


„  ..Mi 

Si  « 


GO  TO  400 
+Y?I)**2)+XSET 


DIF ) +YSET ) 


IF )+YSET ) 


ur^«ir^w>xsET 


IZ1.EQ.-999.0R.IZ2.E0.-999)  GO  TO  402 


. _  _ SET+YSETj/2. 

IF^0RTl?fGlLpt-^l!^LlG2LAST-G2),*2)?GT. 
^fXYHB0L((G!+G1LAST)/2.  ,  (G2+G2LAST) ?2. , 

G1LAST =G1 
G2LAST=G2 
18=  I8-MSYP1*  5 
'SYM=0 

SHADE. NE.’Y’ 

fm\\u 

bkmff 


25. 

SYM 


Wl. 

CALL 


THEN  , 

m] 


CALL  UC 
SQ=SGRT 


EN 


.  _  ,Y(NJ))**2) 

MLjA»iioY;’f?fT((x(Nj)-x(i),‘‘2t(Y(Nj)' 

rAI  I  I  IDDkl T  i  I  *7  M  T  II  ^  DC  41  ’  | 

t) 


CALL  IJPRNT 
CALL  unovr 
CALL  UPEN 
CALL  UMOV 
CALL  OPEN 
ZS=0.0 


REAL’, 

YSET 


THEN 


Y { 1 T ) *  *2 ) 
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Zi=2HIN 
500  Z5=Z5*YGAF 
Zi=Zl+ZPL0I 

CALL  UflOVEI  100,  ,Z$+YAET ! 

CALL  UPRNT  -  Y 

CALL  UflOVEI  .. 

CALL  UPENIXSET 
CALL  UflOVE  “ 

CALL  UPENT 
CALL  UWHER 
CALL  UflOVE 

CALL  UPENT.,  _ 

CALL  UUHEREtR3.R4) 

CALL  unOVETSO+iSET , ZS- YGAP+YSET ) 
CALL  UPENISO+XSET , iS+YSET ) 
IF[Z1.6E.ZHAX)  60  TO  510 
GO  TO  500 
510  CONTINUE 


Zl,  REAL  , 
XSEt-10. .ZS+YSET 
fSET, ZS+YSET I 

yittts.VSEIl 

SQ+XS^T, ZS+YSET) 


CALL  UFLUSH 
“LL  UflOVEI: 


Rl+60. , R2+20. ] 


CALL  UMOVE 
CALL  UPRNT 
CALL  UMOVE 
CALL  UWHER 
CALL  UMOVE 
CALL  UPRNT 
IFtSEC.NE. 

CALL  - 

CALL 


CALL  . 

DPR=0ESC2  , 

CALL  UPRNT 1  DPR,  TEXT 
DPR=SEC2 

CALL  UflOVE  (SORT  (X(l)**2+Yfl) 
CALL  UUHEREIG1.G2 

. . T$i+lo, 

DPR, '.LA 
ORTIX  NJ 
G1.G2 

1-S5..R4-3 
QPR 


*2 ) +XSET , YSET ) 


R2- 

EX' 


f.l2 


'\’,’TE 
,  YSET-70. 

-VY2 ) *  *2) 


XT’] 


2+Y(NJ)**2)+XSET,YS£T) 

. . .  tALL  UPRNT1 

...ll  m  p  it  httv 

I.GT.100. j  CALL  UPRNTljpi+,05  ’INTE’! 
FtDI.LE.JOO.J  CALL  UPRNT 1 (01 REAL 

«  “m* . 


CALL  UHL. 
CALL  UPEN 
CALL  UWHE 
IFtflfl.EO.  1 


CALL 


IF  (fin.  EO.  2)  CALL  UMOV 


G1-35..G2-30. 
61-15. , G2-30. 


600 


CALL  UPRNTltM  ,  TEXT’) 

Inl.EQ. l)  CALL  UPRNTlj’  OV’.’TEXT’), 
fcli)  CALL  UPRNT1(  L\  ■  TEXT  ) 


RX=R3+2Q 

RY=R4-30 


READtYLABEL, ’ (40A1)’ )  (YL(I ) , 1=1,40) 

??(YL?li'N?:l’;|  GO  TO  851 
CONTINUE 


850  ..... 

1-40 

851  CONTINUE 
DY= ( RY-YSET+ 
IFtDY.GT.Sp 
CALL  US£Trv 
CALL  UP' 


SOFT’) 

’  IIC D  T  * 


/ 1 

=30. 


UPSET ( , VERT , , 0. 75*DY*PSFT 
CALL  UPSET (  HORI  .0. 50  DY'PSFT 


f1 


UR 

IF 


Wi&i; 


otW{P 


700 


CALL  UflOVE ( RX ,  RY 1 
QAL^UPRNTltYYO,  TEXT’] 


tY-OY 


CALL  UPSET  VERT  .30. *PSFT 
CALL  UPSET  ’HQRI’. 20. *PSFT) 
IEiSYM.NE. ,  1  THEN 

( ILEG  +  1 ) 


.  (5YM.Nl  , 
IF  LEG.NE  ’ 
YLEG=R4-3&. 
ILEG= ILEG+1 
XXSS=XSET 


THEN 
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IF-'LR.E'J.  'R' )  XXS$=XXSS+<00. 
CALL  XY«BOL f XXSS+50. . YLE6+10. 
CALL  UhOVE • XX5S+75. .  YL£6 ) 

CALL  UPRNT1 (GEL.  TEXT  T 

ENOIF 

ENOIF 

CALL  UMOVE (0. . 0. 

CALL  UFLUSti 
CALL  USETrSOFT’) 

RETURN 

END 


SUBROUTINE  FENFI 


V,&h 


ISKIP.XMINX, 


colllioNl?RlmEM6PCX’PCY’^FT’c°LOR^L^CK’PSI2E'LPEN  8 

88R!l8N!fu?feBo/XWMll5l?XU^AX°Vwft^N^YU«AX 

COMMON/DAREA/XUDMIN.XUDMitX,  YUDnIN,YUDMAX,FSOFT 
DIMENSION  XTI260) . Yt(260T, ZTi260T 


YLABEL’40 


DATA  SLASH/’  V/ 

I IPEN=LPENt 2T , 

JFpKlf.EO.O)  60  TO  1900 

CALL  CLEAR 
PRINT*/ ...WAIT..’ 

IFlj?.N£.2l  THEN 

EY=Y1 

jpIPr 11  ,he" 

Y2=EY 
X1=TX 
Y1=TY 
60  TO  123 
ENOIF 

IFlABS(X2-Xl).LT.O. 1. AND.Y2.LT. Yl)  THEN 

X2=TX 

Y2=EY 

XI =TX 

Y1=TY 


- GE. 0.1)  THEN 
THEN 


123  CONTINUE 
ENOIF 
XB1=X 1 

X2 

Y2  , 

K.NE.2  THEN 
33  1=1,1041.2 
IF  XBl.LE.&.g.AftD,: 

IF  XBl.61,0.)  THEN 
XB  =X1+F*  X2-X1) 

pFYUF,|Y2-Yli 

IFIXEI. 61,0.0)  THEI 
XEUX2-FMX2-X1 
YE1=Y2-F*(Y2-Y1 
ENOIF 

DO  334  3=1, MB. 

CALL  CORCuft ( J ) 
IFlXBl.GT.O. 1  THEN 
CALL  AREA  CX  1 J . CX 


(Ml 


CALL  CORCORf 
IFfXBl.GT  0. 
CALL  AREA  CX 
CALL  AREA  CX 
CALL  AREA  CX 

ifTai.le.q. A 

XB 1  =  - XB  1 


.XBl.CY  1 


.  C Y  j  2 } , YB 1 . A 1 
,CY  3  .YB1.A2, 


<  B  !  =  -■'?  i 

ENDIF 

ENDIF 

IF' XEl . GT .0 
fit  i  as? i  r 


CALL  A»EA  C*  li.CX  2  .XE1.CY  1  .CY  2  , 

CALL  AS? A  CX  2  .CX  3  .XE1.CY  2  .CY  3  ! 

CALL  AREaICX!3).CX(1  .XEi.CYfj]  CYtlK 

IFIA1 . LE .0. AND. A2.LE.0. AND. A3.LE.0. I  T 

XE1--XE1 

YE1--YE1 

END  F 

END  F 

334  CONTINUE 
333  CONTINUE 

335  CONTINUE 

»::»! 

X2= -X£ 1 
Y2= - YE  1 
ENDIF 

CALL  USE T ( COL  1 ) 

IFTnn.EO.lJ  CALL  UCRCLE  X1.Y1.5. 
IFtnn.EO.ll  CALL  UCRCLE I X2. Y2, 5. j 


, YE1.A1 1 
! YE1 . A2 
. YE  1 , A3 1 
THEN 


NJ=!X2-Xll/5.  , 

IF  NJ.EQ.O  NJ-  Y2-Y11/5. 

NJslOl 

DO  100  1=1, NJ 
01=1 

QNJ=NJ  ,,, 

XT[l|=xl(^|x2-Xl| 

100  YT  I  =Y1+FMY2-Y1 
2HIN=99.E30 
ZHAX=-ZniN 

JFND=0 

LFN0=0 

L||=| 

LTIHE=0 
315  CONTINUE 

00,210  J=LB8,LBJ,L8C 

IF  LFND?£0.1)  GO  TO  210 

FIND  TRIANGLE  POINT  IS  IN . 

CALL  UALPHA 
D0200  IG=II,nB+II-l 

IF? IG. GT . HB)  KK=HB 
I=IG-KK 

IF^T.HBj.  GO  TO  200 

m  wl Mim  J|:K 

CALL  AREAlcxl3.CX(l).XT  J  .CY  3 J , CY t 1 ) . YT t J]  A3) 
IFIA1.LE.O.ANO.A2.LE.O.AND.A3.LE.O. I  THEN 
CALL  UALPHA 

IFTj.EO.LBB. AND.LT IHE.EQ. 1)  JFN0=1 
rOUNb  IT...,  ,  ,  n 

CALL  CON$t$(XT( J  ,YT(J).Z  (J)) 

IFIZTIJ  .GT.ZHAX  ZHAX=ZT  J] 

IF  ZTtJ  .IT.ZrllN  ZHIN=ZT  J 
il^HAXT i,I-lHB/lO)T 
GO  TO  210 
ENDIF 

200  CONTINUE 

IFtJFNO.EO.il  LFND=1 

210  iyJzjftf.LE. -998.0. AND. ZTI101).LE. -998.0)  LFND=! 

LB  J=N  J 
LBC=1 

IFlZT <1  j .LE. -995.0. AND. ZT( 101 ) . GT. -998.0)  THEN 
LBB=N] 

LBI=1 

LBC=-I 

ENDIF 

LTINE=LTIHE+1 
IF  ILTIflE.  EO.  1 )  GO  TO  315 
88  CALL  CLEAR  , 

IFTISKIP.EQ.l)  THEN 


.  I - ( HB / 1 


. EO.  1 )  LFND=1 
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00  603  1=1 .NJ 
‘.n'.  URT  TP !  23  ’ '  * 
ENOIF'.  , 

TASK-  0 
ZZ(1IN=Zf1N 
ZZf1AX=Zf1X 
7002  CONTINUE 
1900  CONTINUE  , 

IFINH.EQ.2)  THEN 
CALL, PEND. 

DE V=  HP4  !  ,  , 
CALL  PSTART (2) 
CALL  UERASE, 

CALL  UBACKG  CLRBC 


•  I ' .  VTI I ) . ZT 1 1 ! 


300 


5iPfALL  UCHPENlLPEN(lj) 

CALL  UDaREA( XUDHIN, XUDMAX , YUDMIN, VCIOflAX ) 


XHINX, XHAXX, YHINY, YHAXY, SH) 


IFIHH.EQ. 
CALL  USET 
CALL  UPSE 
‘  "L  UPSET 


,2A.‘PSFT’FS0FT 
, 12. *PSFT*FSOFT 


.Yl-10.  SIN 
,h-30.*SIN( 
),Yl-60.  SIN 


-20.  *PSFT,YS) 


UPRNTlI  ABC TEXT 
FTYPE.EO  ’N’)  THEN, 

=atan5iy5-yi 


CALL  UHOVE 
CALL  UPENI 


X2+ 

2+31 


n'M  S(ANG),Y2+10. *$IN(AN6) ] 
0. ‘COSTANgT. t2+30. ‘SINlANGl  , 
50.  *COsTaNG) , Y2+50. *SIN? ANu I 


CALL  UM0VE(X2+50. fCOSfANG) , Y2+5G 

^XV^oMiH^L  U(l0VE(XS-20. 
CALL,UP?NtllXVC,  INTE'I 

^pUPRNTK  ABC,  ’TEXT’] 

CALL  UH0VE(X1,Y1! 

CALL  UPENJX2, Y2J 
CALL  UHOHE 
CALL  UFLUSH 
QRAU  THE  Pi 
IFJPX.GT.O.i 

XP=0. 

THE=ATAN2  Y2- 
CALL  USET (COL 

K&99 

ZZ 

CALL  UHOVE IX 
CALL  UPEN  " 

CALL  UPEN 
CALL  UPEN 
CALL  UPEN  ,  . 

CALL  UHOVE F XT 
CONTINUE 
ENOIF 

CALL  UHOHE 
CALL  UFLUSH 


'PSFT.YS) 


PX 

URN 


X2-X1) 


eWII^MIIJhe)! 
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v0LUf,£  under  th£  n£sH  ?  !y/ni  ( ’ 

IFIPASK.EG.  V)  then 

455  ELEVATI0N=  ’•znIN 

'  i’  (2.24c)  ENTER  REFERENCE  ’’BASE”  ELEVATION  (DEF=0.0)  > 


AiOT’.REFEL 
L.EQ.  ,)  REFELE=0.0 
L.NE.  I  THEN 

rti  c.’m/ / 71/ ’  \ 


IF  REFEL.EO. 

W 


n«im 

REMIND 

READU3.  ’'■ERR-44$5)  REFEL 
CL0SEI4S,§TATUS=  DELETE*) 


ENDIF 
AT0T=0 . 0 
VOL-O.O 
VOL1=0.0 

D0Ly'l-HBL,«B 

att  m  i:  i,c 

A=Sf1  Sfl*ABS(A)/$. 
AT0T=AT0T+A  .... 


CX(2) ,CX(3)!CY(l),CY(2),CY(3),A) 


ZAV=  (CZ(  1 HCZ  21+CZ  3 
VOL=VOL+A*  t  ZAV-REFELE 
VOL 1=V0L 1+ A*TZAV-ZniN 
V0L2=V0L2+A  l  ZAV-ZflAX 


REF={-V0L1 * (ZhAX-ZHIN) / ( V0L2-V0L1 )  +ZHIN 
ENDIF 

XDE=X0AX-Xf1IN 


yhmim:i*w 

8yiRl!l:188:i;p®:Yx!lI^ 

0XLEN=  I  XtlAX-XfllN)  *$M 


OXLEN=  ( XMAX-XfilN)  *SI1 
OXHT=(YHAX-Y(1INrSh 
IF  IF IMESH.  NE .  HlflESH ) 


Pall  uiACKi(CLRBCK) 

HIM|SH=FIMESH 

PaPlUDAREAI  XUDHIN,  XUDflAX,  YUDHIN,  YUOflAX ) 

IFlttH.EO, 2)  THEN 

CALL.PEN&, 

QEV=  IBHH  (  , 

feta11  , 

Pall  MlMT' 

®fU0UTLN( 

CALL  UDAREAt XUOMIN, XUDflAX, YUDtllN, YUDflAX ! 
DRAU  POINTS 
60  TO  1639 


1938  CONTINUE  _ 

CALL  uIaCK6(CLRBCK) 
CALL  UOUTLN 

1939  CONTINUE  , 

CALL  USET ( C0L2 ) 


CALL  USET(COLOR) 
XPUXMMN 
XPR=XHMAX 
YPL=YMMIN 

ypr=yhHax 

CALL  CLEAR  .  , . 


OOH.EO. 
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PRINT.  :(A  U  AN  T  TO  (CLIP  OR  'ZIOOM  AREA  TOEF=N)  ?(C  OR  Z>  > 
READ  (Air  .ZOOM  ,  ,, 

IF i ZOOM. EQ.  C’ .OR. ZQOM.EQ.  Z  )  THEN 

PRlST^lAf)’  ’  WINDOW  FROM  (F)ILE  OR  (C)URSOR  (DEF=C)  ?  ’ 

fftPoR^^GO  TO  8346 


9367  ^NT^EN 
PR1 NT,  (A\ 

OPENfAiJII 
GO  TO  9$47 
9346  CALL  CLEAR 
PRINT, ,IA\ 

fflt  M 


FILE  WHICH  CONTAINS  COORDINATES  OF’ 
LEFT  CORNER  AND  WIDTH  OF  WINDOW  > 

, ST  ATUS= ’ OLD ’ ,ERR=9346) 


j’^FILE’.FORCFILE,’  not  found  ...press  return  > 


GO  TO  9367 
CONTINUE  . 


REA0U3,*  XPL ,  YPL ,  WIDT 
CLOSE  4$)  ’ 

XPL={  XPL-SUBX-XfllNXl / 
YPL=  YPL-SUBY-YMINY  I 
WIDT=WIDT/SM 

xpr=xpl+widt 

YPR= YPL 
GO  TO  8345 
8346  CONTINUE 
CALL  CLEAR 

PRINr,VUSE  CURSOR  TO 
CALL  UflOVEUOO.  ,400. T 
CALL  KURSINlXPL  YPL, IF 
CALL  UMQVE(XPL-iO.,faL 
CALL  UPENjXPL+lO. .iPLj 
CALL  UMOVEIXPL. YPt-10. 
CALL  UPENlXPL,YPL+10.: 


SM  +100. 
SM  +100. 


:  CURSOR  TO  LOCATE  LOWER  LEFT  CORNER  OF  AREA’ 
itxPL’ YPL' IP) 

XPL, YPd-10. ) 


CALL  CLEAR 

PRINT  V  USE  CURSOR  TO  INDICA1 

P|jt  SWRlP) 

lP?A^sVxPR-XPL).GT.0.001)  YPF 

IF  ABSYPR-YPL.GT. 0.001  XPfl 

XWMIN=XPL 

XWHAX=XPR 

YWniN=YPL 

YUMAX=YPR 

PlSoN.EO  *C’)  THEN 
XUDMAX=  (XPR-I.)M65./1039,  + 

P®  (ft  8:1: 118: : 

CALL  UDAREA  XUDfllN,  XUDMAX ,  YUC 
GO  TO  1938 
ENDIF  ,  , , 

IFIZOOM.EQ. 'Z  )  THEN 

CALL  UWlND6(XPUXPR.YPL,YPR) 

FS0FT=(XPR-XPL)h03i 

GO  TO  1938 

ENDIF 


TO  INDICATE  WIDTH  (OR  HEIGHT)  OF  AREA’ 

.IP) 

0.001)  YPR=YPL+ ( 780. /1039. ) * ( XPR-XPL ) 
0.001)  XPR=XPL+l 1039. /780. !  lYPR-YPL) 


XPR-XPL 

YPR-YPL 


EO.’C’)  THEN 


mmm 

MIN,  XUDMAX,  YUDNIN, 


YUOMAX) 


EallFupset  (;vert;,3o.‘fsoft*psft) 

M 1  Wi X ,fS0F,' Wl1 

IF(KW.EQ.l)  THEN 


RDASK, 0, XMINX, XMAXX, 


ENDIF 

886«  CAi  L  CLEAR 
XBX=QXBX 
YBX=OYBX 
XL£N=OXLEN 
XHT =OXHT 


ANGG=OANGG 
PR  NT‘,;(2.41a)’ 

PR  NT*,  ENtER  PROFILE  GRID 
PR  NT*,  BOX  TO  HORZ.AENGTI 
PR  NT  , 'OR, ENTER  FOR 
READ  tA60)  ,PIF 


BOX  ORIGIN  OF  COORDINATES.  ANGLE 
H  OF  BOX,  &  HEIGHT., U.Y.ANG.L, HI 
SAME  AS  LAST  TIME  (DEMUTO) 
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r I  r  AR 

:  ;j  GO  TO  3828 

ItfuPM'  - 

84. *1  XBX.YBX. 


GO  TO  8868 
ENDIF 
OPEN (85 .FJLE=H33 
URITE(85/[2X.A60 
REWIND  8$ 


u-  .  j.  i  uL’ .STATUS=' OLD’ . ERR=S869 ! 
XBX . YBX , ANGG, XLEN. XHT 


PIF 


OPEN 
WRIT 
CLOr 
CLO. 

GO  T 

8863  OPEN, _ _ 

READ(84  *.ERR- 
CLOSE(8i 


ERRf 8863 j  jf^X, YBX, ANGG, XLEN, XHT 


=  PIFF 


■FlSIirV  XBX, YBX, ANGG. XLEN, XHT 

o  alIiTATUS=’DELETE,) 

J84|FILE-PIF,$TATUS=’0LD’ . ERR=8869 ) 
184.?. ERR=88&9 )  XBX.YBX, ANGG , XLEN , XHT 


8828 


CLOSE185, ST ATUS=’ DELETE’ 
OPEN(84.nLE-  PIFFEL.y) 
WRITE[8i/(5F15.4rI  XBX, 
0SE(84) 


Y8X, ANGG, XLEN, XHT 


♦  100. 

SMj+100. 


CL...,. 

8868  CONTINU, 

CALL  CLEAR 
XAM  XBX-SUBX-XHINX 
YA=l|YBX-SUBY-YniNY 

P2=YA  , 

ANG=ANGG/57. 29578  ,,  , 

P3=(TXBX-SU8X-XHINX+XLEN  /$H  UOO. 
P4=P2 

6=UYBX-SUBY-YMINY+XHT)/SM)+10Q. 


BX, YBX, ANGG, XLEN. XHT 

8855  I  X=,F12.2.3H  Y=,Fl2.2,5H  ANG=,F6.2,3H  L=,F10.2, 


.  REDO  BOX 

.Etype 

THEN 


?  (Y/N)  > 


19 


5578 


GO  TO  8869 

P§fST*,’(2,42)  WANT  PARALLEL, (TO  BASE)  OR  NORMAL 

«bt'  { 1  p " 1  1 

CALL  CLEAR  . 


GRIDS  ? 


READ  !A1., 
IFjEXAG.EO. 
IFlEXAG.NE. 


A 


m 


AND.FTYPE.NE.’N’I 
2.43]  HOW  MANY  PR 


GO  TO  19  , 

OFILES  ON  THE  PAGE  ?(DEF=11)  > 


MANY=11 
THEN. 


rnmvm 

REWIND  4i 
READ  [43.  . ERR= 5578 )  MANY 
CL0SEU3, STATUS:  DELETE') 

px=minIxl!n^h!x^/sm)/many 

CALL  CLEAR  ,  ,  ,  , 

PRINT  MA\  (2.44)  ENTER  SCALE  EXAG.  FACTOR  (DEF=1.0)  > 
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READ  ( A ! 0 j  .  EXAG . , , 

READIEXAG.  F10.ni’  1  EX 
IF l EX . EO. 6.  EX:  1.0 
IFlEXAG.EO.  0V!  EX^O. 

J  F — 0 

IF?EX.GE.O. 1)  THEN 
3366  CALL  CLEAR  , 

PRINT*  MIN  EL=  \ZMIN.  MAX  EL=  .ZMAX  ’  FOR  PROFILE  SCALES 
PRINT.  (A\)  .  ENTER  NEW  MIN  4  MAX  EL  S.  (DEFEND  CHANGE!  ! 


IS  :;  #5E  - 

IFTNEUZS.NE.',  I  THEN 

skihs  (“°”  ,E"2S 

READU3  *  ERR=3366)  ZMIN.ZMAX 
£LgSpJ,$TATUS=mETEM 


mi'n 

IF ( MANY . EO. 0 )  THEN 
CALL  FENFILlZMAX.ZMIN 
+  YMINY. YMAXY,SM,6.FTY 

aw 

ENDIF 


N.MM,X1,X2,Y1,Y2,EX,RDASK,0.XMINX,XMAXX, 


^NT^flHE  6RIDDED  PAGE  PLOT  FILES  ARE  NAMED  WITH  THE  SAME’ 
PRINT*,  ROOT  AS  THE  .MSH  FILE.  EXTENSION  ARE  .OP.  ,1P.S  , 

PRINT,  (A\) '.’ON,. IN,  ETC  PRESS  RETURN  TO  CONTINUE  )  ’ 

pADfFlMESH,  I3A1)’ )  (Ff  J),  J=l,8) 


HaJ-i 
GO  TO  27 
ENDIF 

277  CONTINUE 

278  CONTINUE 


. OR . F ( J ) . EO . ’ 


®I>F '' THEN 

M=M+ 1 


(Ef J) , J=l. 12) 


177  W^iWnME  ’(12A1)’] 


IF  MANY. 
IF  I  MANY. 
XlsXBX-S 
X2=XBX+X 


.  I  I  I1LH. 

1  FO=.XHT 
1  FD=f XHT 


(MANY-1) 


X2=XBX+XLEN-SUBX 

Y1:(YBX-SUByTmM+JF)*FD 

SI  itsp  MM: 


XA=XBX-$08X  , ,  . 

YA=((YA-YHINY]/SMMOO. 

CALL  ROTATEIXI.XA, Yl, YA.ANG 
CALL  R0TATE(X2.XA,Y2,YA,ANG 


IF^TYPE.EO.  N  )  THEN  % 
IF  MANY.NE.l  FD=  CXLENl 
IF ! MANY . EQ . 1  FD= l XLEN J 
XI  =  f XBX-SU6X  +  I-l+JFJ* 
X2=X1 

Y 1 =YBX-SUBY 
Y2=YBX+XHT-SUBY 


(MANY-1) 
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*!  =  •  >i-Xf1INX.t /SMJ  +  100. 
n- f  f  X2-XMINX I /SMJ  + 100. 

Y!=>  i  Yl-YdlNYl/Sfl  +100. 

Y2=  ( >  Y2-YMINY} /Sfl!  + 100. 
ANG=ANGG/57. 29578 

XA=i -  - 

X  A= ! 

Y  A=Y  BX -SUB Y  ,  , 

YA=  ( (YA-YMINY ) /Sfl ) +100. 

CALL  ROTATE  XI, XA,Y1,YA,AN6 
CALL  ROTATE  X2  XA.Y2  YA.ANG 
ENOIF 

VX=SM*  Xl-100.  UXdlNX 

yy=su* * Yi-inn  >avmtwv 


a=I!xa-v8!nx)/sii)+ioo. 


VY=SH*(Y2-100. ) 
WRITE  23/1  VX. 
PORN-  v'  ,  , 
IFlFTYPE.EQ.  N 
IF  FTYPE.EO.  P 
IFlFTYPE.EQ. ’N’ 
WRITE  23.  IA22, 
WRITE  23,*}  4 

Qycj^-f  IIITC6  TAB 


CO.  +YMINY 
....  vy 
+XHINX 
+YMINY 
VX.  V  Y 


PORN: 
PROFI: 
PROF J= 

a. Air) 


GRID  PROFILE 
GRID  PROFILE 
PROFI. 1-1, PORN 


(A22J 


ELEVATION 

GWEL 


HR... 

IVC=I 
CALL  CLEAR  , 

CALL  FtNFIL(2NAX,Z“IN^HN^X| 


YHINY.Y! 
CLOSE  1 2 J ) 
77  CONTINUE 
79  CONTINUE 
IFIPASK.Ei 
IF  ZOOM. Ei 
IF  l  ZOOM ,  E1 
CALL  USE! 


flAXY.SPl,  IVC, 


SOFT 

(  ’  uro 


,X2,Y1,Y2.EX, ROASK , 1 . XMINX , XMAXX , 


(10/PCX,  75/PCX.  10/PCY.75/PCY) 
!l.,1040!,0.,780:j 


CALL  UPSET  ! VERT  ,30/PSFT 
CALL  UPSEU  HuRI  ilS-  PSFT 
CALL  UflOvf--  -1 
CALL  UPRNI 
CALL  UPRN 


sd 


CALL  UPRNI 
CALL  UPRN" 
CALL  UPRN 
CALL  UPRN 
CALL  UPRN 
CALL  UPRNT 
IFIHH.EQ.2 
CALL  UDARE.. 
CALL  UWINOO 
CALL  CLEAR 
PRINT*, ,  VOL 
PRINT  .  AV  Z 
CALL  UBELL 
CALL  UPAUSE 
ENDIF 

CALL  CLEAR 

IF  "  ‘ 


% 

V0La* 

V.L, , 

CALL ’ UsItV’  HARD’  I 
XUOMIN.XUOflAX,  YUDHIN.  YUOMAX) 
XWhIN,  J'WHAX,  YflfllN,  YWFIAX ! 


ABOyE^REFELE,  ,  VOL ,  ’ 


AREA=  ' ,ATOT 


RE 

L 

IFIKW. 

MU 


TO  DRAW  ANOTHER  GRID  BOX  PROFILE  ?  (Y/N)  > 


1456 


WANT  TO  REDRAW 
,  JASK 


IP 
H 

j£SK 
r  ^ 

read2'?ai 

CALL  UERA 
CALL  UBACL, 

IFIjASK.EQ. 

ENDIF 
ENDIF 
ENDIF 

CALL  uERASE, 

CALL  UBACKG  CLRBCK I 

RETURN 

END 

SUBROUTINE  XYMBOL ( X, Y, SYPl) 
CHARACTER  $YH*l 
!F(SYfl.EQ. ,  RETURN 
IFISYH.EQ.’*’  THEN 


THE  CONTOURS  ?  (Y/N) 


GO  TO  1456 
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CALL  LINE' SYM. 
CALL  LINE  SYM. 
CALL  LINE  SYM. 
CALL  LINE  SYM. 
ENDIF 
IFISYfLEQ 
CALL  LINE 
CALL  LINE 
CALL  LINE 
CALL  LINE 
ENDIF  ,  , 
IFTSYM.EQ.  T 


SYM. 

SYM. 

SYM. 

SYM. 


X. Y. 1.6) 
X.Y.J.ii 
X.Y.2.5 
X.Y.7.S) 

THEN  , 
X.Y.1.3 
X.Y.3.6 
X.Y,6,A 
X.Y.4.1) 

THEN 


LINE($YM,X,Y,’5j  i] 

THEN 
Y.2.5) 


_ M.EO.’C’) 

CALLUCRCLEIX, 


s  Yn.EO, ’O’  THEN 
CALL  LINE  SYM. X,Y, 2,7 
CALL  LINE  SYM, X  Y, 2,8 
CALL  LINE  SYM, X, Y, 5, 7 
CALL  LINEISYMjX  Y,5,8 
'"MF 
TURN 


100 


SUBROUTINE  LINE(SYH,X,Y, I, J) 

»iQN  y(8J.Vi8) 

U=1.0 

IfIsYH.EO  ’0’),U=1.25 

CALL  UMOVE(X.Y) 

|  wk»W!t»WlI)’ 

RETURN 

END 

X.XA.Y.YA.ANG)  % 
-  Y-*A  ‘slNfANG 
.  ♦(x-xaIsin(ang) 

X=XN+XA 

Y=YN+YA 

RETURN 

END 
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A P PENDIX  C.  FORTRAN  LISTING  FOR  McGRID 


HinENSIOM  XI 2000 
DIMENSION  IBB ! 40 
CHARACTER  FILLE’24 .  ISG'l.HONOR’l. ZCHAR*20.  ICBOR'l.FF  14j*l 
CHARACTER  I  ASK  *  I ,  T  YP*  4 .  SL  ASH^l .  ASI S*  1 .  DUM 1 200 .  FUfli  1 )  *  200 
CHARACTER  VORflAT4 1S2. CORHAT’182. V0R1*90, V0R2*^0.C6R1*90.C0R2'9 
CHARACTER  CDUMI2Q  ’20.  IFDUM  (200Pl.DUMMl24 
CHARACTER  VOR  20j‘«.C0R  20 ) *9 , L 1  *  1 . L2*6 . L3 * 2 . OPREN* 1 . CPREN* 1 
CHARACTER  TOR ( 20 1 *9 .Ml *4 .M2*  1 . M3*2. 10R1*90. 10R2*90. IORMAT’185 
CHARACTER  0NEX*3, DE^AUL5^ 

CHARACTER  GXT200i*l . KEEPflO  *20 , NKEEP * 20 . PCDEN *20 . C ID*20 
CHARACTER  CD1200Q T * 20 , LLBL^ 10 , T^ASH^ 1 


.  Y 1 2000 1 . Z 1 2000 1 . 1 C ! 20 1 . ROOM [ 20 ! . ZZV 1 2000 1 


CHARACTER  CDI200 
DATA  ZZV/ 2000* J. . 
DATA  GX/200*  .1 


DATA  ONl)X^)X,/ 
DATA  OPREN /,(  / 

DATA  CPREN/  )  / 

DATA  V0R72O* ; 

DATA  C0R/2Q* 

DATA  I0R/20  ’ 

DATA  CDUM/20*  XXXXX 

PRINT^f^**  *  **** 


xxxxxxxxxxxxxxxxxxx’ / 


|  ««««*  ttttttt  ttttt  XOHtlll  > 


CODE”**’ 

(601)  634-3357..’ 
RIMENT  STATION* 


PRINT*  ’’‘THIS  IS  McCON— A  GENERAL  CONTOURING  CODE****’ 

PRINT*.  . TROUBLED....  CALL  JOHN  PALNERTON  [601]  634-3357..’ 

PRINT*.  . OR  URITE  TO:  USAE  WATERWAYS  EXPERIMENT  STATION* 

PRINT*.  ATTN:  WESGS-Ry 

PRINT*.  GEOTECHNICAL  LABORATORY 

PRINT*.  SOIL  AND  ROCK  MECHANICS  DIVISION’ 

PRINT*.  3909  HALLS  FERRY  ROAD  .  , 

prInt*’  *»«*«********Sfim8!ls*»»!!?1»l§i!;vJ}l******»*it*' 
print*!;  . 

PRINT,  (A\  ’.’  PRESS  ANY  KEY  TO  CONTINUE  )  ’ 
readt;(*d;,iask 

PPTWT*  'PflIITTMF  il  I  nuc  Vflll  Tn  U AMP  AM  TMP1IT  hATA  CTI  t  Akin  TUPIl’ 


3909  HALLS  FERRY  ROAD  .  , 

i»*«»«***»***»»***» 


PRESS  ANY  KEY  TO  CONTINUE  ) 


PRINT*. 

PRINT*, 


YOU  WANT  WAS  Jl 

60  Big&TWJP 


V  IWVIAHL^I  XI  mil  I  ini  uni 

,« tsa'Wfis  ® 

. . NG* 


’’CONTOURING 
NOW  TO  GO  TO  C 


PRINT*  ,  PRESS  C  NOW  TO  GO  TO  CONTOURING’ 

PRINT*,  OR  PRESS  "V"  TO  GET  MORE  VERBAGE/ 

PRINT,  (Aj)'.  OR  PRESS  ”  RETURN”  TO  BEGlA  INPUT  PROCESS.  ) 

IFtlASK.EO'cf)  STOP  ’..GOING  TO  CONTOURING...’ 

TC(T»tY  crt  ’  M  rr\  ta  n 


IF  IASK.EQ, 

PR  NT,  ’ 

PR  NT*,  ..A  WORO  ABOUT  THE  INPUT  FILE.  THIS  ROUTINE  WILL  ASK’ 
PRINT*  'YOU  FOR  THE  NAME  OF  A  DATA  FILE  WHICH  YOU  HAVE  PRE¬ 
PRINT*,  yjOUSLY  PREPARED, ,  THIS  DATA  FILE  MAY  EITHER, BES  , 

PR  NT*.  ”FREEFORMvs  OR  COLUMNAR”.  IF  ”FRtEFORM  ,,  THE  , 
PR  NT*,  DATA  MUST  ALL  BE  NUMERIC  (NO  LETTERS  ANYWHERE)!  EACH’ 
PR  NT*.  LINE  OF  THE  FILE  MUST  CONTAIN  THE  X  i  Y  COORDINATES 
PR  NT’,  Pf  f  DATA  POINT.  THE  LINE  MAY  ALSO  CONTAIN  MULTIPLE 
PRINT*  ,  Z  (OR  ELEVATION)  VALUES  FOR  EACH  DATA  POINT.  IT 
PR  NT’,  MAY  AL$p  INCLUDE  f  POINT  IDENTIFIER  COLUMN  TNUMERIC) 

PR  NT*.  AND  A  yvSELECTOR  ”  COLUMN  (USED  TO  SKIP  OVER  CERTAIN', 
PR  NT*  ’DATA  POINTS).  THE  ROUTINE  PRINTS  THE  FIRST  FEW  LINES’, 
PR  NT*,  AND, THEN  ASKS  YOU  TO  SELECT  THE  MEANING  OF  THE  COLUMNS’ 
PR  NT,  (A\)’.v  PRESS  ANY  KEY  TO  CONTINUE  >  v 

pIint*  ^  ASK 

PRINT*!, ’..IF  THE  INPUT  FILE  IS  ’’FREEFORM”:’ 

PRINT’  THERE  MUST  BE  THE  SAME  NUMBER  OF  ENTRIES  ON  EACH  LINE 


GO  TO  17 


PR  NT’, 
PR  NT*! 
PR  NT*. 
PR  NT*. 
PR  NT*. 
PR  NT*. 
PR  NT*. 
PRINT*, 
PR  NT*. 
PR  NT*. 
PR  NT*. 
PR  NT*. 
PR  NT*. 
PR  NT*. 


..IF  THE  INPUT  FILE  IS 


’FREEFORM’ 


..IF  THE  INPUT  FILE  IS  ’’COLUMNAR”:’ 

Sb)p8lET^k?‘)SEc8XY:N2oSF!TL!EKE.R!GK,-JUS,IFIEO 

,  Hi  ili !?  l5‘  ?  ,8* 

IF  THE  FILE’ 


CONTAINS  F 
LINE  ABOVE 


IELDS  THAT 
.  THEN  THE 


IGNORED  IF  YOU  ARE 


T  ARE  BLANK  !5TH  COLUI 
E  INFORMATION  IN  THAT 
DOING  AN  OPERATION  TH< 


ION  THAT  INVOLVES' 


PRINT',.  T 
PRINT’.,  t 
PRINT’, 
PRINT,  I A\ I 


THAT  PIELD  i LIKE  SPECIFYING  THAT. THE  STH  COLUMN  I?' 
JHE  ELEVATION  TO  EE  CONTOURED'. 


read  ’(*1 mask 

PR  NT’.;  ’ 

PR  NT’,  ...ONE  MORE  THING . ’ 

PR  NT*.  IN  THE  EVENT  THAT  THE  IDENTIFIER  (WHEN  SPECIFIED)  IS’ 
PR  NT;,. REPEATED,, ONLY  THE  LAST  LINE  WITH  THE  SAHE  VALUE  I$? 

PR  n}*-’R"AINE0' 

PR  NT  '  (A\)’ . ’  PRESS  ANY  KEY  TO  CONTINUE  >  ’ 

f  i  " 

PRINT*  ! 

PPlMT*  ’ 

PRINT  '  (A\)’ 

Ii  rutrn  T  rs  *  t  a  rn  r  tnur  l  r.  r  r  _  u  xr  u  a  *  »  *  \  \  * 


PRESS  ANY  KEY  TO  CONTINUE 


PRESS  ANY  KEY  TO  CONTINUE 


+R|T9llSlfeRIGIN?L!  °ATA  , 


FILE  NAME  ( DEF=MCDEMO. DAT )  > 


6  PRINT*  ’FILE  ’.FILLE.’  NOT  FOUND. ...TRY  AGAIN.’ 
GO  TO  1 7 

7  CONTINUE  .  .  .  . 


OPEN!  10,FILE=’  TEtl’ ,$TATUS=’ OLD’ . ERR=7390) 

CLOSET l6, STATUS:' DELETE  T  ,  ,, 

D  OPEN|lPjRM:vTEHV,STATUS=’NE«’ ) 

*[ Jl.Olbl'ls’FILE  ’FREE-FIELD’  OR  'COLUMNAR'  (F(DEF)  OR  C)  ?  ) 

W-H\  * >  f'EF0,=1 

+  11.02]  UANt  INVERSE  POWER  GRIDDING  ?(Y,N  OR  ?)  >  ’ 

?eyM*-!s* 

ffiP;50,:?  1 THEN 


R;Bi: 

PR  NT* 

a  SI ; 

H  n:  ; 


THIS  OPTION  WILL  CREATE  A  RECTANGULAR  GRID  OF  DATA’  , 
POINTS  WHICH  WILL  ENCOMPASS  THE  ORIGINAL  POINTS  (WITH  , 

A  SMALL  BORDER],  THE  CONTOURABLE  VALUE  (ELEVATIONJ  OF’ 
THE  GRID  POINTS  IS  DETERMINED  BY  THE  INVERSE  POWER' 
AVERAGE  OF  ALL  ORIGINAL  DATA  POINTS.  THE  ORIGINAL  DATA’ 

IS; , 

IT  IS  NOT  RECOMMENDED  FOR  CASES  ilHERE  YOU  HAVE  GOOD* 
SPATIAL  COVERAGE. 

THE  FORflULA  FOR  INVERSE  POWER  AVERAGING  IS:’ 

Zs.UJrSUMt  (Zo(i)/(1.+P.(ij])’*0] . /.SUMll.  /(l.+D(ij))*;0]’ 


PRINT* 

Bin; 


z8S1i!esz2}! 

ml 


*0J  /  SUM 
LEGATION  i 


0  IS  THE  POWER 


»LEVAlI 

W 

LY  2) 


ION  AT  GRID 
AN  ^RIQINAL 


flUEf  5’ 


KBOR=’N’ 

0=3 

A$I$=’N’  ,  , 

IF  ISG.pQ.’Y’.OR 
PRINT*.  Q3]  u 
PRINT,  (A  V  0 

£  *5  S.IE.'y'.* 


wan?g^?giHa^  data°p 

OUT  TO  EDGES  OF  PLOT 


T^jAS  JS'  ?  (NO  FILLING  ’ 


G.EQ.’Y \6R.ISG.EQ.V)  THEN 


TER  VALUE  FOR  Q  (POWER  FOR  INVERSE  FITTING)  DEF=2  >’ 
, )EFAUL 

tDEFAUL.NMBER.2) 


PRINT  ’ (A\ I ’ , 

*REA^?I  )^?ASKAPR0N  ^  AR°IJN°  °ATA  P°INTS  ?  !Y/N!  } 
ICBOR=’r  '  K 

IFIIASK.EQ. ’Y’ .OR. I'SK.EQ. ’v’ )  ICBOR=’N’ 


C2 


I80R=0 

IF ! ICBOR, EG.  Y  .OR. ICBOR. EG.  y  )  IB0R=2 
HONOR;  Y 

IF >  ISG. £0.  Y  .OR.  ISG.EG.  v  )  THEN 
18  PRINT  ’ [A\ I  . 

*  11.06)  DO  NOT  HONOR  ORIGINAL  DATA  POINTS  (Y.N  OR  ?  )  ’ 

READ  ' [All. I ASK 
HONORS  Y*  ,  , 

IF  IASK.EO.  Y ’OR, IASK.EO. ’ y  !  HONOR:  N 
IF  IASK.EO.  ?  )  THEN 
PR  NT*.  ’ 

PRINT*,  ORIGINAL  DATA  POINTS.  IF  HONORED.  WILL  BE  ADDED, TO  THE' 

PR  NT*’, ;tnOO?H?NG°OFT^AKSKANDf1v|LLiYSIUHICHL^FT^NARESULTS  FROtl’ 
PR  N|*, J INyERSE  POUER  AVERAGING. 

GOTO  i8 
ENDIF 
ENDIF 
ENDIF 

VORMAT:  * 

CORflAT:  (*]  , 

IF(fREFOR.EQ.O  GO  TO  8447 


CORf1AT=  ( 
IFIFREFOR 
REWIND  0? 
8870  READ109, ' 
READ  DUfl, 


DO  8878  1=1,200, 
pUfl(I)  NE. 

BNUE 
8879 
NUE 
D  09 


1,1=1,200) 


^Dunlf):Gx!l|° 

IMAM  ,  gx  i  ,1=1,80 


8450  CONTINUE 

8451  lLASt=I+l 


DO  8460  1=2, IL AST 

ibSM-i  i  E0'  •AND-IFDUf1(I-1,-N£- 


ENDI^ 

8460  CONTINUE 


I 

L2=’.0.1X 

NUliBMBBl 


7  1=2, J 


,L2:’,0 


NUHB= IBB (}) -IBB t 1-1) 
IFfI.GT.2)  NUf)6=NUHB-l 

ur1teTl3,  (12)  r  numb 

V0K[I;1):L1//L3//L2 

I  1  -  f  A  * 


in  i, 

?  Mm 

lltf-2'  * 


iiiiiFi 


8777  CON 


V0R1=V0RH 

COR  1 =COR 1 1 
+  COR 


COR? 
♦  COR! 


j4j[/yoR(5)//V°R(6)// 
VgR(14)//VOR(15)//VOR(16)// 
U  //C0R(5)//C0R(6)// 
/CMll4)//C0R(15)//C0R(16)// 


m m  nm  imp  $  5 1!6  1 

OR  H WM (l? WM  1J  ({IOB  16 
IORf 17  //IOR  181//IQR  l^l//IORi 20 
=OPREN/  V0R1//V0R2//CPREN 
:OPREN//COR1//COR2//CPREN 


f/I0R(5l//I0R!6)// 


V0R(1AT=0PREN//V0R1//V0R2//CPREN 
CORHAT:OPREN//COR1//COR2//CPREN, 
I0RMAT=0PREN//0NEX//I0R1//I0R2//CPREN 
8447  CONTINUE 


.  //IQR(15)//I0R(16)// 


C3 


PRINT’ , 

PRINT’.  ..FIRST  FEU  LINE?  OF  DATA  FII  E  .FILLE.  LOOKS  LIKE... 
PRINT’. 

IFiFREFOR.EQ.il  PRINT  IQRflAT.  1 1. 1=1.  J-i l 
PRINT  M1X.79A1  MGX  I).I  =  1.7<!l 
REWIND  J9  ^ 

READ^n?  "1a2GQ),’ , END=667 )  DUH 


6S6  RW5Mi‘:EBln 

667  CONTINUE 


8899  KEEn??)^"1’ 10 


ffifV;1’  THEN 

IS  'I0ENTIFIER'  ■  (0  IF  N0NE’  DEpr°'  )  ' 


iwteL(ota’N"BER'o) 

0EN=O 


M=9(a\]’ 

V  (1.07b). DO’ NOT  REMOVE  DUPLICATE  LABEL  ENTRYS  ?  (Y/N)  )  ’ 


|llIAS^Ei.;^SOR.IASK.EO.’Y’)  IELIf1=l 

08 rujltH  COL  IS  "SELECTOR'  ?  (0  IF  NONE,  DEF=0)  )  ’ 

READ  IA5 J  ,DEFAUL 

!ALL  DEFALL  t  DEFAUL , NMBER , 0 ) 

S  LsNNBER 

8  S*7TIa\)?!  ,Tj21?o?)  enter  ist  selector  "keep  value" (def=cr)  > 


EAO  ’ ( A20  j ’ ; KEEP ( 1 


ii  i ixul.1 


READ  (A20J  , KEI 
GO  TO  6693 
END  F 
lh  J-l 
END  * 

END  F . 


1,09)  ENTER  NEXT  SELECTOR  "KEEP  VALUE" (DEF=CR) > 


BST-IJlj'ltolk"1  ""ICH  “l  ,s ?  (“EFs1’  ’ ' 

'MICH  C0L  IS  'Y‘  ?  (0EF=2)  )  ’ 

READ  f A5 J  , DEFAUL 
CALL^FaCL [DEFAUL, IY, 2) 

*1  114?  WlfU,  COL  IS  "21"  ?  (-)SIGN  IF  DESIRED  (DEF=3)  )  ’ 
READ  lA5j  .DEFAUL 


1.10)  WHICH  COL  IS  "X"  ?  f DEF=1 )  >  ’ 


L ill' JnlCH  COL  IS  "Y"  ?  (0EF=2)  >  ’ 


W  If  15  ‘2;‘ 

^#(U™4Z1,3) 


CAUDf  A  UDEFAUL.^U) 

4  (1.13  UHlfcH  COL  IS  "22"  ?(FOR  21+0R-22)0R  OfDEFlWHEN  DONE  > 
READ  LA5j  .DEFAUL 

?(2§  W'IZMI 

PRINT  ’(a\7\  .  . 

4  l.)3l  WpitH  COL  IS  23  ?(21+-Z2*-23)0R  O(DEF)  WHEN  DONE  >  ’ 
cIlL  Dt^ALLllipAuL IZ3,0) 


PRINT  ’lA\?  .  fl.  13)  WHICH  COL  IS  "24"  ?  (DEF=0)  (ETC) 

BF 1A  1  I  i C  i  '  AFFiln 


^UfOEF  ALL  l  DEFAUL  .124,0) 

K  WHICH  COL  IS  '25'  ?  ( DEF=0 )  ETC)  >  ’ 

READ  (A5  .teFAUL 
^A^DEFALLIdEFAUL.  125.0) 

1 2f1AX=f1AX  ( ABS  ( IX) ,  ABS ( IY ) .  ABS ( 121 ) .  ABS ( I Z2 1 .  ABS T IZ3! .  ABS ( IZ4 ) , 


A88:  !Z5' . AB?I  I SEL  1 .  AE f 1  IDEM  > 


REUIND  0« 
X  C  0  E  N = 0 
U  HINTING 


XCP£N=XCCiEN+l 

IPfpREFOR.EQ.l)  RE  AD  ( 0« .  CORHAT .  ENQ= 1 S .  ERR =3E- » (f  PUm  1 1 .  M.  IZHAX  i 
IF  FREFOR.EQ.O  READCO^.  * . ENO=  1 5 . ERR=3S I  lEDUnr  I  ).!  =  !.  IZflAX  I 


IF  f  FREFOR. EO. i 
00  8855  1=1 . 12 

MMix 


THEN 

MAX 


END 
6655  CON 


IF!I.ES,ABSlIX).OR.I.EOtABS(IY  .OR.I.EQ.ABS(IZ1!.OR 
+ABSrrZ2) .OR,I,EO. ABS  JZ3).0R,I,EQ.A8S!IZ4].0R.  LEG. 
READICDt)mi).’lF20.0r)  FDUI1  I 


OR. IRQ.  ,, 

‘  ABSTlZ5))THEN 


ffir 

60  TO  39 

38  LCDEN=XfDEN 
*  n 


PRINT*.,. TOATA  ERROR  LINE  ’ , LCQEN, ’  PR08ABLY  NON-fHMERIC  DATA’ 
PRINT*.’,, LINE  V.LCDEN,S  IGNORED . 

skis*:1  ^ 


GO  TO  14 
39  CONTINUE 
CID=CDUn  IDE 


NK«P=CDUNlllELfN 
IF(NK^p!1o!^P(II))  60  TO  U03 


7733  mm 

1403  CON - 

END 

cx= 


GO  TO  14 


GO  TO  14 


GO  TO  14 


GO  TO  14 


GO  TO  14 


2000 


0)  CX, CY, CZ, CID 
) 1  \ti\U ( { 0?20bo )  CX ,  C Y .  CZ ,  PCDEN 

,F19.9, 1X,F19.9,1X,F19.9, 1X.A2Q) 


ffjjplil'1* 

19  CONTINUE 
00  77  JrK,14 


C5 


77  FF 1 1 '  = '  ' 
r.n  in  s-, 

END  IF 

58  CONTINUE 
53  ^ONTINUE 

UR  If  8,  (J4A1)’!  (FF( I i . 1=1. 14 ) 

PRINT*  '  r 

FR  NT*',’.  *********** IMPORTANT ;*******J*************************’ 
PR  NT”  OUTPUT  DATA  FILE  IS  MAP. DAT 

PR  NT*.  1 

PR  NT*  'DATA  POINTS  UHICH  ARE  CLOSE  TOGETHER  MAY  BE  ASSIGNED 

PR  NT*.  TO  A  SINGLE  POSITION  (THE  ELEVATIONS  OF  THOSE  POINTS 

PR  NT*.  ARE  AVERAGED  AND  ASSIGNED  TO  THE  SINGLE  POSITION]. * 

gg  ST*.*  • 


PIMP: 


JGRID-I2 

xfliN=99.E30 

XMAX=-99. E30 

YMIN=XMIN 

YMAX=XMAX 

2MIN=XMIN 

2MAX=XMAX 


auit  CHECKING  DISTANCE  BETWEEN  DATA  POINTS  (TO  BE 

B  M  lim  lh 

7  SCREEN  UNITS  OF  EACH  OTHER  WILL  BE  COMBINED. 

NEW  CHECKING  DISTANCE  IF  DESIRED  ( DEF=7 )  >  ’ 

AUL 

AUL.MSET.7) 

WAIT’ 


137  READ110.2000.END=139)  X(I  ,Y(I) 
iFfx  iT.gt.xmaxT  XMAX=X  I 


MUM  SMI 

il.LT.YMINJ  YMIN=Yt I 
.137 


139 

slx=?XMAX-XMIN) /800. 

SSY=lYMAX-YMIN) /800. 

SM=S$X 

IF(SSY.GT.SSX)  SM=SSY 

IDUP=0 

LDUP=0 

23  READ (10, 20QO, END=44)  X(J),Y(J),Z(J),CD(J) 


DO  29  1=1.  - 
DIST=S0RTt(X[ 


H 


Y(J))**2)/SM 


.Lc.DSET)  THEN 


ir I  j 


.NE.CDjl) 


HIM 


=J+1 
0  TO  23 


PRINT  4455, IDUP+l.MSET 
F0RMAIIT12J.5H.PPINT5 


PRINT*. lELEVATIONS 
PRINT*.  POSITION  0 
PRINT,  A\Jy  ... 
READ  !Al  .TRASH 
PRINT*.  ’ 

fnLDUP.GT.O)  THEN 
PRINT  V  ’ 


OINTS  WITHIN  .I4.27H  SCREEN  UNI 


7H  SCREEN  UNITS  OF  EACH  OTHER) 

iMBP  o«  WW* 

RESS  RETURN  TO  CONTINUE* 


C6 


?PINT  4456.LQUP+1 

44  5^  FORMAT : |l2. 17H  DUPLICATE  LABEL?) 

PRINT*.  INFO  FOR  LAST  DUPLICATE  LABEL  IN  PILE  UAS  USED.' 
PRINT*  S 
PRINT. 


-  ,,IA\  V  .. 
READ  ’(Al!  .TRASH 
PRINT*.'  ’ 

ENDIF 

Xf1AX=-9<? .  E30 
Xf1IN=-XHAX 
Yf1AX=Xf1AX 
YfllN-XHIN 
-,33|:p 


.PRESS  RETURN  TO  CONTINUE’ 


ffl 

IF 

II 


ill 


XHAX 
.  YHAX 
.  XMIN 
.YfllN 
.  ZMAX 
.zniN 


33  CONTINUE 
SUBX=0. 

SUBY=0 

IFIABS  XHIN  .GT.lQO 


Xf1AX=X 

YMAX=Y 

XM1N-X 

YHIN=Y 

2MAX=Z 

2hIN=2 


*  ( XPIAX-XMIN) 


+  ABSC YHINl .67. 100. *  YHAX-YHIN 
L0^=L0GipiXriI(4  j. 


,0R. 

I  THEN 


IF  Xf1IN.Lt.*,. 
XHANTIS=XH!n/L. 
IM=XHANTIS*1000 

xn=in, 

XMIN=(Xf1 
L06-LQG1 


=LOG- 1 
LOG) 


*(10. “LOG) 


11/1000)  * 
...lOlYHIN 
IFtYfllN.LT.  1,1  L0G=L0G-1 
YflANT IS=Yf1IN?  { 10 .  *  LOG) 
If1=YHANTIS*100O 


YM=Ih 
Yf1IN=(Yf1 
SUBX=216 


8B8B  DATA  POINTS;  ,14) 


T  (A\J 

^R^LTIPLIER  FOR  Z  (RETURN  IF  HAPPY) 


ZflULT 


THEN 


923 


fRINT 

+REr 

^A^pC(F20.0)’] 

swh 

GO  TO  $20 

B8Ei. 

«■ 

PRINT  81. . ,  — ... — . . , 

+  -flIN(XHIN,  XHAlt)  ,Hlfl(YHAX.YNIN)  ,HAX(YHAX,  YHIN 
+  MAX  f  YMAX,  VHIN) -flINf  YHAX,  VniN) 


FACY=1 .  , 

820  PRINT*. ’  ’  ,  ,  , 

PRINT, 8889, HIN ( XHIN, XHAX ), HAX(XHIN, XHAX ),HAX( XHIN. XHAX) 

+  -MIN ( XMIN, XHAX] .HINiYHAX.YHIN) ,HAx(YHAX, YHIN! , 

+  HAX(YHAX,$HIN)-fl!N(  YHAX,$HIN) 

8889  F0RHATI6H  XHIN=, F15. 5, 6H  XMAX= , F 1 5 . 5, 6H  QIFF=, F15. 5./, 

+  ,j,H  ?PIN=  F15.5.6H  YHAX=  F15.5.6H  DIFF=  F15.5 

•'S'fSlo  (lt|RHHj|LTIPLIERS  FOR  X  AND  Y  (RETURN  IF  HAPPY) 

iRzCHCNE  „  ’)  THEN 

mwwaa. 


321 


822 


XHULT.YHULT 


READ143  *.ERR=821) 

print* ,  Inter  two  nuhbers  separated  by  a  cohha.’ 

CLOSE ( i3 ! 

GO  TO  820 


CONTINUE 

CLQSEU3. STATUS; 
F  acx=Facx*xmult 

FACY=FACY  YHULT 
DO  824  1  =  1  .  J 


DELETE’: 


C7 


kit:;*!: i -vnm  T 
?24  Y f 1 1 =YfT 1  * YnOr I 
xn  i  n = v  n  i  n  *  x  fiUL  t 
Yf1IN=YrtIN’V(1ULT 
xnAX=>:(iAx*xnuLT 
YHAX=Yf1AX*YHULT 
60  TO  820 
ENDIF 

Xf1IN=<59.  £30 
YniN=xniN 

mm* 

00  553  1=1. J 

IF  X  .LT.XMIN  XMIN=X  I 

IF  Y  .LT.YfllN  YflIN=Y  I 

IF  X  . GT.XNAX  XMAX=X  I 

IF  Y  .GT.YHAX  YHAX=Y  I 

553  CONTINUE 


WRITE  t  Si ’ [2A1]’ 


HONOR. SLASH 
ICBOR, SLASH 


WRITE  8,:  PQ„ 

WRITE  8,  A4J  )  TYP 
WRITE  8. *  FACX,FACY 
WR^T|  8  *)  SUBX.SUBY 

Sll  !lT|‘«!i:W5F,,Msl'“) 

WRITE  11,  *1  SUBX.SUBY  ,  . 

IFIASIS.EQ.  Y\0$.ASIS.EQ.  V  )  THEN 
DO  87  1=1. J 

|Si|li[¥l(W(i),2(I),I+1000,’  ’  ,LLBL 


EO, ' TR I \ ’ )  P0=0. 


iS  i  !t  ! 

m.15S 


MhwOe 


■  OR.ISe.EQ.  V] 


FUBQR ' EQ  2  J  IGG-I66-2 

RINr  ’(i;j6)  WANT  NORMAL ( N)  OR  FINE(F) , INVERSE  POWER’  , 
RINT,  (ill  ’  IASK  0IS™N(-E  GRIDOING  (0EF=N)  ?  (F  OR  N) 

6tiut»  "  '  HUT  > 


...WAIT... 

IGRI0=IGG 


PRINT*.’...,.,,.. .WAIT...’ 

IFlIASK.EO.  F  )  IGRI0=IGG 
ENDIF 

XGR0=IGRI0 

DHAX=(TXf1AX-XMIN}**2+(  YHAX- YfllN)  *  *2 )  *  *0. 5 
HLAB=1 

00,30. 12=1, IGRID+l  +  IBOR 
K= I  - 1 ) *  *12 

IFTK.lt. 0]  K=0  ... 


xb=Tk-q.  5*  1 1  xtiAX-xniNT /xgro)  )  ♦xiiim*  ( 1 1-1- 


cm 

ST=0. 


i  4w!£PP!,Iii’fl8'FACYI' 

F  6.Lt.0.007$‘Df1AX)  GO  tO  36 


IBOR/2 ) * ( ( XMAX-XHIN ) ) 


*2))**. 5 


lF?H0N0R.E0.’Y’.0R. 

IF(0AT.0.0075’Df1AX 


HONOR. EO.’ 
I  GO  TO  3C 


ifiMii,rtfristKE-,'!T<E" 

WRITE! 11 . * j  XB.YB, ST/SB, HLAB, ’ 

GO  TO  28 

5f|  l,EO. l.ANO.K.EO.Q)  THEN 
WRI  El  11. M  XB. Y&, St/SB, MLAB. 

GO  TO  28 


C8 


ENDIF 

IF! 12, EQ. IRRIO+l+IBORI  THEN 

up 7 TP i i  <  ' I  ya  va  cT/ca  m  »a  ' 


WRITE ’ ! I . * T  XB.  YB.ST/58.hLA8.  ‘ 

GO  TO  28 

Ilii!?;if!8:}S!W:SrSP.“»  »!  ’"E“ 

enotf  28 

28  £0N|JNUE 

IF?  SG.NE  ’Y’.AND.ISG.NE.  V),G0  TO  31 
Kpk  XB.YB,  ST/SB.  fILAB, 

31  CONTINUE 
30  CONTINUE  ,  , 

IF (HONOR . EQ . ’ Y ’ . OR . HONOR . EQ . ‘ y ’ )  THEN 
DO  55  1=1 . j 

Md  ,1,  fL.BL=* 

URITEhl.M  *(lT,Y(I).Z(I). HL AB+IOOO .  ’  ’.LLBL 
I1LAB=HLAB+1 
55  CONTINUE 

158  fe]1  , 

OPEN(22,nLEp  IJCPICK.DAT’) 

|  ;i  I  I 

15S  railffil#” 

l=[YB/FACYltSUBY 
!ITE(22, 160)  XB,YB.ST,DUf1f1 


157  CONT 
160  FORK 


mw 


;16.0,F16.2,A24) 


!?r’:...N0U  GOING  TO  CONTOURING  ACTIVITY.... 
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